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Construction of a Diffraction Grating Spectrograph 


R. L. PurBrIck 
Willamette University, Salem, Oregon 


(Received August 11, 1952) 


The construction and operation of a simple Wadsworth-mounted 21-ft 10-in. diffraction 
grating spectrograph is described. The instrument covers the range from 9000A to 1800A in 
the first order with a linear dispersion of 5A/mm. 


A pine for mounting is one of the problems 


which presents itself to a group planning the 
construction of a grating spectrograph with 
sufficient dispersion for both research and class 
use. After considerable investigation we decided 
to construct a Wadsworth-type mounting for a 
21-ft 10-in. diffraction grating. This mounting is 
compact, and a normal stigmatic spectrum is 
produced at the normal to the grating. Various 
workers have described their versions of this type 
of mounting.! Our modification yields a rela- 
tively simple instrument, both from the stand- 
point of construction and operation. 


DESCRIPTION 


The arrangement of optical elements for a 
Wadsworth-type mounting is shown in Fig. 1. 
Two equations express the conditions which must 
be satisfied : 


n\ =d(sini+ sinr), (1) 
s=R/(1+cos2), (2) 

where » is the order number, A the wavelength, 
d the grating space, 7 the angle of incidence, 


1W. F. Meggers and K. Burns, Bur. Standards Sci. 
Paper 411, 18, 185 (1922); R. F. Jarrell, J. Opt. Soc. Am. 
32, 666 (1942); J. W. Forrest and H. W. Straat, J. Opt. 
Soc. Am. 32, 669 (1942). ; 


r the angle of reflection, s the distance between 
grating and photographic plate, and R the radius 
of curvature of the grating. Provision must be 
made for changing 7 and s. 

Figure 2 shows our version of the mounting. 
The main support is formed by 6-in. channel 
irons welded together by }-in. steel plates. Five 
4-in. diameter steel pipes, 3 ft in length and with 
flanges at both ends, constitute the legs. Another 
6-in. channel iron has the grating at one end and 
the plateholder at the other. A Timken 366 
bearing mounted directly beneath the grating 
supports one end of this beam, and a similar 
bearing near the plateholder supports the other 
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Fic. 1. Arrangement 
of optical elements for 
a Wadsworth-type 
mounting. The arrow 
indicates the direction 
the plateholder must 
be moved for higher 
wavelengths. 
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end. A steel pipe extends from the slit towards the 
mirror. This prevents light from the slit from 
directly striking the photographic plate. 

The concave mirror is 6 in. in diameter and 
has a radius of curvature the same as that of the 
grating. It renders the light incident on the 
grating parallel and is the primary factor con- 
tributing to a stigmatic image. 

The grating was ruled by Dr. John Strong’s 
group at The Johns Hopkins University, It is 6 in. 
in diameter, has a 21-ft 10-in. radius of curvature, 
14 400 lines per inch, and high intensity in the 
second order. The grating is mounted in a 
dustproof cover, and provision is made for 
rotation about three mutually perpendicular axes 
to facilitate alignment. 

A round straight steel bar at the plate end of 
the instrument serves as a rail for moving this 
end of the channel iron on which the grating 
and plateholder are mounted. The rail is straight, 
but this end of the channel iron follows the arc 
of a circle when moved; hence it was necessary 
to mount it as shown in Fig. 3. Three motions 
are allowed: (1) movement along the rail, (2) 
rotation of the channel iron about its point of 
support, and (3) movement along a line connect- 
ing the grating and plateholder. This arrange- 
ment has one advantage over a steel bar ground 
to the arc of a circle; the channel iron here can 
be rotated for any position of the grating on its 
steel plate support. 

The in-and-out motion of the plateholder, 
which must be carried out to satisfy Eq. (2), is 
provided for as shown in Fig. 3. Set screws at 
the ends of the 2 in. X10 in. photographic plate 
are used to curve the plate for proper focus. 


Fic. 2. The spectrograph mounting described in this 
paper. The case for the instrument is not shown 


PURBRICK 


Unfortunately this curvature varies slightly with 
spectral region. Focusing can be visual with a 
ground glass plate, or a scale may be mounted on 
the channel iron and photographic alignment 
used to calibrate the scale. 

The mounting is encased in a cabinet formed 
from 3-in. plywood from the top of the legs 
upward. The case follows the outline of the 
mounting and is 2 ft 5 in. in height, 8 ft 5 in. 
wide, and 13 ft long. Two 1 ft X4 ft hinged 
doors are located at the plate end. 


PERFORMANCE 


The instrument is aligned with Eqs. (1) and 
(2) in mind. Figure 4 illustrates how critical the 
position of the plateholder is for proper focus. 


|i 
cs ae 


Fic. 3. The plateholder end-of the movable beam. 


The slit width was held constant; the mercury 
2536A line was photographed, and its width was 
measured for various positions of the plate- 
holder. The curve indicates that the plate should 
be positioned to millimeter accuracy. 

The slit and grating rulings must be parallel. 
This is best achieved photographically. Two 
neon lines at 6143A and 6164A are shown in 
Fig. 5; (A) was taken with a small angle between 
slit and rulings while in (B) they were parallel 
and the images much sharper. 

The spectrograph is mounted on a concrete 
floor laid on the ground, and to date no diffi- 
culties have been encountered as a result of 
vibration. During long exposures an effort is 
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made to keep the room temperature reasonably 
constant. However, this has posed no serious 
problem in exposures up to 12 hours. 
The instrument covers the region from 9000A 
to 1800A in the first order with a linear dis- 
persion of approximately 5A/mm. The disper- 
sion varies about 0.2A/mm between the red and 
ultraviolet. 
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Fic. 4. The variation of spectral line width as a function 
of plateholder position on the moving beam. Negative 
values indicate displacement from the point of sharpest 
focus toward the grating. 


The author is indebted to the Research (a) (b) 


Corporation for a Frederick Gardner Cottrell 
grant, which aided in this project, and tostudents 
Robert Bennett, James Carpenter, Lawrence 


Fic. 5. Photographs of two neon lines showing the 
necessity of having slit and grating rulings parallel. (a) 
was taken with a small angle between slit and grating 
rulings; in (b) they are parallel. 


Cherry, Alan Goudy, Frederick Gruver, and 
Philip Phipps, who helped in the construction 


Carpenter, Sr., contributed greatly to the con- 
as part of a special problems course. James 


struction of the welded portion of the instrument. 


The chief objects of the [Royal] Institution were to provide a speedy and general diffusion of 
knowledge of all new and useful improvements, to teach the application of scientific discoveries to 
the improvement of arts and manufactures in England, and increase domestic comfort and con- 
venience. Under Rumford’s guiding hand a house was purchased in Albemarle Street and laid 
out according to his plans. There was a lecture theatre, a model room to exhibit mechanical inven- 
tions and improvements, workshops and kitchens. A school for mechanics was established to extend 
the benefits of science to the working poor. Professors and lecturers were appointed to give the 
nobility and wealthy an opportunity to share in the advantages of scientific instruction. Various 
committees were established to supervise scientific investigation; a printing press was bought, and 
a journal was to be published. The Institution was to be a ‘‘clearing-house’’ of scientific information, 
a source from which great good could flow. As Rumford said: ‘‘I am only desirous that science and 
art should once be brought cordially to embrace each other, and to direct their united efforts to the 
improvement of agriculture, manufactures, and commerce, and to the increase of domestic comfort." 
—GeorceE A. Foote, Isis 43, 7 (1952). 





On a Routine Analytic Method for the Solution of Problems in Statics 


A. W. Srmon 
Alabama Polytechnic Institute, Auburn, Alabama 


(Received September 15, 1952) 


A routine analytic method for the solution of problems in statics is given. The method re- 
quires no geometry other than that required for the identification of the coordinates of certain 
cardinal points involved in the problem. The method is illustrated by the detailed solution of 


an actual case. 


HE methods of solution of problems in 
statics given in most textbooks are charac- 
terized by the fact that they depend largely on 
the use of special tricks and artifices, which are 
applicable in one case but not in another, with 
the result that the student is usually at a loss 
how to proceed and often wastes a great deal of 
time before he does find the proper method of 
attack. Furthermore, the methods of solution 
given usually make the mechanics subservient to 
the mathematics, rather than vice versa. In this 
paper a method is presented, which is straight- 
forward and directly applicable to any problem 
in statics, and in which the equations set up on 
the basis of the mechanics involved serve to indi- 
cate the best method of mathematical solution. 
Furthermore, the method to be presented re- 
quires no geometry (or trigonometry) beyond the 
ability to identify the coordinates of points, the 
entire solution being analytic, that is algebraic, 
after the coordinates of certain cardinal points 
involved in the problem have been determined. 


I. THEORETICAL BASIS 


A method for the solution of problems in 
statics, and one particularly convenient in the 
case of forces in space, can be based on the ana- 
lytic forms of the fundamental equations 


> X=0, (1) <=M.=0, (4) 
LUY=0, (2) 2M,=0, (5) 
~Z=0, (3) ~M,=0, (6) 


in which, as usual, X, Y, Z represent the axial 
components of the force F, and M., M,, M,, 
represent the axial components of the moment, 
with respect to the origin, of the force F. 


In analytic form, Eqs. (1), (2), and (3) become 
> F(x2—x1)/r=0, (7) 
LU F(y2—-91)/r =0, (8) 
LF (22—21)/r=0, (9) 


where x1, 1, 21, and x2, V2, 22, represent the co- 
ordinates of an (arbitrarily selected) first point 
P, and second point P2, respectively, on the line 
of action of the force F, and r represents the dis- 
tance between the two points. The distance r is, 
as usual, given by 


r? = (x2 —X1)?+ (v2 —y1)? + (22 —21)?. (10) 


Furthermore, in analytic form, Eqs. (4), (5), 
and (6) become 


XU(yZ —zY) =0, (11) 
> (2X —xZ) =0, (12) 
L(«xY —yX) =0, (13) 


in which x, y, 2, on the right-hand side now 
represent the coordinates of a point on the line 
of action of the force F. For convenience this 
point may be taken to coincide with either one 
of the points 1 and 2, the coordinates of which 
have previously been determined. 

In the following we shall write the coordinates 
of the points P; and Py» in the form P(x, 1, 21) 
and P2(x2, ye, 22), and the vector P,P, extending 
from P; to P: as 


P,\P.2=[(x2—x1), (ye—y1), (22 —21) J. 


While any point may be chosen as the origin 
of coordinates, the moment equations can ordi- 
narily be simplified if a point of concurrency of 
two or more forces is selected as the origin, pro- 
vided, of course, such a point exists, as is usually 
the case. Furthermore, if part of the data is in 
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trigonometric form, a convenient artifice is to 
choose two points which are unit distance apart. 
The direction cosines of the corresponding angles 
may then be used directly. 

With regard to the direction of the unknown 
forces, an alternative of the above scheme is to 
assume all unknown forces to be in such direc- 
tions as to correspond to tensions; a positive 
result then indicates tension, a negative result a 
compression. 


Il. PRACTICAL PROCEDURE 


(1) Identify the coordinates of a first and a 
second point on the line of action of each force 
involved, both known and unknown. 

(2) Compute in order for each force the differ- 
ences (x2—%1), (ye—1), and (z2—2;), the magni- 
tude of r, and the ratios (x2—x1)/r, (ye—91)/r, 
and (z2—2;)/r; that is, the direction cosines of 
the vector P,P2, or, alternatively, the compon- 
ents of a unit vector in the direction of P,P». 

(3) Insert the values so found in Eqs. (7), (8), 
and (9). 

(4) Insert the same direction cosines and the 
coordinates of either one of the points P; and P: 
(or, if desired, other convenient point on the line 
of action of the force) in Eqs. (11), (12), and (13). 

Solve the resulting six equations simultane- 
ously to determine the magnitude of the unknown 
forces. A positive result means the corresponding 
force is actually in the direction from P, to P2, 
while a negative result means it is opposite 
thereto. 


Ill. ILLUSTRATIVE EXAMPLE 


The method can best be illustrated by working 
out in detail an actual case. 

For this purpose consider the structure pic- 
tured in Fig. 1, in which a horizontal boom is 
supported by a ball-and-socket joint, or equiva- 
lent, at one end and by two cables at the other 
end, while a force of 100 lb is applied in the 
direction shown. Let the dimensioning be that of 
the figure. It is desired to determine the force 
which each of the cables exerts on the boom, as 
well as the three components of the reaction at 
the ball-and-socket joint. If the various cardinal 
points are lettered as in Fig. 1, the desired forces 
are those along (but not necessarily in the same 
direction as) the vectors AB, AC, AD, DE, and 
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Fic. 1. System of forces analyzed in the illustrative 
example. The directions of the forces indicated are arbi- 
trarily assumed. 


DP, with the known force in the direction of the 
vector GH. 

If now the origin is taken at A and a right- 
handed system of coordinates with the x axis 
extending horizontally to the right and the y axis 
vertically upward is employed, the coordinates 
of the cardinal points involved are 


A(0,0,0),  B(—12,10,4), C(—12,3, —6), 
D(-12,0,0), E(—12, —5,0), F(—12, 0, 4), 
G(-2,0,0), H(—4, —5, 6). 


It is to be noted that further reference to the 
figure is now unnecessary, the rest of the compu- 
tation following automatically. 

For each of the vectors involved we now write 
the two end points in inverted order and subtract 
corresponding coordinates. Thus, for AB, we 
write 

B(—12, 10, 4) 
A( 0, 0,0) 


AB=[-—12, 10, 4]. 


The corresponding value of r is 


= 4/260 =16.1, 


while the corresponding direction cosines or com- 
ponents of a unit vector in the direction of 
AB are 
—12/16.1= —0.744, 10/16.1=0.620, 
4/16.1=0.248. 
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The results may be conveniently written in 
the form 


AB,=[ —0.744, 0.620, 0.248], 


where the subscript 1 indicates a unit vector. 

If these computations are carried out for each 
of the vectors involved, and the results gathered 
in a table, we have 


AB,=[ —0.744, 0.620, 0.248], 
AC,=[ —0.873, 0.218, —0.436], 
AD,=[-1, 0, 0], 

DE,=[0, —1, 0], 

DPF,= [0, 0, 1], 

GH,=[ —0.248, —0.620, 0.744]. 


In the case of such forces as AD, DE, and DF, 
which are in directions parallel to the coordinate 
axes, the corresponding direction cosines can, of 
course, be determined by inspection without 
recourse to calculation. 

Substitution in Eqs. (7) to (9), inclusive, yields 
(the two letter symbols now represent the magni- 
tudes of the corresponding forces) 


—0.744A B —0.873AC—AD 
—0.248 X 100 =0, 

0.620AB+0.218AC —DE 
—0.620 X 100 =0, 


0.248A B —0.436AC+DF 
+0.744 X 100 =0. 


(7a) 
(8a) 


(9a) 


We turn next to the moment equations. Select- 
ing first Eq. (11), we note that if we choose the 
points A, D, and G, for each of which y=z=0, 
the corresponding equation reduces to the form 
0=0. This is merely the mathematical expression 
of the fact, obvious from the figure, that none of 
the forces involved have a moment around the 
chosen x axis. 

Turning next to Eq. (12), we note that if for 
the forces AB, AC, and AD we select the point A, 
for which z=x=0, the corresponding terms 


SIMON 


vanish; while, if for the force DF we select the 
point D for which z=0, the only term left is 
that corresponding to -—xZ, which becomes 
—(-—12)DF; further, if for the force DE we 
select again the point D, the first term vanishes 
for the same reason as before, while the second 
term vanishes because for this force Z =0; hence 
DE makes no contribution. This leaves finally 
the known force GH, for which we select the 
point G, and have as its only contribution 
the term corresponding to -—xZ, which is 
—(—2)(0.744)100. Hence, Eq. (12) reduces to 


~(-—12)(1)DF —(—2)(0.744)100=0. (12a) 


Similarly, by a judicious choice of points in- 
volved in Eq. (13), the latter reduces to 


~12(—1)DE+(—2)(—0.620)100=0. (13a) 


A glance at Eqs. (7a) to (13a), inclusive, now 
indicates the method to be followed in complet- 
ing the algebraic part of the solution. 

The solution is DE= —10.3, DF= —12.4, 
AB=27.9, AC=158, AD= —184. The minus 
signs indicate that the corresponding forces are 
opposite to the directions originally assumed. 

It will be noted that an alternative choice of 
origin is the point D, in which case (axes parallel 
to the original directions) the coordinates of the 
points A and G become 


A(i2,0,0) and G(10,0,0), 


and the corresponding moment equations re- 
duce to 


—12(0.248)A B —12( —0.436)AC 
~10X 100 X0.744=0, 
+12(0.620)AB+12(0.218)AC 
+10 100( —0.620) =0. 


The direction cosines are, of course, unchanged 
by a shift of origin provided the axes of the new 
coordinate system are taken parallel to and in 
the same direction as the original. 





The Constant Battery Wheatstone Bridge 


T. TOWNSEND SMITH 
University of Nebraska, Lincoln, Nebraska 
(Received July 31, 1952) 


The article discusses the sensitivity of a constant battery Wheatstone bridge and shows that 
under certain conditions the bridge sensitivity will be independent of the magnitude of the 
resistance to be measured. A chart is given, which indicates for a given resistance the bridge 


ratio which provides the highest sensitivity. 


The bridge is insensitive for measuring resistances which are large in comparison with the 


resistance of the ratio arms. 


HE question of the sensitivity of a con- 

stant battery Wheatstone bridge is inter- 
esting and of some importance in an electrical 
laboratory. 

Suppose that a battery of emf E and of 
negligible resistance is connected to a Wheat- 
stone bridge as indicated in Fig. 1. Suppose 
further that the bridge as shown in Fig. 1 is in 
balance and that the resistance R; is then in- 
creased to R,+AR,. The sensitivity S of the 
bridge may be defined in terms of the galvanom- 
eter deflection D and the relative change in Ri, as 


(1) 


The sensitivity thus defined can be expressed in 
terms of E, of the resistances of the various arms 
of the bridge, and of the resistance G and the 
sensitivity o of the galvanometer, where o is, by 
definition, 


S=D(R/AR)). 


o=D/Ig. (2) 


As a first step one may write for the current 
through the galvanometer, in terms of the cur- 
rent through Ri, 


I,(AR,) 


¢=———_—____—_—__-. (3) 
G(1+R2/Ri)+Rit Re 

One may now express J,, to a first approximation, 
neglecting terms in Ig/JI,, since Ig¢/I; will be 
small for a bridge only a little out of balance. To 
this approximation the value of I; may be 
written 

I,=E/(Rit+R3). (4) 


Equations (1), (2), (3), and (4) may be combined 
to give an expression for the sensitivity of the 


bridge. Thus, 


Ri 
S=D(R,/AR;) = Eo 


it Rs; 
1 
G(1+R2/Riy)+RitRe 


One may express S as a function of a parameter 
(Ri+R.)/R, where R is used as a notation for 
R3+ Ry, with the aid of the relations 


Ri/Rs= R2/Ra= (Rit R2)/R. 
So expressed, Eq. (5) becomes 
(Ri+R2)/R 
"T+ (Ri+R:)/R 


(5) 


(6) 


1 
“GE1+ (Rit Rs)/RI+R(Rit R:)/R 


If p is used as a symbol for (Ri+R:2)/R, Eq. (7) 
becomes 


(7) 


p 1 
ere 
(1+p) Giit+tp)+Rp , 


Equations (7) or (8) show that the sensitivity 
of the bridge does not depend upon R; directly 
but only on the ratio (Ri+R:)/R. As long as this 
ratio is kept constant, the sensitivity of the 


(8) 


Fic. 1. A “‘center- 
connected” Wheat- 
stone bridge, with 
the battery con- 
nected between the 
junction of the ratio 
arms (R;—R,) and 
the junction of the 
known and unknown 
resistances (R,— Re). 
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Fic. 2. The sensitivity of a constant battery, center- 
connected Wheatstone bridge. This chart may be used to 
obtain the resistance ratio which should be used to get high 
sensitivity in measuring a particular resistance. The ordi- 
nates are S/Smax. The abscissas are marked in terms of 
R,/R, and also in terms of Ri, for the case R= 1000 ohms. 


bridge will be unchanged and any resistance may 
be measured with the same precision, in so far as 
the precision depends on the sensitivity of the 
bridge. 

The preceding paragraph should not be inter- 
preted as asserting that with this type of bridge 
all resistances can be measured with the same 
precision, for that is not true, though it is true 
that all resistances of less than a certain critical 
value (which depends on R and G) may be meas- 
ured with the same precision. If Ri is greater 
than this critical value, the sensitivity is less. If 
R,/R is large, the bridge is very insensitive. 

The bridge sensitivity is a function of p and G. 
It is a maximum when p is $ and when G is 3R. 
To find the value of p for which the sensitivity is 
a maximum, G remaining constant, it is conveni- 
ent to write Eq. (8) in the form 


Ea 
Se ann, 
[G(1+p)*/p]+R(1i+p) 


and to find the value of p for which the denomina- 
tor is stationary. This value is given by 


p°=G/(R+G). (10) 


To find the value of G which makes S a maxi- 
mum, one needs first to express the galvanometer 
sensitivity o in terms of the galvanometer re- 
sistance G. For a given type of galvanometer this 
relation is expressed to a fair approximation by 
the equation! 


(9) 


o=0G'. (11) 


1T,. T. Smith, Electricity for Students of Science and 
Engineering (International Textbook Company, Scranton, 
1949), p. 128. 


By combining Eqs. (9), (10), and (11), one gets 
an expression for S: 
Gi 
= Eop—__—_—_. 
[Gi+ (G+R)*} 
The value of G for which S is a maximum is found 
by setting equal to zero the derivative of the 


expression on the right side of Eq. (12). This 
procedure gives 


(12) 


G=1R. (13) 


If now this value for G is inserted in Eq. (10), the 
result is 


e=(RitR2)/R=}. (14) 


Thus, as indicated above, the bridge is most 
sensitive when p is $ and when G is 3R. 

Since R; and R: are both positive, it follows 
from Eq. (14) that the largest resistance which 
can be measured with the maximum sensitivity 
is 3R. An.expression for this maximum sensi- 
tivity can be obtained from Eq. (9) by substitut- 
ing 3R for G and 3 for p. The result is 


Smax = Eo/3R. (15) 


The curves of Fig. 2 show the ratio S/Smax as a 
function of the bridge ratio (r= R:i/R2=R3/R,) 
and of the ratio R,/R. The formula by which 
these curves are calculated is obtained by ex- 
pressing p in terms of Ri/R and 17, substituting 
this value in Eq. (9), and then combining the 
resulting equation with Eq. (15). The expression 
for p is 


Rit Re R, 1 Ry 1+r7 
Ry Set gg 
R R r R r 


and the resulting formula for calculating S/Smax is 


S 9 


Seen 1+7r\Ri r\R « 
(EE 
r JR r+1/ R, 
The curves are drawn for various values of r from 
0.001 to 10, with an indication of the position of 
the curve for the limiting bridge ratio of infinity. 
We have in our electrical laboratory a bridge 
of the type here discussed with a resistance R of 
1000 ohms and with a dry cell for E. The gal- 
vanometer has a sensitivity of 45 10° mm/am- 
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CONSTANT BATTERY WHEATSTONE BRIDGE 


pere and a resistance of 265 ohms. The peak 
sensitivity is about 2600, which is almost exactly 
the value which would be given by Eqs. (9) and 
(10), or by Eqs. (15) and (11). A chart showing 
the curves of Fig. 2 hangs by the bridge. As soon 
as the approximate value of an unknown .re- 
sistance is found, a glance at the chart enables 
one to choose the bridge ratio for the highest 
sensitivity. 

Some comments seem to be pertinent: 

(1) A change in the galvanometer resistance G 
shifts the curves of Fig. 2 without greatly altering 
their form. Take, as an example, a bridge ratio of 
0.1. With the best galvanometer (G=}R=333 
ohms for an R of 1000 ohms) the peak sensitivity 
is at 45 ohms. In our bridge G is 265 ohms, for 
which the peak sensitivity is at 41 ohms with a p 
value of 0.46 [Eq. (10) ]. When G=500 ohms the 
peak sensitivity would be at 52 ohms with a p 
value of 0.58, and for G=1000 ohms the peak 
would be at 71 ohms with a p value of 0.71. 

(2) A switch for interchanging galvanometer 
and battery would add little to the usefulness of 
the bridge. If the battery and galvanometer were 
interchanged so that the battery would be across 
the ends of the ratio arms, the sensitivity with a 
1-1 bridge ratio would be somewhat increased 
for small values of R:. With ratios of 0.1 or 10 
the change would decrease the sensitivity for 
small values of Ri. For large values of Ri both 
bridges are insensitive and there is little to choose 
between the two arrangements, though the end- 
connected bridge may be made a little more 
sensitive. For general use, considering both sensi- 
tivity and convenience, the center-connected 
bridge is to be preferred. 

(3) This note deals only with a Wheatstone 
bridge used as a testing set with a constant 
battery, perhaps a single dry cell. For resistance 
comparisons of high precision one would nor- 
mally increase the sensitivity by using an electro- 
motive force that provided currents close to the 
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tolerance limits of the bridge conductors. For 
such work the conditions for high sensitivity and 
for best galvanometer resistance are entirely 
different from those discussed here.? For example, 
for high precision measurements a desirable 
bridge arrangement is the equal arm bridge, 
Ri=R2=R;=R,, and in this case the best gal- 
vanometer resistance is Ri(G=R,). In such work 
bridge sensitivities of the order of a million 
may be obtained with fairly insensitive gal- 
vanometers, if R, is large enough to use a Wheat- 
stone instead of a Kelvin bridge. With precision 
measurements the sensitivity should be inde- 
pendent of Ri, if Ri is large enough for Eq. (11) 
to be applicable. With a constant battery the 
sensitivity is proportional to (1/R:)* for large 
values of Ri. 

(4) A bridge of the sort here discussed is not 
suitable for large bridge ratios. For example, 
assume that R is 1000 ohms, R:/Rz is 0.001, R: is 
about 1 ohm, and a dry cell is applied directly to 
the junctions Ri;—R2 and R;—R,. With these 
values the power dissipation in R; and R; will 
be somewhat over a watt each. This power may 
be above the tolerance limits of the coils. 

(5) For a galvanometer of a given type the 
sensitivity of the bridge changes very little for 
relatively large changes in G. For the center- 
connected bridge the best value of G is 3R. If G 
is changed by a factor of 3, to R or to $R, Eq. 
(12) shows that S changes by a factor of 0.9. 
Galvanometer resistance G would need to be 
changed by a factor of nearly 10 to reduce the 


sensitivity to 3 of the peak value. 


2T. T. Smith, Electricity for Students of Science and 
Engineering (International Textbook Company, Scranton, 
1949), pp. 54-56; Glazebrook’s Dictionary of Applied 
Physics (Macmillan Company, New York, 1922), Vol. 2. 

p. 713, 714. Other approaches to this question are given 
in Page and Adams, Principles of Electricity (D. Van Nos- 
trand Company, Inc., New York, 1931), pp. 175; Harnwell, 
Principles of Electricity and Electromagnetism (McGraw- 
Hill Book Company, Inc., New York, 1949), second edi- 
tion, pp. 133, 134. 
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Spin Waves in Ferromagnetic and Antiferromagnetic Materials* 
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A brief outline is given of the spin-wave approximation to the near ground states of ferro- 
magnetism and antiferromagnetism. Simple pictorial models of spin waves are introduced. 
These models clarify the striking difference between ferromagnetic spin waves, which obey the 
dispersion law w « k*, and antiferromagnetic spin waves, which obey the law w «<k. The proper 
normal modes for antiferromagnetic spin waves are presented. The nature of ferromagnetic and 
antiferromagnetic resonance is illuminated from a spin-wave point of view. 


I. INTRODUCTION 


UR objective in this paper is to give a 

simple pictorial representation which dis- 
plays the striking difference between ferromag- 
netic and antiferromagnetic spin waves. Since 
the theory of spin waves is not too widely known, 
however, we shall begin with a few words about 
this approximation to the nearly ordered, low 
temperature states of ferromagnetism and anti- 
ferromagnetism. 


A. Some Remarks about Ferromagnetic 
Spin Waves 


At each lattice site in a ferromagnet let us 
assume a localized spin with spin quantum num- 
ber S and with z component of angular mo- 
mentum 4S. In the ground state of the spin 
system, which will be realized at 0°K, each spin 
has the maximum allowable value of S, namely, 
S. A spin wave may now be described as a sinu- 
soidal disturbance of the spin system with 
amplitude at each lattice site proportional to 
S-S?. 

It was first shown by Bloch? that the states 
near the ground state of a ferromagnet can be 
approximated by superpositions of these sinu- 
soidal spin waves. We shall indicate the direction 


* This research was assisted in part by the U. S. Office 
of Naval Research. 

t Now at Department of Physics, University of Pitts- 
burgh, Pittsburgh, Pennsylvania. 

t Now with Solid State Group, Massachusetts Institute 
of Technology, Cambridge, Massachusetts. 

§ Now at Department of Physics, University of Illinois, 
Urbana, Illinois. 
1For an excellent introductory review of the quantum 
theory of ferromagnetism the reader is referred to J. H. 
Van Vleck, Revs. Modern Phys. 17, 27 (1945); for an 
excellent review of ae to J. H. Van Vleck, 
J. phys. et radium 12, 262 (19 Yi 

F. Bloch, Z. Physik 61, 206 (1930). 


of Bloch’s quantum-mechanical proof and then 
go on to a semiclassical approach, introduced by 
Heller and Kramers,* which allows simple pic- 
torial interpretations of spin waves. 


Quantum- Mechanical Approach 


The exchange Hamiltonian of a system of 
spins is 


H= -—2J Les S;-S;+ const, (1) 


where the summation is taken over all combina- 
tions of nearest neighbors; J is the exchange 
integral between nearest neighbors. We shall 
limit our discussion to spin S=4, and shall intro- 
duce the Pauli matrices ¢o+=0,*+10,"; 0%. The 
eigenfunctions of o,* are a; (eigenvalue +1) and 
B; (eigenvalue —1). We assume either a small 
applied field or some anisotropy in order that z 
be the preferred axis; but for simplicity we will 
omit such terms from the calculation. The 
Hamiltonian (1) becomes (apart from the 
constant) 


HR’ = —3J Duy (Goitoj-—+ 40 -0;++0%07), (1’) 


which has the ground-state eigenfunction, if N 
spins are present, 


Yo = aa2013° + - ay. (2) 


This is the state at 0°K; it represents the maxi- 
mum alignment of spins, or saturation mag- 
netization. As the temperature increases the 
system will be excited out of the ground state. 
The next state may be thought of as one in which 


*G. Heller and H. A. Kramers, Proc. Roy. Acad. Sci. 
(Amsterdam) 37, 378 (1934). For a simplified reformulation 
of the semiclassical theory see C. Herring and C. Kittel, 
Phys. Rev. 81, 869 (1951), especially Appendix A. 
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a single spin is reversed; however, the function 
Pj = A1Q2°° * 0318 j00j41° *"Qn, (3) 


is not an eigenfunction of Eq. (1’). An eigen- 
function ¥, may be formed by taking a linear 
combination of ¢;, each member of the com- 
bination containing a reversal at a different 
lattice site: 

We = Dos C# 93. (4) 
It can be shown that each such proper combina- 
tion is equivalent to a spin wave of wave number 
k, and that the allowed values of & can be deter- 
mined from periodic boundary conditions.‘ 

We shall demonstrate the wave-like properties 
of the solutions (4) for the case of a linear chain. 
Requiring periodic boundary conditions is equiva- 
lent to bending the chain around into a ring so 
that the first spin is also the N+1 spin. Each 
spin 7 has two nearest neighbors, j7+1 and j—1, 
with which it interacts through the exchange 
integral J. The Schrédinger equation, 


RY. = EX, 

allows us to equate the inner product 
(5 | | Yu) = —ILe jap cj st +3(N —4) cj") 
with the inner product 
(93 | E| vx) = Ec;*. 
Therefore, dropping primes, 
LE+3(N —4) J]cP+JI(cj41*+0;-1*) =0. (5) 

This equation has solutions of the form 

c# = N-te*is, (6) 
where @ is the distance between spins. The 


periodic boundary conditions require c# =c;,y*, 


and hence, the allowed values of k are given by 
kaN/2x=0, +1, +2, +--+, &3(N—1), +3N. 


The wave-like property of the solutions (4) is 
shown by Eq. (6). On inserting Eq. (6) into 
Eq. (5), we obtain an equation relating the 
energy E and the wave number k: 


E+}3(N —4)J+2J coska=0. 
For small values of & this reduces to 


E-E\=Jk*a’, 


‘ This demonstration, which was basic to Bloch’s theory, 
was first given by J. C. Slater, Phys. Rev. 35, 509 (1930). 
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where Ey= —3NJ is the energy of the ground 
state. It is seen that E—Ep is the energy of the 
spin wave, E;. If we quantize this energy, we 
arrive at the dispersion law 


E,.=ho.&(const)Ja*k? (ferromagnet). (7) 


For a linear chain with S=} the value of the 
constant is unity; for a three-dimensional lattice 
of cubic symmetry and with arbitrary S the 
constant is ZS/3, where Z is the number of 
nearest neighbors. The extension of Bloch’s 
theory to arbitrary S was made by Meller.® 


Semiclassical Approach 


It is well known that the classical vector model 
of a precessing spin gives results in agreement 
with quantum mechanics if the length of the 
vector, in units h, is taken to be [.S(S+1) ]# and 
the values of S,; limited to —S, —S+1, ---, 
S—1, S. We may think of our lattice as composed 
of N such spins, each spin precessing about a 
small applied field or some effective anisotropy 
field determining the axis of precession. The 
ground-state function (2) is equivalent to all 
spins precessing in phase, each with S#=S. In 
the function (3) one of the spins has $= S-—1. 
What is the equivalent classical picture of the 
functions (4)? Heller and Kramers* showed that 
here the spins divide up the unit loss of z com- 
ponent, each S# being equal to S—(1/N). The 
reason this classical picture works is that the 
quantum-mechanical expectation value of S¥ in 
the states (4) is just S—(1/N). 

For large quantum number S we may take the 
classical picture quite literally. The many differ- 
ent functions (4) are then seen to correspond to 
the possible different phase relationships among 
the precessing spins. If in a given direction in the 
lattice each successive spin differs in phase from 
the preceding one by a certain amount, @, then 
as the spins precess a continuous wave will be 
traveling through the lattice. The wave ampli- 
tude might be measured, for example, by the 
value of S;*(#). This is a slightly different version 
of a spin wave than is Bloch’s propagated spin- 
reversal, but it gives nearly equivalent results. 

It is seen that as the phase angle @ gets larger 


(increasing k), more energy is required to estab- 


5C, Moller, Z. Physik 82, 559 (1933). 
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lish the spin wave; this is because neighbor spins 
are more out of alignment and, hence, more work 
must be done against exchange forces. Thus, the 
dispersion law (7) arises. A more detailed picture 
of a spin wave in a ferromagnet and a semi- 
classical derivation of Eq. (7) will be given in 
Sec. II. We shall see that the semiclassical ap- 
proach gives the same results as the quantum- 
mechanical in the usual limit S(S+1) ~ S?. 


The Principle of Superposition 


As the temperature is increased, more and 
more spin waves are excited. The basis of Bloch’s 
approximation is that the higher eigenstates of 
Eq. (1’)—that is, states in which several re- 
versals of spins are present—are very nearly 
combinations of functions (4), at least as long as 
the number of reversed spins is small compared 
to N. This is the principle of superposition, 
fundamental to all wave theories. It breaks down 
at a temperature of the order of one-tenth the 
Curie temperature. Beyond this temperature the 
spin-wave picture is of little merit. 

The superposition principle allows any number 
of spin waves of a given k to be excited. This 
means that spin waves, in common with photons 
and phonons, obey Bose-Einstein statistics. 
Thus, the number of spin waves excited at a 
temperature T will be 


n= Die m= Dx Lexp(Ex/kaT) -1}", 


where kz is the Boltzmann constant. To evaluate 
n we insert the dispersion law (7) into Eq. (8) and 
perform the summation over the possible values 
of k allowed by the periodic boundary conditions. 
In this way it is found that m, and hence the 
negative change of total S*, are proportional to 
T'. Or, if M is the magnetization, 


AM(T) = M(0°K) —M(T)=(const)T3. (9) 


This is Bloch’s T? law for the low temperature 
dependence of magnetization of a ferromagnet. 
Experimental measurements by Fallot® verify 
this law and may be said to confirm the spin- 
wave approximation. Further details, including 
the value of the constant in Eq. (9) for a number 
of cases, may be found in Herring and Kittel.’ 


®M. Fallot, Ann. phys. 6, 305 (1936). 
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B. The Case of Antiferromagnetism 


The exchange integral J is positive for ferro- 
magnetism, negative for antiferromagnetism. 
This means that nearest neighbor spins in an 
antiferromagnet will tend to align themselves 
antiparallel to one another. The simplest version 
of this is the two-sublattice picture, in which at 
0°K all the spins of one sublattice point up, all 
those of the other point down. For example, in a 
body-centered cubic all the corner sites might 
have spins up, all the body-center sites have spins 
down. This is not quite correct since the function 


Vo = 0182038405 - «By, (10) 


is not an eigenfunction of the Hamiltonian (1’). 
However, Anderson’ has proved that the exact 
ground-state energy eigenvalue must be close to 
the energy of the approximation (10). So we 
shall assume the two-sublattice picture to be 
approximately correct. 

One can now imagine, on the semiclassical 
approach, a spin wave being excited in an anti- 
ferromagnet. It would appear at first glance that 
such a wave would be very similar to its ferro- 
magnetic analog. Neighbor spins would be point- 
ing antiparallel, to be sure, and the total mag- 
netization would always be zero in the absence 
of an applied field ; but the precessions and phase 
relationships should be as in ferromagnetism. 
Thus, one might expect the dispersion law given 
by Eq. (7) also to hold for antiferromagnetism. 
That this is not the case was first shown by 
Hulthén,® who found 

E,=hwo,.=(const) k (antiferromagnet). 


(11) 


Recent work by Anderson® extending the spin- 
wave theory of antiferromagnetism to include 
the zero-point spin-wave energy has incidentally 
confirmed Hulthén’s dispersion law and has re- 
newed interest in the difference between ferro- 
magnetic and antiferromagnetic spin waves. The 
question arises: Why should it take so much 
more energy (for small k) to excite an anti- 
ferromagnetic spin wave than to excite a ferro- 
magnetic one? Since none of the mathematical 


7P. W. Anderson, Phys. Rev. 83, 1260 (1951). 
( — Proc. Roy. Acad. Sci. (Amsterdam) 39, 190 
193 
®P, W. Anderson, Phys. ae a 694 (1952); see also 
R. Kubo, Phys. Rev. 87, 568 (1952). 
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treatments of the problem seem to answer this 
question in a clear-cut physical way, we believe 
that the simple pictorial models given in this 
paper will be of considerable aid. From these 
models we shall derive the. two dispersion laws 
(7) and (11) in a very elementary fashion. 

Further, these models clarify the theory of 
antiferromagnetic resonance,!° which becomes, in 
a sense, a branch of spin-wave theory—at least 
in the low temperature region where spin-wave 
theory is applicable. 


II, SPIN WAVES IN A FERROMAGNET 


The dispersion law w « k? holds for one-, two-, 
and three-dimensional ferromagnets, only the 
proportionality factor changing. For simplicity 
we shall consider a linear ferromagnetic chain. 

Such a chain of N spins is shown in Fig. 1(a). 
We may imagine a small field Hy applied in the 
—z direction. Since the spin of an electron is 
directed opposite to its magnetic moment yp, the 
spins will all point up, as shown. We have 


S:=ui/g8=ui/hy, (12) 


where B=eh/2mc is the Bohr magneton; y =ge/ 
2mc is the magnetomechanical ratio; g is the 
spectroscopic splitting factor, approximately 
equal to 2 in ferromagnets. Owing to the negative 
sign of the electronic charge e, S; is oppositely 
directed to wi. 

At 0°K all the spins will precess in phase about 
z at the Larmor frequency, w=yHo. The z com- 
ponent of the total spin, S*, will be a maximum. 
That the individual spins cannot, however, ac- 
tually point in the z direction, but must have x 
and y components, is a well-known consequence 
of quantum mechanics; this has been shown by 
Klein and Smith" to lead to a zero-point spin- 
wave energy. 

If now we excite a spin wave of wave number 
k in our linear chain we have the situation as 
shown in Fig. 1(a), which might be considered a 
“snapshot” of a spin wave taken at a given time. 
The amplitude of precession, measured by R, is 
so increased on the excitation of a single spin 
wave that the z component of total spin, S*, 
decreases by one unit. The spins are also no 
longer in phase, the phase angle between succes- 


10F, Keffer and C. Kittel, Phys. Rev. 85, 329 (1952). 


1M. J. Klein and R. S. Smith, Phys. Rev. 80, 1111 
(1950). : 
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Net |). 


b= 1/2 k*® R/S, 


--} R 
Ho 
(b) TOP (c) SIDE 


Fic. 1. Semiclassical picture of a spin wave of wave 
number & in a linear ferromagnetic chain of N spins. 1(a) 
The N precessing spins as seen in perspective at a given 
time. The spins are precessing clockwise as seen from the 
top, and each successive spin is slightly out of phase with 
the preceding one. 1(b) A top view of three spins collapsed 
onto the same precessional circle. Arrows indicate direction 
of precession. 1(c) A side view of three spins, also collapsed. 
Owing to this spin wave the precessional frequency of the 
spins is increased over the Larmor frequency by an amount 
proportional to the angle @ (see text). Hence, the ferro- 
magnetic dispersion law w « k?. 


sive spins being equal to ka, where a is the dis- 
tance between spins (lattice constant).Periodic 
boundary conditions require the N+1 spin to be 
at the same phase as the first; hence, ka must be 
an integral multiple of 27/N. This is the same 
condition as was deduced below Eq. (6). 

In Fig. 1(b) we have collapsed three spins onto 
the same precessional circle, as seen from the top, 
in order to exhibit the effect of the phase relation- 
ships. A side view, also collapsed, is given in 
Fig. 1(c). It is to be noted that the angle ka is 
very much exaggerated in the drawings; for 
small k we may well approximate the sagitta by 
3k?a?R. The angle ¢ shown in Fig. 1(c) is approxi- 
mately equal to $k?a?R/So, since R must be small 
compared to the length So of a spin vector in 
order that the semiclassical approximation ob- 
tain. Here Sp=(S?+S)!, and in the semiclassical 
approximation So~S. 

We now show by classical and by quantum- 
mechanical arguments that, owing to the excita- 
tion of a spin wave, the increase in w over the 
Larmor frequency is proportional to the angle ¢ 
shown in Fig. 1(c), that is, to the angle between 
a spin vector and the resultant vector of its two 
neighbor spins. 


Classical Argument 


In an effective magnetic field Herp each spin, 
because of its associated magnetic moment, will 
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experience a torque uiXHers=vhS:XHerr. This 
torque is equal to the time rate of change of 
angular momentum, or 


d(hS;) /dt = yhS; X Het. (13) 


We may identify Her: from the Hamiltonian of a 
single spin, including the Zeeman term, 


K;= —gBS;- [Hot (2J/g8) >; S;] 

= — gBS;- Herr. (14) 
The j summation is over the nearest neighbors 
to 7. From Eqs. (13) and (14) we have 


d(hS,) /dt = yhS;X[Ho+ (2J/g8) £5 S;]. (15) 


2,4 
d= 1/8 kK PR/Sy 


(b) TOP (c) SIDE 


Fic. 2. Semiclassical picture of spin waves in an anti- 
ferromagnet. 2(a) The two sets of spin waves 1, 3, 5--: 
and 2, 4, 6--+ as seen in perspective at a given time. 
2(b) Top view collapsed showing phase relationships. The 
arrows indicate precessional directions in absence of Larmor 
precession. 2(c) Two side views, the first in the plane made 
by 2 and the resultant of 1 and 3, the second in the plane 
made by 3 and the resultant of 2 and 4. Since the exchange 
forces make the down spins precess oppositely to the up 
spins (see text), this is not the correct picture for an anti- 
ferromagnetic spin wave. 


Quantum- Mechanical Argument 


We follow Herring and Kittel.* The quantum 
equation of motion for §; is 


th(dS;/dt) =[S:, #] 
=(S:, —28 ©; 8;-Ho —2J Dy,2)8;-Sk] 
=[g8(X0S;-Ho) +2J(XS;-S,) 8S; 
—S,[g8(XS;-Ho)+2I(XS;-Sx)] 
=D; {[¢8(Si-Ho)+2J(S;-S,) JS; 
—S,[g8(S;-Ho)+2J(S;-S;) J} 
=(S;xS,] x [g8Ho+2J ¥; S;]. 


12 In obtaining the total exchange Hamiltonian of Eq. (1) 
from summing Eq. (14) over i, one must be careful to sum 
only over combinations of nearest neighbors in order not to 
count mutual energy terms twice. 
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On using the commutation relation SXS=7S and 
the relation g8 = hy we obtain the classical result, 
Eq. (15). The classical and quantum-mechanical 
equations of motion are the same due to the 
well-known relationship between commutators 
and Poisson brackets. 

The right-hand side of Eq. (15) is the torque 
acting on a spin 7. If the spins precess together 
about the z axis, always maintaining the same 
phase relationships, this torque remains constant 
in magnitude throughout the precession. Such 
a mode of motion must exist due to symmetry of 
the Hamiltonian about z; we call this motion a 
normal mode. From Eq. (15) and Fig. 1(c) the 
magnitude of the torque is easily found to be 


T= hySolLHo sin (S;, Ho) ats (2J/gB)2So sing | 
= hyLRHo — (2J/g8)Sca*k*R). (16) 


It should be noted that both Sp and R as used 
here are dimensionless. From Eq. (15) one finds 
the angular velocity of precession due to this 
torque, 

hwR=T; (17) 
and hence, 


wyHo —(2I/h) Soak’. (18) 


We have used the right-hand-screw sign conven- 
tions: T positive when into the plane of Fig. 1(c); 
likewise w positive when clockwise about +3. As 
y is negative the precession will be clockwise 
about —z, or when seen from the top as in 
Fig. 1(b). 

Since the excitation of a spin wave reduces S*, 
and since Hy is along —z, it is seen that the energy 
in the spin wave (18) is —hw. 

We have thus given a simple derivation of the 
dispersion law, Eq. (7), to which we have added 
a Larmor term. We note that the factor S in-the 
constant of Eq. (7) is replaced by Sp in Eq. (18), 
as was to have been expected from a semiclassical 
approximation. 


Ferromagnetic Resonance 


In case a spin wave with k=0 is excited, the 
amplitude of precession will change but the phase 
relationships will remain unchanged. From Eq. 
(18) it is seen that the frequency of such a spin 
wave is the Larmor frequency. Excitation of 
k=0 spin waves takes place in ferromagnetic 
resonance, in which an rf field is applied per- 
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pendicular to Hy and of such wavelength that all 
spins remain in phase. Our pictorial representa- 
tion makes it clear why the exchange interaction 
has no effect on the ferromagnetic resonance 
frequency, which in the absence of demagnetizing 
fields is simply the Larmor frequency. 


Ill. SPIN WAVES IN AN ANTIFERROMAGNET 


Spin waves in an antiferromagnet are shown 
in Fig. 2(a). We may begin by setting up two 
sets of spin waves, one in the sublattice pointing 
up (spins 1, 3, 5, etc.), and one in the sublattice 
pointing down (spins 2, 4, 6, etc.). These sets, 
however, are not independent, owing to the 
exchange interaction between nearest neighbor 
spins. 

In Fig. 2(b) a top view is given of several spins 
collapsed onto the same precessional circle. Side 
views, also collapsed, are given in Fig. 2(c). 

We apply Eq. (15) to an up spin and a down 
spin, obtaining with the aid of Fig. 2(c) two 
equations analogous to Eq. (16): 


T= hyLRHo+ (J/2g8)Soa?k*R], 


Ta=hy[ —RHo+(J/2g8)Soa*k?R]. (19) 


Here 7, and T, are the magnitudes of the torques 
on an up spin and on a down spin, respectively. 
Again the sign convention is T positive when into 
the plane of Fig. 2(c). Maintaining the conven- 
tion w positive when clockwise about +2, we have 


hw,R=T,, 


—hwagR= Tu. (20) 


Thus, the angular velocities of precession are 


= yHo+ (J/2h)Soa*k’, 


wa = ¥Ho —(J/2h) Soak’. (21) 


But the exchange parts of these precessional fre- 
quencies are in opposite directions! The spins 
will not precess together about the z axis always 
maintaining the same phase relationships, as 
was postulated in assuming the magnitude of 
torque (19) was a constant of the motion. How- 
ever, as in the ferromagnetic case, the Hamil- 
tonian has symmetry about g; and hence, normal 
modes must exist in which phase relationships 
are maintained. 
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The Normal Modes 


It. is easiest to visualize the normal modes in 
the absence of Ho. Then, as seen from the top, all 
spins will precess clockwise (mode 1), or all spins 
will precess counterclockwise (mode 2). As we 
shall see, these modes can be achieved only by 
having the up spins precess in a larger (mode 1) 
or in a smaller (mode 2) circle than do the down 
spins. The normal mode 1 is illustrated in Fig. 3. 
The normal mode 2 would be illustrated by Fig. 
3(a) turned upside down, except that the pres- 
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Fic. 3. Semiclassical picture of normal mode 1 of spin 
waves in an antiferromagnet. The up — precess in 
circles of radius R, the down spins in circles of radius P, 
with R>P. 3(a) As seen in perspective at a given time. 
3(b) Top view collapsed showing phase relationships. Posi- 
tions marked —1, —2, —3, —4 are those of imaginary 
vectors directed opposite to spins 1, 2, 3, 4, respectively. 
Arrows indicate precessional directions, which are now all 
clockwise. 3(c) Two side views. Note that 2 is on the oppo- 
site side of the line made by 1 and 3 as compared with 
Fig. 2(c). This is why the down spins, 2 and 4, are now 
precessing clockwise. This mode leads to antiferromagnetic 
dispersion law w« k, 


ence of Hy causes the two modes to be slightly 
different. That is, the Larmor precession’ | is al- 
ways clockwise about Ho, seen from the top, and 
hence will add to the clockwise precession of 
mode 1 but subtract from the counterclockwise 
precession of mode 2. 

In the absence of Hp the two normal modes are 
degenerate in energy and various linear combina- 
tions could be considered as ‘‘normal’’ modes. 
One such combination was obtained by Ander- 
son® when he used a canonical transformation to 
put his Hamiltonian into normal coordinate form 
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for two sets of harmonic oscillators. The canoni- 
cal transformation is a mathematical way of 
making the transition from Fig. 2 to Fig. 3. That 
Anderson’s transformation is not the correct one 
in the presence of Hp has been demonstrated by 
one of us, and the correct transformation has 
been given.!* 

In Fig. 3(b) we show a top view of the preces- 
sional circles. Shown on the circles are the phase 
positions of the spins 1, 2, 3, 4 and of the corre- 
sponding oppositely directed —1, —2, —3, —4. 
Side views are given in Fig. 3(c). It is clear from 
the figures that 1 and 3 act on 2 through an angle 
¢; determined by x, whereas 2 and 4 act on 3 
through a slightly larger angle ¢2 determined by 
y. It is also to be noted that the vector cross 
product of the resultant of 1 and 3 with 2 is 
oppositely directed to its analog in Fig. 2(c). 
This accounts for the fact that the down spins 
now precess in the same direction as the up spins, 
and hence, we have a normal mode. 

We have a normal mode, that is, if the angular 
velocities of precession of the down spins is the 
same as that of the up spins. We now find these 
velocities. The torque equation (15) yields 


T= hyLRHo+ (2J/g8)2ySo], 
where R and P are the radii of the circles of pre- 
cession of the up and down spins, respectively. 
Instead of Eq. (20) we now have 


ho,R=T., 


—hwgP = Tu. (23) 
On equating w, and w, as obtained from inserting 
Eq. (23) into Eq. (22), we have as the condition 


13F, Keffer, thesis, Berkeley, January 1952. We repeat 
here the transformation. Anderson’s Eq. (15) (see reference 
9) is replaced by 
Py =—n(i-—2)“prt+(1—2)“" da, 
Qn=n(1+n)“gr+(1+7)“g2r, 
Ry=(1+2)"ga+n(1+n)1ga, 
Sy= (1-2) pn-al(l—1)“pa, 
where 
n=—a—|(a?—1)4], 
with, using Anderson’s notation, 
e=[1+(K/DJ)]/n. 


The normal modes in the presence of Ho have also been 
derived (using a different formalism) by J. M. Ziman, 
Proc. Phys. Soc. (London) A65, 540 (1952); and by T. 
Nakamura, Prog. Theoret. Phys. 7, 539 (1952). 
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for a normal mode 


y/R=x/P. (24) 


Since the phase angle between spins 1 and 3 is ka, 
we have from Fig. 3(b) 


x= R—P —(k’a?R/8) ; (25a) 


similarly, 
y=R-—P+(k’a?P/8). 


From Eqs. (24) and (25) 
R-P=ka[(R?+P?)/8]}, 
and after a bit of algebra 


ka(P+R) 


(25b) 


(26) 


+ x 


R P (8R+8P) 


In this expression we may set R=P, for small k, 
by virtue of Eq. (26); thus, | 


y/R=x/P=ka/2. 
Finally, from Eqs. (22), (23), and (27) 
@u” = wa =yHo+(2J/h)Soak (mode 1). 
Similarly, we obtain 
wy = wa =yHo —(2J/h)Soak (mode 2). (28b) 


In determining the direction of these precessions 
one must recall that both y and J are negative. 
Sign conventions used here are given below 
Eq. (18). 

Thus, we have demonstrated the dispersion 
law for antiferromagnets, Eq. (11). 

It should be noted that the energy of a spin 
wave in mode 1 is —iw“ whereas in mode 2 it is 
+hw®. This will be clarified in Sec. IV. 


(27) 


(28a) 


IV. ANTIFERROMAGNETIC RESONANCE 


On exciting a spin wave of mode 1 we reduce 
the z component of the total spin of the up sub- 
lattice by more than we increase the z component 
of total spin of the down sublattice. It can be 
shown!* that the net effect is to decrease the z 
component of total spin of the entire lattice, S*, 
by h. Similarly the net effect of exciting a spin 
wave of mode 2 is to increase S* by h. Thus, 


S?¥= Die (Mx —Mux), (29) 


where my, and my, are the number of spin waves 
of wave number k excited in mode 2 and mode 1, 
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respectively. This equation may be obtained by 
an extension of Anderson’s work® on use of the 
correct transformation in the presence of Hy.!*:14 

The statement about the energies of the spin 
waves made at the end of Sec. III becomes clear 
in the light of Eq. (29). Since Ho is directed 
along —z, energy is required to increase S*. 

As in ferromagnetic resonance, discussed in 
Sec. II, an rf field applied perpendicular to Ho 
will excite only k=0 spin waves. It will excite 
equal numbers in mode 1 and mode 2; hence, S* 
cannot change, and no energy will be required. 
That is, an indefinite number of spin waves of 
k=0 will be excited in the presence of no rf 
energy whatever! The spins will constantly be 
turning over. This has been pointed out by 
Anderson,® who introduces some anisotropy to 
keep the sublattices from interchanging. 

Probably anisotropy is always present in anti- 
ferromagnetism. In MnF2, for example, one of 
us'§ has calculated an effective anisotropy field 
of 8800 oersteds tending to keep the spins lined 
up. We need to modify our spin-wave picture to 
take account of this. We introduce an effective 
anisotropy field in the negative z direction acting 
on the up spins (whose magnetic moments are 
down), and one in the positive z direction acting 
on the down spins. Just as with the exchange 
fields of Fig. 2 these fields will tend to make the 
members of the two sublattices precess in opposite 
directions. In order that the spins precess in the 
same direction they must precess in different size 
circles. Thus, whether or not the spins are all in 
phase, i.e., even for k=0 spin waves, the spins 
will be inclined to one another. To perform this 


14 It is to be noted that Eq. (29) cannot be obtained from 
Anderson’s transformation. 


15 F. Keffer, Phys. Rev. 87, 608 (1952). 
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inclination requires work against the exchange 
field Hz as well as against the anisotropy field H.. 
We use the notation Hg=|(2J/g8) +; S;|; see 
Eq. (14). It has been shown by one of us using 
Anderson’s formalism!* and by one of us using 
the semiclassical formalism" that this amount of 
work (for small k) is equal to ” 


Ey = Fhe, = FhyHy —hy(2HeH a 
+Ha4?+const a?k*)'. (30) 


In + the — is for mode 1, the + for mode 2. The 
value of the constant for a three-dimensional 
simple cubic is H,z?/12; for a one-dimensional 
chain it is Hz?/4. For Ha=0, Eq. (30) reduces 
to Eq. (28). 

To excite a spin wave of k=0 in a resonance 
experiment we must match the precessional fre- 
quency given by Eq. (30), 


wo= yHoty(2HeHst+H,?)?. (31) 


This resonance frequency has been derived from 
a molecular field approach by Kittel'* and by 
Nagamiya.'® The normal modes for resonance, 
very similar to the normal modes shown here in 
Fig. 3, are given in Fig. 1 of Keffer and Kittel." 
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16H. Kaplan, Phys. Rev. 86, 121 (1952), and thesis, 
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17 See the remarks at the end of Sec. III relating fre- 
quencies to energies in the two modes. 
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Cerenkov Radiation in a Dispersive Medium 


H. Motz anp L., I. ScHirr 
Microwave Laboratory and Department of Physics, Stanford University, Stanford, California* 
(Received October 10, 1952) 


Attention is called to the disagreement between a recent statement of Sommerfeld that the 
direction of propagation of Cerenkov radiation in a dispersive medium is determined by the 
group velocity, and the generally accepted result that the phase velocity must be used. The 
latter is shown to be correct, and the role of the group velocity is discussed in a way similar to 
that given in Tamm’s 1939 paper on the theory of the Cerenkov radiation. 


N a recently published book on optics, Som- 
merfeld! states that the direction of propaga- 

tion of Cerenkov radiation in a dispersive me- 
dium is given by the usual formula? if the 
velocity of light in the medium is taken to be 
the group, not the phase velocity. This is con- 
trary to the generally accepted result, according 
to which the phase velocity must be used. The 
difference between these two results can be of 
considerable physical significance in an experi- 
ment like that recently performed by Mather.’ 
That the use of the phase velocity is actually 
correct can readily be seen in the following way. 
For simplicity, we restrict our attention to an 
infinitesimal range of radiated frequencies; this 
might be achieved in practice by use of a narrow- 
band filter. Then the usual calculation shows 
that the Poynting vector for each monochro- 
matic component has the direction calculated 
from the phase velocity. All of these Poynting 
vectors are infinitesimally close together, so that 
the direction of propagation of the entire range 
of frequencies is infinitesimally close to that 
calculated for any one frequency component in 
the range by use of the phase velocity. 

It is still of some interest, however, to see as 
clearly as possible the role played by the group 
velocity. This point has already been discussed 
by Tamm‘ in his classical paper on the theory of 
the Cerenkov radiation. He shows that for a 

* The work reported herein was performed in part under 


a contract between Stanford University and the U. S. 
Office of Naval Research. 

1A. Sommerfeld, Vorlesungen tiber Theoretische Physik, 
~~ IV, Optik (Dieterich Verlag., Wiesbaden, 1950), 
p. . 

2 See, for example, reference 4, or L. I. Schiff, Quantum 
Mechanics (McGraw-Hill Book Company, Inc., New York, 
1949), P: 264. 

*R. L. Mather, Phys. Rev. 84, 181 (1951). 

‘I. Tamm, J. Phys. U.S.S.R. 1, 439 (1939); see espe- 
cially Section 5. 


narrow range of frequencies, the radiation at any 
instant is concentrated in a thin conical shell 
whose vertex is at the moving charge, and whose 
half-angle of opening @ is given by 


coté =[(mov/c)? —1 }! 
+-wono(v/c)?(dn/dw) ol (nov/c)*-1]}-#, (1) 


where v is the velocity of the charge, c the 
velocity of light in vacuum, 2(w) the refractive 
index of the medium at angular frequency w, 
and the frequency range is centered at wo with 
no=Nn(wo). If the radiation were to propagate in 
directions perpendicular to the planes tangent to 
the surface of this cone, it would agree with the 
correct theory only for a nondispersive medium 
(dn/dw=0). Actually, the radiation propagates 
in directions that make the angle 4 with the path 
of the charge, where 


tan6y = [(nv/c)? —1 ]}, (2) 


so that the cone may be thought of as ‘“‘side- 
slipping”’ as it moves along with the charge. 

It seems worth while to recast Tamm’s argu- 
ment in a slightly different form, which brings 
out the analogy with the well-known derivation 
of the group velocity of a one-dimensional wave 
packet.5 For simplicity, we regard the field as 
varying only in a plane through the path of the 
charge, so that only two dimensions need be used. 
A general wave packet solution for any field com- 
ponent is then 


u(x, y, t) = f f A (kz, ky) 


Xexp[i(kax+kyy —wt) ]dk.dk,, (3) 


where w=c(k2+k,?)'/n(w); as before, the inte- 


5 See, for example, M. Born, Atomic Physics (Blackie 
and Son, Ltd., Glasgow, 1946), p. 295. 
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gration is over a narrow range of frequency or 
propagation vector. We expand the integrand 
about the center of the range wo, Ko, and rewrite 
Eq. (3) as the product of a slowly varying ampli- 
tude and a rapidly varying phase factor: 


u(x, y, t) = B(x, y, t) 
Xexpli (Rox + kay —wot) ], 


B(x, y, 1) = f f a(kis, Ry) m 


Xexp[t(k izt+ Ruy —w),t) ]dk izdk ly» 
k=k)+ki, w=otw1, 
a(kiz, Ry) =A (kz, ky). 


The amplitude B is a slowly varying function of 
x, y and ¢ since k, and w, are restricted to narrow 
ranges of values. The energy is concentrated in 
regions where B is large; the motion of these 
regions may be found by writing down the condi- 
tion that B=constant, which has the differen- 
tial form 


Bax+B,dy+Bdt=0, (5) 


where B,=0B/dx, etc. Equation (5) describes a 
series of parallel lines of slope —B,/B,; the inter- 
cept of any one of these with the x axis moves 
with the velocity —B,/B., and the y axis inter- 
cept with velocity —B,/B,. From Eq. (4), we 
have 


a f f pee hud 


B,=if fea exp( )dkidkwy, 


(6) 


B,= -if fow exp( )dkidkyy. 


So far, the wave packet treatment is quite 
general. In the Cerenkov case, A(kz, ky) vanishes 
unless ky =k.[(nv/c)? —1]!, so that w=kw. This 
follows from the usual theory, or may be seen 
physically from the requirement that the radia- 
tion pattern be stationary with respect to the 
charge moving in the x direction, which means 


IN A DISPERSIVE MEDIUM 


Fic. 1. A charge moves in the direction OX, and a line 
of constant radiation amplitude moves from AC to BD in 
time At. The direction of energy flow OR is such that the 
distance EF is equal to At times the group velocity of the 
radiation. 


that the phase velocity component in that direc- 
tion, w/kz, must be equal to v. Then in the first 
of Eqs. (6), Riz can be replaced by w:/v, and in 
the second ki, can be replaced approximately by 


wild{ (w/v) [(mv/c)? —1]*} /dw]o 
= (w1/v) {L(mov/c)? —1]! 
+ wono(u/c)*(dn/dw) of (nov/c)? —1 }-3}. 


These substitutions lead to 


p= —B,/v, 
B,= —(B./v){L(mov/c)? -1]}! 
+ wonto(v/c)?(dn/dw) of (mov/c)? —1]-3}. 


Thus the lines (actually generators of cones) of 
constant amplitude are inclined to the x-axis 
with the angle @ given by Eq. (1). 

The direction of flow of the energy may be 
found from the construction of Fig. 1. The 
charge moves in the direction OX, and a line of 
constant amplitude has the position AC at one 
instant of time and the position BD at a time 
At later. From the discussion of Eq. (5), we 
have AB= —(B,/B,)At, and CD= —(B,/Bz)At. 
The radiated energy must move in the direction 
OR, which is so determined that its point of inter- 
section with the line of constant amplitude 
moves with the group velocity u of the radiation; 
that is, the distance EF must be equal to uAt. 
Now wn(w)=ck, where k is the magnitude of 
the propagation vector, and u=dw/dk; therefore 
u=c/[no+wo(dn/dw)o]. It can then be shown 
quite readily that the angle ROX in Fig. 1 is just 
equal to the angle 0 given by Eq. (2). 





Behavior of Magnetic Materials 
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This is a review of recent work in which the atomic theory of ferromagnetism and the domain 
theory of magnetization are applied to new materials. 


oe two new types of magnetic ma- 
terials have come into prominence, which 
illustrate very well two important problems in 
ferromagnetism that have existed for some time. 
These problems are: 

(1) What is the atomic structure of magnetic 
materials, and what relation does this structure 
have to magnetic properties? 

(2) What is the mechanism of magnetization 
whereby the magnetization of a material is 
changed by application of a magnetic field? 

These will be discussed in relation to two new 
kinds of materials: (a) ferrites, and (b) fine par- 
ticle magnets, both of which are of considerable 
commercial importance as well as of scientific 
interest. 

First, some general remarks on atomic struc- 
ture are needed. For a material to be ferromag- 
netic its atoms must possess a permanent mag- 
netic moment; that is, each atom must be a small 
magnet. Also, there must be some sort of inter- 
atomic force which maintains the moments of 
neighboring atoms parallel; otherwise, the ma- 
terial would be paramagnetic, that is, only 
weakly magnetic. 

Figure 1 shows the structure of an iron atom. 
It is well known to be built up of electron shells, 
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Fic. 1. Electron shells in an iron atom, showing location 
of unfilled 3d shell responsible for the magnetic moment 
of the atom, and for the ferromagnetism of the metal. 


as shown, and in each shell the electrons have 
both orbital motions around the atomic nucleus 
and spins about their own electron axes. The 
first two shells are complete. They contain as 
many electrons as these shells can ever contain, 
and the magnetic moments of the electrons in 
each shell are balanced so that the net magnetic 
moment is zero. However, in the third shell one 
of the subshells is completely filled, whereas its 
counterpart is only partially filled, so that we 
have a net magnetic moment corresponding to 
the spin of 5—1 or 4 electrons. This atom then 
has an atomic moment of 4 Bohr magnetons. 

This is the well-established structure of a free 
atom of iron. When this atom becomes part of a 
solid metallic structure, the distribution of elec- 
trons is disturbed, and we know by experiment 
that each atom has on the average a magnetic 
moment of 2.2 Bohr magnetons. We do not 
understand the structure of a metal well enough 
to know why this number should be what it is. 

Relatively little is known about the inter- 
atomic forces that maintain the moments of 
neighboring atoms parallel. However, Heisenberg 
showed in 1928 that it is a quantum-mechanical 
exchange force which depends on interatomic 
distance. In 1930 Slater found that the distance, 
D, between neighboring atoms, divided by the 
diameter, d, of the incomplete shell, in which the 
magnetic moment resides, is larger in the ferro- 
magnetic elements than in the nonmagnetic ones. 
Bethe suggested that the atomic moments in a 
crystal are maintained parallel by quantum- 
mechanical forces only when they are separated 
by more than a minimum distance, and that 
when the distance increases beyond this the 
forces become first larger and then smaller again. 
However, the form of this interaction curve, 
relating the force of ferromagnetism to atomic 
distance, is not at all well established. 
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DOMAINS 


Now, as neighboring atoms are aligned parallel, 
why is not any piece of ferromagnetic material 
completely and permanently magnetized in one 
direction? In 1907 Weiss suggested an explana- 
tion which is illustrated in Fig. 2(a) and which 
in principle has been very well substantiated 
since its proposal in 1907. According to this idea 
small regions or domains are magnetized to satu- 
ration in various directions so that the net mag- 
netization of the whole specimen is zero, before 
a magnetic field has been applied. 

Since Weiss’ time we have come to know that 
there are certain directions in a crystal which are 
preferred, and in iron we know that these direc- 
tions are the six directions which lie parallel to 
the cubic axes of the crystal. A revised domain 
picture then is shown in Fig. 2(b), where the 
directions of the domains in a single crystal are 
limited but still the net magnetic moment of the 
whole specimen is zero. 

Recently H. J. Williams! and his colleagues 
have been able to see under the microscope the 
actual domains in a crystal specimen and detect 
the direction of magnetization in each domain. 
A photograph of this kind is shown in Fig. 2(c), 
and other specimens are shown in some of the 
following figures. The technique used is similar in 
principle to the detection of cracks or flaws in a 
large piece of iron by sprinkling iron filings over 
it. When a colloidal suspension of Fe;O, is placed 
on a crystal whose surface is carefully prepared, 
the particles of Fe;O, collect along domain 
boundaries where there are small stray fields, and 
they appear under the microscope as fine light 
lines. 

When a field is applied, the domain structure 
of Fig. 2(c) changes to that of Fig. 2(d); that is, 
the zig-zag domain wall shifts in position. This 
is the well-known mode of magnetization in weak 
fields—a change in magnetization takes place by 
movement of the boundary separating domains 
so that the domains which are aligned more 
favorably with respect to the field (those above 
the line) are enlarged at the expense of other 
domains (those below). Such changes have been 
shown in a striking way in moving pictures. 


1 Williams, Bozorth, and Shockley, Phys. Rev. 75, 155 
(1949). 


SCHEMATIC 


(a), 
CRYSTAL AXES XX H—> 


Fic. 2. Domain diagrams (a) and (b), and actual 
photographs (c) and (d), of domains in an iron crystal 
(Williams). 


Many domain structures have been observed 
under the microscope, and we can interpret many 
of them, although not all. We are, however, quite 
confident that we know the principles which 
govern domain formation, and these are illus- 
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(d) 


Fic. 3. Possible domain structures, showing large mag- 
netostatic energy in (a), and successively lower energy in 
other structures in which there is wall energy (b) and 
magnetostrictive energy (c); (d) represents the kind of 
structure actually observed (Fig. 4). 


trated in Fig. 3.2 If we have a block of magnetic 
material as shown in Fig. 3(a) and if this is 
magnetized as a single domain, then we will have 
a large external field as indicated by the curved 
lines. Of course, a considerable amount of energy 
is stored in this field. However, if we divide this 
block into two portions by a domain wall so that 
the two halves are magnetized in opposite direc- 
tions, the energy in the external field is obviously 
greatly reduced. We can also make the domain 
structure shown in Fig. 3(c), where the exterior 
field is completely eliminated and the domains 
are always magnetized in a direction of easy mag- 
netization parallel to a crystal axis. But here an 
additional energy is present. Besides the energy 
associated with a domain wall, which in iron is 
about 1 or 2 ergs per square centimeter, we have 
energy associated with magnetostriction. When 
a domain becomes magnetized it changes its di- 
mensions slightly by a few parts per million, and 


?R. M. Bozorth, Ferromagnetism (D. Van Nostrand 
Company, Inc., New York, 1951). This book may also be 
consulted for further information about various topics 
discussed in this article. 
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the domains which are magnetized in a horizontal 
direction in the figure do not quite fit into the 
spaces allotted to them. As a result the structure 
is slightly strained, and we can calculate the 
amount of energy associated with this strain. In 
the structure shown in Fig. 3(d) this magneto- 
striction energy has been reduced considerably 
by reducing the amount of strained material, but 
the additional wall area has increased the energy 
associated with it. 

The stable structure which is produced in 
nature is, of course, that for which the total 
energy is minimum. This requires a calculable 
balance between the volumes with horizontal 
magnetization and those with vertical magnetiza- 
tion. The scale of the structure derived by theory 
is actually very close to that observed. Figure 4 
shows the powder pattern of the domain struc- 
ture observed in an actual crystal, and its re- 
semblance to Fig. 3(d) is obvious. 

Another type of domain pattern is shown in 
cobalt. This material has a hexagonal crystal 
structure, and magnetic measurements show that 
the hexagonal axis is a direction of easy mag- 
netization. Unlike cubic iron, which has six direc- 
tions of stable magnetization, cobalt has just two 
(antiparallel) directions. Figure 5(a) shows a sur- 
face of a cobalt crystal which has been cut so 
that it contains the crystal axis. The domain 
structure of this material is like a bundle of 
needles with the axes of the needles all parallel 





Fic. 4. Observed magnetic domain pattern corresponding 
to Fig. 3(d). Edge of crystal shows near bottom of figure. 
Length of specimen shown is about 0.6 mm. Retouched; 
arrows added to show local direction of magnetization. 
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Fic. 5. Domain patterns of cobalt on surfaces parallel (a) and perpendicular (b) to hexagonal axis, 
showing rod-like domain structure. 


and the magnetization of neighboring domains 
lying in the two antiparallel directions. A domain 
pattern, Fig. 5(b), on a surface cut so that it is 
normal to the axis of the crystal, shows the ends 
of the needle-like domains, as one would expect. 


FINE PARTICLES 


An interesting application of domain theory is 
concerned with the size of the individual speci- 
men. If we take a cube of the material, we may 
assume either of two configurations as shown in 
Fig. 6. The first diagram shows the material as 
a single domain magnetized vertically. The sec- 
ond shows the cube composed of several domains 
with walls separating them. In the first case we 
have a large amount of magnetostatic energy 
associated with the exterior field arising from the 
north and south poles on the surfaces. The second 
diagram has no exterior field but has an energy 
associated with the domain walls. One can readily 
calculate the relative energies of the two con- 
figurations, and it is found that the magneto- 
static energy E; of the first cube is 


E, = 212d, 
and the wall energy E2 of the second is 
E.=2v2yd?, 


I, being the intensity of magnetization at satura- 
tion, d the length of the edge of the cube and 
the energy of 1 cm? of wall. The first is propor- 
tional to the third power of the cube edge, the 


second to the square of the edge. Normally in 
large particles the second configuration is more 
stable, and several domain walls are present. 
However, if the particle size is small enough, FE; 
will be less than E2 and the first configuration 
will be more stable, and so each particle will be 
a single domain. In iron the two energies are 
equal when the particles are about 300A in size. 

Now consider what happens to these different 
domain structures when a field is applied. When 
domain boundaries are present, application of a 
field will move the boundaries as shown in the 
upper half of Fig. 7, and this process will proceed 
in relatively low fields with great ease. However, 
if we apply a field to a single domain particle, 
this process cannot occur because no boundaries 
exist. Therefore, we must apply a field strong 
enough to rotate the magnetization out of the 
direction of easy magnetization. The field re- 
quired for this is determined by the nature of the 
crystal. In a cubic crystal like iron the energy 
required to rotate the magnetization is shown 


€, * 213d° E> + 2Vey7d" 

Fic. 6. Two possible domain structures for a single 
particle of iron: (a) stable in small particles, (b) in particles 
larger than about 300A, when E, ~ Ez. 
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CRYSTAL AXES: a 


Fic. 7. Top portion shows 
mechanism of change of mag- 
netization by displacement of 
domain wall and then by rota- 
tion of magnetization. Lower 
diagram shows mechanism when 
walls do not exist and mag- 
netization must be rotated over 
energy hump caused by crystal 
anisotropy. 


qualitatively by the clover-leaf pattern at the 
lower left of Fig. 7, which is consistent with the 
cubic symmetry of the system. The scale of the 
pattern is measured by the magnetic anisotropy 
constant, K, which can readily be determined 
experimentally. In iron the anisotropy is such 
that about 500 oersteds are required to rotate 
the magnetization from an easy direction to a 
hard direction. In a specimen containing many 
domains, on the contrary, a change in magnetiza- 
tion can be produced by a shift in the domain 
walls, and these move readily in a field of about 
1 oersted. 

The practical result of all this is that if we form 
iron into very fine particles, of the order of 200 
to 300 angstroms in diameter, we have made 
tiny permanent magnets of it. Such magnets 
have been made on an experimental basis by 
careful reduction of the metal oxide, and these 
have a coercive force of about 500 oersteds. Also, 
magnets made of pressed iron powder have been 
made on a commercial scale with quality that 
compares favorably with the best permanent 
magnets that were known 20 years ago. So far 
they cannot compete with our best modern alloys, 
but other materials in powder form have the 
possibility of exceeding our best commercial ma- 
terials in quality. This field is developing rapidly, 
and it is not impossible that powder magnets will 
soon be produced that will be superior to any- 
thing that we know of today. An important step 
in this direction is the recently reported Bis- 
manol,? a compacted powder of manganese bis- 


3 Adams, Hubbard, and Syeles, ‘‘Bismanol, A New Perm- 
anent Magnet.” Navord Report 2440, U. S. Naval Ord- 
nance Laboratory, Silver Spring, Maryland (May 20, 
1952), pp. 1-15. 


muthide which was found by Guillaud‘ to have 


a high crystal anisotropy and a high coercive 
force. 


FERRITES 


Now let us turn our attention to the ferrites 
which have recently come into commercial 
prominence particularly in applications where 
high frequency is important. These are refractory 
oxides composed of iron oxide, Fe2O3, chemically 
combined with one or more of the oxides of 
nickel, manganese, zinc, cobalt, etc. Usually they 
are pressed in powder form to the required shape 
and fired, after which they are hard, strong, and 
brittle. 

When any magnetic material is magnetized in 


108 


3 
ve 


0.001" 4-79 
Mo PERMALLOY 


ee 


VALUES OF 4% AND yQ 


2 
10 
10? 107 104 10° 
FREQUENCY IN CYCLES PER SECOND 


Fic. 8. Effective permeability u, and its product with the 
quality factor Q (reciprocal of energy loss) as dependent on 
frequency. Note difference between laminated metallic 
material and unlaminated ceramic (ferrite), and superiority 
of the latter at high frequencies (1 Mc/sec). 


a high frequency field, the power losses resulting 
from eddy currents become an increasingly im- 
portant factor. Ordinarily eddy currents are re- 
duced by laminating the material, and Fig. 8 
shows the way that the losses vary in a well- 
laminated structure of a high quality metallic 
material, Molybdenum Permalloy, in which the 
laminations are only 0.001 in. thick. The losses 
are given in terms of Q (the quality factor), 
the reciprocal of the loss. More specifically 
Q=2zfL/AR, L being the inductance at the 
frequency f, of a coil of wire containing the ma- 
terial, and AR the effective resistance of the coil 
caused by the presence of the material. 

4 C. Guillaud, “‘Ferromagnetism of binary alloys of man- 


ganese,” Thése, l’Université de Strasbourg (France) 1943, 
pp. 1-129. 
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In addition to the increase of loss at high 
frequencies, (AR), we have a decrease in the 
permeability, un, which arises because the whole 
flux cannot follow the rapid alternations of the 
field. The best single criterion for a material for 
use at high frequencies is the product of these 
quantities, uQ. The uQ curves for Molybdenum 
Permalloy and for manganese zinc ferrite are 
shown by the solid lines in Fig. 8. One sees that 
in 1-mil Molybdenum Permalloy the curve falls 
rather rapidly with increasing frequency, whereas 
in the curve for the ferrite (not laminated) it 
holds up well to over one megacycle. 

The efficacy of the ferrites is due to their very 
high electrical resistivity; they do not require 
lamination. Although ‘the resistivity of the Per- 
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Fic. 9. Theoretical and observed atomic moments of 
various ferrites (MnFe.O,, etc.), showing general validity 
of Néel theory based on atomic moments of Table I. 


malloy is about 55X10-* ohm-cm, that of the 
nonmetallic refractory ferrites is about 10? to 10°. 
The ratio for the two materials is then more than 
a million. Although the permeability of the 
ferrites at low frequencies is not high, their uQ 
is relatively large at frequencies of 1 megacycle 
and above. There have obviously been many 
applications of this material in communication. 

We are concerned here primarily with the 
structure of the ferrites and its relation to their 
magnetic properties. The chemical formula of the 
nickel ferrite (for example) may be written, 
NiO-Fe.O; or NiFe,O,. The crystal structure of 
these materials is determined first of all by the 
oxygen atoms, which are much larger than the 
metal atoms, and which form a close-packed 
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structure like that of a pile of billiard balls. In 
this structure there are two kinds of interstices, 
in which the metal atoms lie. One kind, the ‘‘A’”’ 
position, is equidistant from four oxygen atoms 
(tetrahedral holes); the other, or ‘‘B”’ position is 
surrounded by six oxygen atoms (octahedral 
holes). In a normal ferrite the bivalent atoms 
are in the ‘‘A”’ positions and the Fe*+** atoms in 
the ‘‘B” positions. These ferrites are nonmag- 
netic; the magnetic ferrites all have the inverse 
spinel structure, with the Fet**+ atoms equally 
distributed between the ‘‘A”’ and “B”’ positions, 
and all of’ the bivalent atoms, e.g., Nit* in 
NiFe.O,, in the “B’’ positions. This may be 
represented schematically as follows: 


A B 
Nit*+ (Fett, Fet**) O,-- 
Fet++ (Nit+, Fet++) O,-- 


~ —- &— 


(normal spinel), 
(inverse spinel). 


It was proposed by Néel® that in the magnetic 
inverse spinel the magnetic moments of the ions 
in the “A” positions are held by quantum- 
mechanical forces so that they are antiparallel 
to those in the “B”’ positions, as shown by the 
arrows below the symbols. The net magnetic 
moment of the ferrite molecule is then the mo- 


TABLE I. Atomic structure and magnetic moments 
of some iron-group ions. 
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® Net moment is 4 in a free atom, 2.2 in the compact metal. 


5L. Néel, Ann. phys., series 12, 3, 137 (1948). 
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ment of the bivalent metal ion, i.e., of Nit+* in 
N iFesOx,. 

The magnetic moments of the various ions 
that occur in the common ferrites are shown in 
Table I. These are derived from the atomic 
structure in much the same way as the moment 
of a free iron atom is derived from the structure 
shown in Fig. 1 and are simply the differences 
between the number of electrons in the third 
shell having plus spins and those having minus 
spins. In NiFe,O, the moments of the Fet*+ 
atoms cancel each other and the net moment of 
the molecule is that of the Ni*+ ions or 2 Bohr 
magnetons. 

The atomic moments of various ferrites, as 
predicted by atomic theory, can be compared 


with experiment by measuring the saturation 


magnetization of the materials at low tempera- 
tures and extrapolating to 0°K. Measurements 
have been made® for manganese, iron, cobalt, 
nickel, and copper ferrites, and Fig. 9 compares 
the number of Bohr magnetons deduced from 
the magnetic measurements with the predictions 


*For references to the sources of data see L. Néel, 
Comptes rend. 230, 190 (1950). 
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of Néel’s theory, as indicated by the straight line. 
The small discrepancy between the calculated 
and observed moments has been attributed to the 
magnetic moments caused by the orbital motions 
of the electrons, which are not taken into account 
by theory. 

We can see then that we have a very good 
basic understanding of the atomic structure of 
the ferrites and its relation to one of the prin- 
ciple magnetic quantities, the magnetization at 
saturation. 


DOMAIN MOVIES 


Further description of the processes that occur 
in magnetic materials, when a magnetic field is 
applied, are illustrated in the film “Action Pic- 
tures of Magnetic Domains,’ prepared by H. J. 
Williams, C. Kittel, and F. M. Tylee. A limited 
number of copies of this film are available for 
loan for educational purposes. Requests should 
be addressed to the Publication Department, 
Bell Telephone Laboratories, 463 West Street, 
New York 13, N. Y. The normal time for projec- 
tion (without sound) is about 30 minutes. 


Practical Aids for Physics Teachers 


The American Journal of Physics plans to introduce a 
new feature as soon as possible under the title, ‘Practical 
Aids for Physics Teachers.”’ At first, it may begin in a 
small way, but we hope that it will grow to occupy several 
pages in every monthly issue. 

The aids that go under such a title must be contributed; 
they cannot be prepared in the Editor’s office. Material 
of the following kind would be most welcome: 


A. Examination papers that seem to have been 
successful ones at all levels of difficulty, even into the 
beginning graduate level, but especially at the ele- 
mentary level. 

B. Problems, particularly at the elementary level; 
also at the upper undergraduate level. If problems are 
very difficult, solutions should be provided. They will 
be printed with the problems. 

C. Useful analogies that may be used to clarify 
some of the basic concepts of physics; or useful illus- 
trations of fundamental ideas. 

D. Neat ideas for lecture demonstrations. 

E. Useful material in books, magazines or trade 
journals that are not likely to come to the attention 
of the average teacher of physics but which offer 
illuminating discussions of simple matters. Give the 


reference completely and specifically; summarize in 
one or two sentences. 


The contributor of each item will be recognized in the 
following way: At the end of each item (or it may be at 
the beginning) will be placed the following statement: 
Contributed by John Doe, Penpoint College. This form will 
be used if the contributor wishes to be identified. 

Should the contributor not wish to be identified, a form 
such as the following should be employed, for example, in 
the case of a contributed examination paper: Used in a 
first-year physics course for liberal arts students in a private 
college of approximately 800 total enrollment. Or else it 
might run: Used in a physics course for sophomore engi- 
neering students in a state university of total enrollment 
15 000. 

Now you see what I mean. Please send me contributions. 
The only stipulation that I must make is that everything 
must be typed double-spaced. When I say everything, I 
mean everything. Otherwise, there is no room for editorial 
markings. A mimeographed copy that is already double- 
spaced and of good quality would be acceptable. 

This is your Journal, your Association. Won't you help? 

Tuomas H. OsGoop 
Editor 
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A ring stands, with axis horizontal, on a horizontal table and a cylinder is placed, also with 
axis horizontal, resting within it. If motion in the same horizontal direction be given to both, 
they both roll and the motion is of no particular interest. If, however, motion in opposite direc- 
tions be given to the ring and the cylinder ‘“‘rocking” occurs. This motion is interesting because 
its theoretical solution involves the use of Lagrangian equations, there being two degrees of 
freedom. Theory and experiment are here given for rings (thin and thick, complete and semi- 


circular) and cylinders (hollow and solid) of different sizes. Good agreement is observed. 


I. CYLINDRICAL SHELL 


HE experiments described in this paper 

were prompted by previous experiments 
on rocking! and by a problem set in Lamb’s 
Higher Mechanics.’ It reads, 


A thin cylindrical shell of radius R and 
mass M rests with its axis horizontal on a 
horizontal table. Inside it is placed a circular 
cylinder of mass m, radius 7, and radius of 
gyration k. If the system is set rocking, form 
the equations of finite motion. If the motion 
be small, prove that the length of the equiva- 
lent simple pendulum is given by 


2M(1i+k?/r?) 
L=(R-r)————_ : (1) 
2M+m(1+k?/r?) 


We assume no slipping between the shell (or 
ring) and the cylinder, nor between the ring and 
the table. The centers of gravity of the ring and 
cylinder move in opposite directions. 

Figure 1 illustrates the motion. In it, O, the 
center of gravity of the ring, is moving to the 
right; and P, the center of gravity of the cylinder, 
is climbing upwards on the left. In Fig. 1 T is 
the point of contact of the ring with the table and 
Q the point of contact of the ring with the 
cylinder. In the central position O, and P; are 
the centers of ring and cylinder, respectively, and 
T, the point of contact of ring with table and 
also of ring with cylinder. We see that in the 
motion which has ensued since the bodies last 
passed through the central position the point 


ass * eee Am. J. Phys. 18, 405 (1950) and 19, 
$11 (195 

2H. Lamb, Higher Mechanics (Cambridge ee 
Press, Cambridge, 1929), Ex. XIV, No. 6, p. 255. 


on the ring which was the point of contact of 
ring and table has moved up to T2 and the point 
on the cylinder which was, in the central position, 
coincident with 7; has now reached 73. Since 
there has been no slipping, the straight line TT, 
equals the arc TJ, on the shell and the arc 
QT: on the shell equals the arc QT; on the 
cylinder. 

Denote the angles TOT, TOQ by a and 8, re- 
spectively. Draw PS a vertical through P. The 
angle QPS also equals @. Denote the angle T;PS 
by ¢. The radius of the cylinder which was 
vertical in the central position is now inclined at 
angle @ to the vertical. The rotation of this 
radius is a measure of the rotation of the cylinder. 

We have more than one degree of freedom, 
and the simple method used earlier by me to get 
the frequency of vibration will not suffice here. 
Lagrange’s equations for “‘double freedom”’ are 
more convenient. These are® 


d oT oT 


dt 8g. aq: 


OV 


’ 


0q1 
aor or 


dt 0q2 —0q2 


0g2 


where q; and qe are the two variables, and T and 
V are the kinetic energy and potential energy, 
respectively, of the system. 

In Fig. 1 we have designated three angles 
a, 0, ¢, but a geometric relation between these 
serves to eliminate ¢ leaving the two variables 
a and @ which will take the place of qi and qe in 
the above statement of Lagrange’s equations. 


3H. Lamb, Dynamics (Cambridge University Press, 
Cambridge, 1923), 104-106, p. 317. See also other texts. 
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} Fic. 1. A cylinder rolling and rocking within a thin 
circular band or ring which is itself rolling and rocking on 
a horizontal table. At any instant the cylinder and band 
are rolling in opposite directions; the figure portrays an 
instant when the band is moving to the right and the 
cylinder to the left as indicated by the arrows. The dotted 
circles indicate the central positions of the band and 
cylinder. The insert on the right is the velocity diagram 
for the center of the cylinder. 


The geometric relation between a, 0, @ is ob- 
tained from the fact that since there is no slipping 
the length of the arc QT2 whose radius is R 
equals the length of the arc QT; whose radius 
is r. Therefore 


R@ —«a) a r(6+¢), 


$= {(R—r)@—Ra}/r. (3) 


To calculate the kinetic energy of the cylinder 
we need the velocity of its center. Consider the 
velocity diagram of P shown as an insert to 
Fig. 1. Relative to the center of the shell, the 
velocity of the center of the cylinder is PiV; 
=(R-—r)6 at right angles to OP and therefore 
inclined at the angle @ to the horizontal. The 
velocity of the center of the shell is Pi: V2=Ra 
horizontally. Therefore the velocity of the center 
of the cylinder with respect to the table-top 
(assumed to be at rest) is the vector sum P,V; 
of Pi V; and P,V2, and by geometry 


(P; V3)? = (P,; Vi)?+ (Pi V2)? —2P, Vi ‘Py V2 cosé 
= {(R —r)6}?+ (Ra)? 
—2{(R—r)6}(Ra&) cosé. 
If we limit the vibrations to those of very small 


amplitude, we may replace cos@ by unity, whence 
we get 


(velocity of center of cylinder)? 
={(R—r)6—Ra}*. (4) 


We now build up the Lagrangian equations: 


The kinetic energy of the shell, which is the 
energy of rotation about T may be written 


Tu =3M(K°+R’)e, 


where K =radius of gyration of the shell about 
its central axis. Since K=R for a thin cylindric 
shell, this reduces to Ty = MR?c?. 

The kinetic energy of the cylinder T,, equals 
the sum of the energy of rotation about its center 
and the energy of translation of the whole con- 
sidered concentrated at its center; Therefore 


Tm = 4mk?¢?+-3m{(R—r)6 —Ra}?, by Eq. (4), 
= 4m(1+k?/r’){(R—r)6 —Ra}?, by Eq. (3). 


Therefore the whole kinetic energy T, being 
equal to Ty+-T», is given by 


T = MR°a2+-3m(1+k2/r2) {(R—r)6 —Ra}?. (5) 


The potential energy V is given by Vu+Vm, 
in which Vy does not change and V,, is the 
weight of the cylinder times the change of 
elevation of its center. Therefore 

V= —mg(R-—r) cos#+ constant. (6) 


We now form the Lagrangian equations for a: 


aT 
— =2MRa+-m(1+k/r) 


0a 


xX {(R—1r)6—Ra}(—R), 
d oT 


— —=2MR%i—m(1+k/r)R{(R—-1)b—Ra} 
dt 0« 


= {2MR?+m(1+k?/r’)R?}a 
—m(1+k?/r*)R(R—r)6, 


Therefore the Lagrangian equation for a is 


(2MR?+-m(1+k/P) Ry a 
—m(1+k?/r?)R(R—r)§=0, (7) 


which enables us to express a in terms of 8. . 
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We now form the Lagrangian equations for 6: 


T 
“= 0+ m(1+H/1)((R—1)0-Ray(R-N), 


, 1+k?/r?){(R §6—Ra}(R 
ee +k?/r*){(R—r)6—Ra}(R—r), 


oT 


06 


, 


OV 
[nn sind =mg(R—r)6, 


as @ is limited to very small values. 
Therefore, the Lagrangian equation for @ is 


—m(1+k?/r?)R(R—r)a+m(1+k?/r’) 
X(R—r)*6+mg(R—r)0=0. (8) 


By use of Eq. (7) we substitute for a in Eq. (8) 
thus getting an equation involving 0, 6, and 
constants. It is 


m(1+k?/r’?)R(R—r) 
2MR?+m(1+k?/r*)R? 
+m(1+k?/r?)(R—r)*§+mg(R—r)0=0, 
which reduces to 
(R=2) 2M(1+k?/r?) 
2M+m(1+k?/r’) 


—m(1 + /A\RR-1)| 


+20=0, (8) 


whence we see that the motion in 6 is simple 
harmonic and that the length of the equivalent 


TABLE I. Details of experiments with models of form 
(Fig. 2). The outer body: A thin brass ring (or band or 
shell) of mass M=301g, radius R (mean) =10.05 cm, 
thickness 0.2 cm, length 2.9 cm. In all cases and in many 
others not detailed here there is good agreement between 
theory and practice. 


The inner Thin-walled Solid brass Solid wood 
body brass tube cylinder cylinder 
> (Fig. 2A) (Fig. 2B) (Fig. 2C) 


Steel ball- 
bearing 
(Fig. 2D) 


Mass m (g) 
Radius r (cm) 
Thickness (cm) 
Length (cm) 

L cale from 
Eqs. (10) to 
(12) (cm) 

L 


Period, calc 


1400 


92 
2.17 (mean) 2.51 
0.1 eee 


(sec) 

No. of vibra- 
tions timed, 
and repeti- 
tions thereof 

Period, mean 
obs (sec) 
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Fic. 2. Four —— of a rolling and rocking brass 
band containing a rolling and rocking body. In (A) the 
inner body is a thin-walled metal cylinder, in (B) a solid 
metal cylinder, in (C) a larger solid wooden cylinder and 
in (D) a steel sphere (a ball-bearing). Arrows indicate the 
direction of motion for one half the period; in the other 
half the directions are reversed. 


simple pendulum is given by 
2M(1+k?/r’) 


L= (R-r)—— ’ 
2M+m(1+k?/r?) 


(9) 
which is Lamb’s expression. The period of the 
vibration =2zx(L/g)}. 


Corollaries 
(1) An alternative expression for L is 
L=2(R-—r)—————. 
2/(1+-k?/r?)-+m/M 
(2) If the cylinder is hollow and thin walled, 


k? =r’; therefore 


L=(R-r)————_ 
2M+2m 


1 
= 2(R-—r)——_.. 
( im/M 


(3) If the cylinder is solid, k? = 4r?; therefore 


=2(R—r) 


(10) 
M+m 


L=(R—-r) 2(R—-r) 


2M+3m/2- 4M+3m' 


1 
R-r)—————_. 
2¢ "/34m/M (a) 


M 
= 2(R—r)————-= 
4M/3-+m 


(4) If the cylinder is replaced by a sphere of 
radius 7, k? = 2r?/5, and therefore 


2M X7/5 


L=(R-—r)——————_= 
2M+mxX7/5 


2(R -r)————- 
10M+7m 


=2(R—r) (12) 


10/7-+-m/M 
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Sometimes one form is useful for calculation, 
sometimes another. In CGS units the period, 
calculated from the value of LZ given by the 
masses and dimensions (Eqs. (9) to (13)), is, in 
Toronto, very nearly equal to }-L}. The error 
=0.3 percent. 

Equations (9)—(12) have been checked with 
many rings (thin and thick) and cylinders (hollow 
and solid). I quote in Table I a few typical 
results (see Fig. 2). 


Il. THE THIN OUTER RING IS REPLACED BY A 
THICK RING 

The proof given above for the thin-walled 

shell may be modified to suit the case of a thick 

cylindric ring of inner and outer radii R; and Re, 

respectively. Figure 3 represents this case. The 


L=[R:-r] 


which, of course, reduces to Eq. (9) if we let 
Ri =R.=R. 

Experimental check: Table II, first line, shows 
the results of a test in which the outer body is a 
thick brass ring of rectangular cross section, hav- 
ing a radius of gyration given by K?= 148.5 cm?. 
For the inner body a solid brass cylinder was used 
giving for the length of the equivalent simple 
pendulum L=8.73 cm, and therefore the calcu- 
lated period equaled 0.59 sec. Observations 


Fic. 3. A cylinder rolling and rocking within a thick 
circular ring which is itself rolling and rocking on a hori- 
zontal table. At any instant the cylinder and the ring are 
rolling in opposite directions; the figure portrays an instant 
when the ring is moving to the right and the cylinder to 
the left as indicated by the arrows. The insert on the right 
is the velocity diagram for the center of the cylinder. 


M(A+k?/r?) (K2+ Re?) +m(1+k2/r) (#R2/P?+ Re) —m(RRi/r?+ Ro)? 
M(K?+R2)+m(#R?/r+R2) | 
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letters have the same significance as in Fig. 1 
except that 7, is now on the inner surface of the 
thick ring. In the central position this point was 
the point of contact of cylinder and ring and was 
vertically above T;. In this case 


rp = (Ri —7r)6 —Ria, 
(velocity of center of cylinder)? 
= {(Ri —7r)6 —Rod}?, 
Tu =3M(K°?+R?2)e, 
and 


Vm = —mg(R; —r) cos#+ constant, 


where K is the radius of gyration of the thick 
ring about its axis. 

With these substitutions the work proceeds 
as above with the result that 


(14) 


showed 10 vibrations in 6.0, 6.0, 6.0, and 6.0 
sec, whence the observed period is 0.60 sec. 
The agreement is good. 

The rest of Table II refers to the experiments 
of the subsequent sections. 


Il. THE COMPLETE RING IS REPLACED BY A 
THIN-WALLED SEMICIRCLE 


Figure 4 shows the semicircle and the cylinder 
in an off-the-center position. The letters are as 
in Fig. 1 with the addition of G, the center of 
gravity of the semicircle. Point G is at a distance 
2R/x from the center of the circle along its 


TABLE II. Details of experiments with cylinders rocking 
in rings and semicircles. In all cases the objects were of 
brass. Observed periods of vibration were obtained by 
successive measurements with a stopwatch of the time 
required for 5, 10, or 20 vibrations. 


Motion of the two bodies 
in opposite directions in the same direction 

Period Period Period Period 
calc obs L calc ebs 
(sec) (sec) (see"!) (cm) (sec) (sec) 


Motion of the two bodies 


Outer Inner n L 
body* body* (sec™!) (cm) 


A 3 
B 


3.12 101. 
4.05 


8.7 
3.50 
2.40 
3.34 


59.8 


®A Thick ring: M—2508 g, Ri—10.65 cm, R2—13.55 cm, thickness 
1.35 cm. B Thin-walled semicircle: M—177.5 g, R(mean)—15.35 cm, 
thickness—3 mm, width—1.35 cm. C Thin-walled semicircle: M—115 g, 
Ri—4.90 cm, thickness—1.8 mm, width 5.08 cm. D Thick-walled semi- 
circle: one-half of ring A, M—1i253 g. E Solid oe m—1400 g, 
r—2.51cm, length 8.3cm. F Solid cylinder: m—328 g, r—2.00cm, 
length 3.05 cm. G Thin-walled tubing: m—41 g, r(outside)—1.20 cm, 
thickness—0.5 mm, length 13.6 cm. 
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Fic. 4. A cylinder rolling and rocking within a thin-walled 
semicircular band which is itself rolling and rocking on a 
horizontal table. The figure portrays an instant when the 
cylinder and band are rolling in opposite directions as 
indicated by the arrows. The dotted semicircle indicates 
the central position of the band. The insert on the right 
is the velocity diagram for the center of the cylinder. 


middle axis. Call this distance },. In the central 
position the moment of inertia of the semicircle 
about the point 7; is equal to 


M(R? —b2+(R —b;)?) = M(2R? —2Rd)). 
This is also, for small observations very nearly, 
the moment of inertia about T. In motion the 
center of gravity of the semicircle rises. In the 


a-position the elevation is b;(1 —cosa). Therefore 
the expressions for the energies are 


T=3M(2R?—2Rb;)e2 


+4m(1-+#/r){(R—r)6 — Ray’, 
V= — Mgb; cosa 


—mg(R—r) cos#+constant, 


(14a) 


and the Lagrangian equation for a is 
{M(2R(R —b1)) +m(1+k?/r*) R234 
—m(1+k?/r?)R(R —r)6+ Mgbia=0, 
and for @ is 
—m(1+k?/r?)R(R—-r)a 
+m(1+k?/r?)(R —r)*6+mg(R—r)0=0. (16) 
The presence of the term Mgbia in Eq. (15) 
blocks our previous method of expressing a in 


terms of 6 and substituting in Eq. (16). We now 
have to use the method for small oscillations 


about equilibrium given in Lamb’s Dynamics, 
Sec. 109. Here he writes our Lagrangian equa- 
Aa+Hé+aa+ho=0, 
Hai+Béb+ha+b0=0. 
a=Ccos(nt+e) and 6=AC cos(nt+e), (18) 
where C, \, and ¢€ are constants. The equations 


(15) 


tions in the form 
} (17) 
To solve Eqs. (17) assume 
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} (19) 


(20) 


are satisfied provided 


(n?A —a)+(n?H —h)X\=0, 
(n?H —h)+(n?B —b)\=0. 


Eliminating \ we have 


nA—a n*H—-h 
vwH—-h n?B—b 


a 


which is a quadratic in n?. The formulas (18) 
therefore constitute a solution of the Eqs. (17), 
provided m satisfies Eq. (20) and provided the 
coefficient \ be determined in accordance with 
either of the Eqs. (19). Obviously the length of 
the equivalent simple pendulum is given by 
L=g/n’, and the period = 22/n. 

This method of solution was applied to the 
case of the rocking of a large thin brass semicircle 
B carrying a solid brass cylinder F. When the 
two bodies are rocking in opposite directions 
observations give, for the observed period, 
0.38 sec. 

The values of Lamb’s coefficients are 


A=2MR(R—b:)+m(1+k?/r?) R? 
= 14.6410‘ g-cm?, 
H= —m(1+k?/r’)R(R —r) = —10.08 X 104 g-cm?, 
B=m(1+k?/r?)(R —r)? =8.77 X10 g-cm?, 
a= Mgb,=170.1 X10! dyne-cm, 
b=mg(R —r) = 429.3 X10‘ dyne-cm, 
h=0. 


The quadratic equation in n? is therefore 


(14.64? — 170.1) (8.77? —429.3) = (10.08? —0)?, 


whence 


#22 16.7 or 3.12. 


Fic. 5. A cylinder rolling and rocking within a thin- 
walled semicircular band which is itself rolling and rocking 
on a horizontal table. The figure portrays an instant in the 
motion when cylinder and band are rolling in the same 
direction as indicated by the arrows. The dotted semicircle 
shows the central position of the band. The insert on the 
right is a velocity diagram for the center of the cylinder 
at the instant portrayed. 
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Taking the first value of m, the calculated 
period, =27/16.7=0.376=0.38 sec gives good 
agreement with the observed value. 

But what does the value »=3.12 signify? Let 
us see what happens if instead of the cylinder 
and the semicircle rocking in opposite directions 
they rock in the same direction. In this case 
(see Fig. 5) the following geometrical equations 
hold instead of Eqs. (3) and (4). 

For the equality of arc lengths, QT: and QT; 
R(a+6)=r(6+¢9), 


whence 
ro=(R—-r)é+Ra, (21) 


and for the velocity of the center of the cylinder 
(see insert to Fig. 5). 


PiV3 =P,)V?+P1V2+2P1Vi:PiV2 cos@ 
={(R—r)6+Ra}*, (22) 
when @ is very small. 


For this case the Lagrangian equations become 


{2MR(R —b;)+-m(1+k?/r)R*} a4 
+m(1+k?/r?)R(R —r)6+ Mgbia=0, 
m(1+k?/r?)R(R—r)a 
+m(1+k?/r?)(R —r)?6+-mg(R —r)o=0. 


The only difference from those of Eqs. (15), (16), 
and (17) is that Lamb’s coefficient H is now 
positive. 

Therefore, it follows that the quadratic in n? 
is now the same as before and that the »=3.12 
corresponds to the present case of the two bodies 
rocking in the same direction at the same time. 
From =3.12 the calculated period is 27/3.12 


(23) 


Fic. 6. A cylinder rolling and rocking within a thick 
walled semicircular —— is itself rolling and rocking 
on a horizontal table. This figure portrays the case when 
the cylinder and semicircle are, at the same instant, moving 
in opposite directions, as indicated by the arrows. The 
insert on the right is the velocity diagram for the center of 
the cylinder. 
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A fair agreement between experiment and 
theory was observed when the last quoted bodies, 
rocking in the same direction, gave an observed 
period 2.06 sec. Of course, the semicircle is not 
theoretically thin. 

On a second trial a smaller, wider semicircle 
was taken, and the cylinder was a piece of thin 
brass tubing. These objects are designated as 
C and G in Table II, where dimensions and 
masses are tabulated. 

Lamb’s coefficients are A =4137 g-cm?, H= 
—1558 g-cm?, B=1181 g-cm’, a=366g dyne-cm, 
b=156g dyne-cm, and h=0. In coefficients a and 
b the acceleration of gravity is represented by g. 
The quadratic equation in n? is, therefore, 


(n?4137 —366g) (m*1181 —156g) = {n?( —1558)}? 
or 


n*X4.89 X 10° —n? X 10.80 X 108 
+5.48 X 10!° = n* X 2.43 X 108, 


whence »=19.5 and 7.65 and the calculated 
periods 0.32 sec and 0.82 sec. The observations 
were (1) when rocking in the same direction, an 
observed period of 0.80 sec, and (2) when rock- 
ing in opposite directions; an observed period 
of 0.33 sec. The agreement is good. 


IV. THE THIN-WALLED SEMICIRCLE IS REPLACED 
BY A THICK-WALLED SEMICIRCLE 


But no ring or semicircle is as thin as a 
theoretical thin ring or semicircle. If it were so 
it would not keep its shape. So let us now take 
a thick semicircle and see how the equations are 
modified. 

Let the internal and external radii by R; and 
R, (Fig. 6). Let the center of gravity G of the 
semicircle be at a distance b; from the center of 
curvature O. By geometry 


by = (4/32) (R —R;*)/(R? —R,’). 


(1) Rocking in Opposite Directions 


As before (see Eq. 13 for the thick circle) we 
have 


ro= (Ri —1r)d —Ria 
and (velocity of center of cylinder)? = 


{(Ri —r)6 —Ree}?; 
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also 


Tu =}4M{K?+(R2—5)*}, 
Sel 


Ta = imit| 
r 


+$m{(Ri—r)6—Raa}?, 


and 


V= —Mgb, cosa —mg(R; —7r) cosé+ constant, 
and the Lagrangian equations are for a 


[M{K*+ (R2—)b:)"} +-m(k?/rPR2+R-) Ja 
—m(R, —r)(k?/r?Ri + R2)6+ Mgbh,a=0 
and for 6 © 
—m(R, —1)(k?/r?Ri+ Ro) a+m(Ri —71)? 
X (1+k?/r?)6+mg(Ri —r)0=0. 


Experiment: The thick semicircle D used was 
one-half of the thick ring A used in Sec. II. The 
cylinder was also the same as in Sec. II: 


b,=7.74 cm (calc), 
K?=}(R?’+R,’) —b;? = 88.6 cm?, 
R,-r=8.14cm, (Ri —r)?=66.3 cm?, 
R2—b,=5.81cm, (R2—b;)?=33.8 cm?. 


(24) 


The coefficients in Lamb’s equations (17) are 


A = M{K?+ (R2—b;)?} +m(k?/rP’R?+ R2*) 
=4.90 X 10° g-cm?, 
H= —m(Ri—17)(k?/r?Ri+Re2) 
= —2.15X10° g-cm?, 
B=m(R,—1)?(1+k?/r?) =1.39X 105 g-cm?, 
a= Mgb,=95.0X 10° dyne-cm, 
b=mg(Ri —7r) = 111.7 X 10° dyne-cm, 
h=0. 


The equation for n? is 
(4.90? —95.0) (1.39m? —111.7) = ( —2.15n? —0)?, 


whence ~=17.1 or 4.05. The calculated periods 
are 0.367 and 1.55 sec, respectively. 

Observations: For semicircle and cylinder rock- 
ing and moving in opposite directions at the same 
time the observed period was 0.37 sec. 


Fic. 7. A cylinder rolling and rocking within a thick- 
walled semicircular ring which is itself rolling and rocking 
on a horizontal table. The figure portrays the case when, 
at the same instant, the cylinder and the semicircle are 
moving in the same direction as indicated by the arrows. 
The insert on the right is the velocity diagram for the 
center of the cylinder. 


(2) Rocking in the Same Direction 


For the bodies rocking and moving in the same 
direction at the same time (Fig. 7) the equations 
follow the same pattern as in the case of the thin 
semicircle. The mean observed period was 
1.55 sec. 

When moving in the same direction it was 
observed that a jiggling motion corresponding 
to the mode of greater frequency tends to grow 
and superpose itself on the mode of greater 
period. The existence of two modes of rocking 
with the semicircle makes its problem a little 
more interesting than that of the complete ring. 
Both, however, are good fun. The lulling effect 
of the rocking cradle on the loosely packed-in 
baby provides another example of happiness 
easily obtained in this way. These experiments 
give a student good exercise in the use of the 
Lagrangian equations in a problem of two degrees 
of freedom. The materials for experimentation 
are very simple and are available in any good 
workshop. A good stopwatch is, however, 
necessary. I use one where the complete revolu- 
tion of the seconds hand takes 10 seconds. 

I have much pleasure in thanking a colleague, 
Dr. D. G. Ivey, for directing my attention to 


Lamb’s problem and for help in the first part of 
the theory. 












Geometrical Solution for a Nuclear Moment in a Magnetic Field* 


























































































rate of change of angular momentum, becomes! 
dy/dt=g(uXH), (1) 


where g, the gyromagnetic ratio, is the ratio of 
the magnetic moment to the nuclear angular 
momentum. If H is constant, the solution of 
Eq. (1) is a steady precession about the direction 
of H with uw making a constant angle to this 
direction.? The (circular) frequency of precession 
is the Larmor frequency 


















































or=g|H|. (2) 








For simplicity, we shall assume that g is positive, 
so that the projection of u on a plane perpendicu- 
lar to H rotates clockwise when viewed in a 
direction opposite to that of H. If g is negative, 
the sense of precession is reversed. 

In resonance experiments, H is made to de- 


























* Supported in part by the U. S. Office of Naval Research. 
1W. J. Archibald, Am. J. Phys. 20, 368 (1952). 
2G. E. Pake, Am. J. Phys. 18, 438 (1950). 
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Thus, in the rotating system, the behavior of u 
is determined by an equation of the same form 
as Eq. (1) with an “‘effective”’ field in this system 
given by H+(a,./g). The advantage of this for- 
mulation is that by making a judicious choice of 
the rotating system we can give the effective field 
such a simple form that the solution of Eq. (5) is 
very easy to obtain. 

Our first transformation is to a primed system 
whose x’ axis lies along H; and rotates with it as 
illustrated in Fig. 1(a). This means that we set 
w’ = —wk, so that we find from Eq. (5) the field 
in this system to have the components 


H, =H, H, =0, H,, = Ho —(w/g). (6) 
If we introduce the abbreviations 
A=gHy—w=ao-w, wi1=ghMh,, 


3G. Joos, Theoretical Physics (Blackie and Son, Limited, 
London, 1934), p. 144; H. Goldstein, Classical Mechanics 
(Addison-Wesley Press, Inc., Cambridge, Massachusetts, 
1950), p. 132. 


RoaLtp K. WANGSNESS ™ 
U. S. Naval Ordnance Laboratory, White Oak, Maryland 
(Received October 3, 1952) 
TI 
The general solution of the classical equation of motion du/dt=g(uXH), when H consists 
of a constant field plus a perpendicular rotating component can easily be obtained after trans- fie 
forming the problem to a properly chosen rotating coordinate system. The value of this method m: 
for teaching purposes is discussed. we 
N a recent article in this journal, Archibald! pend on the time and generally consists of a con- ar 
pointed out that the behavior of an otherwise stant component H> and a circularly polarized 
free nuclear magnetic moment y in an external component rotating in a plane perpendicular to 
magnetic field H provides a useful source of class- Hp» and in the same sense as the Larmor preces- 
room problems, particularly for a class in quan- sion. The components of H can therefore be by 
tum mechanics. The solution of the classical written as 
equations of motion for the fields of greatest s t 
practical importance is also of interest in that it toed, Lys - Rte, yet fi 
proves to be a good nontrivial example of the use If we let Du/Dt be the rate of change of y as 
of the general transformation of time derivatives _ seen in a coordinate system rotating with angular 
into rotating coordinate systems which students velocity w, with respect to the laboratory system T 
otherwise tend to feel is useful only for the deriva- in which Eq. (3) is applicable, then we have® 
tion of Euler’s equations for the motion of a rigid c 
body. As we shall see, this elegant method, which dy/dt= (Dy/Dt) + (a. Xw). (4) a 
was first suggested by H. C. Torrey, enables us to Combining Eqs. (1) and (4), we find that ¢ 
write down the general solution of the vector t 
equation of motion almost at a glance. Dy/Dt = g{uX (H+ (@./g) J}. (5) 3 
The classical equation of motion, torque equals c 
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the components in the primed system can be 
written as 


Hy=o/g, Hy=0, Hy =A/g. 


These are indicated in Fig. 1(b). 

We also see from Fig. 1(b) that the resultant 
field in the primed system lies in the x’z’ plane, 
makes an angle @ with the z’ axis, and has the 
constant magnitude a/g, where 


a= [A?+w? |}, 
and 


cos#=A/a, sind=w;/a. 


Next, we transform to a doubly primed system 
by a simple rotation of axes through an angle @ 
about the y’ axis so that the 2” axis is parallel to 
the resultant field a/g. The components of the 
field in this system are then 


Ay =H, =0, Hy =a/g. (7) 


This system is also shown in Fig. 1(b). 

Since the field in the doubly primed system is 
constant, the motion of u is a simple precession 
about the 2’’ axis while making a constant angle 


¢ with this axis. According to Eqs. (2) and (7), 
the frequency is a. Therefore, if we take into 
account the sense of precession as previously 
discussed, the general solution in this system is 


Mz =p sing cos(at —e) 
My = —p sing sin(at —e) (8) 
Mz =p COS@=const, 


where »=|w| and ¢ and e are constants giving 
the orientation of yw at t=0. 

We need not have stopped with this last trans- 
formation, of course, but could have made one 
more to a system rotating about 2” with angular 
velocity w’’’ = —ak’’. The resultant field is then 
zero by Eq. (5), so that w is a constant in this 
new system. Transforming back to x", y’’, 2’ 
axes, we then get Eqs. (8). 

To obtain the general solution in the labora- 
tory coordinates, we see that we need only use 
the equations below which are easily obtained by 
inspection of Fig. 1: 


Mz’ =z CosO+ py," sind 
My’ = My’’ 
Mz) = —ymz sind+y,- cos 


and 


Mz= Mz COSwh+ py sinwt 
My = —z Sinwt+ py Coswt 
Uz Hz! 


Finally, combining these equations with Eq. (8), 
we have 


uz=pu{Lcos¢ sind+sing cos# cos(at —e) ] 
Xcoswt —sing sin(at —e) sinwt} 

uy = u{Lcos¢ sind+sing cos# cos(at —e) ] (9) 
Xsinwt+sing sin (at —e) coswt} 

Lz = uLcos@ cosé —sing sin@ cos(at —e) ]. 


Fic. 1. The coordinate systems used to simplify the equa- 
tions of motion. The heavy arrows are the magnetic field 
components in the different systems. 


Equation (9) is actually the desired general 
solution since it contains two constants which are 
sufficient to specify the initial orientation of yu 
and, as can easily be verified by substitution, 
satisfies Eq. (1) with the field components as 
given by Eq. (3). 

The general motion is thus seen to be a pre- 
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cession of frequency w with a superimposed nuta- 
tion of frequency a. 

Explicit expressions for two cases of practical 
importance, which were also discussed by Archi- 
bald, can easily be obtained from Eq. (9). 

First, suppose that at ¢=0, uz=u, and yp, 
=py=0. We can see from Fig. 1(b), or verify 
from Eq. (9), that this determines the constants 
to be ¢=0@ and e= +7. Equations (9) give at once, 


Hz=p siné[cos#(1 —cosat) coswt+sinat sinwt |, 
Hy = —p siné[cosé(1 —cosat) sinwt —sinat coswt |, 
Ms =p (cos’8+sin’6 cosat). 


These agree with Archibald’s equations after one 
corrects his Eqs. (11) by dividing them by a and 
also noting that his A is the negative of ours. 
The usual equations which are used to discuss 
resonance experiments correspond to the special 


case 4,=const and therefore ¢=0, which leads to 


Hz=u(wi/a) coswt, py= —A(w:/a) sinet, 
u2=u(A/a), 


for simple precesssion without nutation. 

Lest this procedure seem too much like magic 
to the students, it will be instructive to have 
them start with the three simultaneous equations 
obtained from Eqs. (1) and (3), and make the 
analytical substitutions required by the succes- 
sive transformations. At each step, one can 
inspect the resulting equations and verify that 
they do have the form of Eq. (1) but with the 
various field components Eqs. (6) and (7). At the 
same time, one sees how each step simplifies the 
equations until the final transformation puts 
them into the usual harmonic oscillator form, 
which is immediately soluble. 


The Tape Recordings of Important Speeches 


BY 


Fk. U. Connon, K. K. Darrow, E. Fermi, J. C. SLATER 


Free tape recordings are available on loan to members 
of the Association who wish to use them for presentation 
to classes, seminars, or science clubs. The subjects are 
four of the six lectures presented as a Symposium on 
Physics Today at the Twentieth Anniversary meeting of 
the American Institute of Physics, Chicago, Illinois, on 
October 25, 1951, as follows: 


The atom. E. U. Conpon. 

Physics as science and art. K. K. Darrow. 
The nucleus. E. Fert. 

The solid state. J. C. SLATER. 


Recordings are now available in two forms: 


1. Double track. Each speech is complete on one reel, 
which must be played at 3} inches per sec. No re- 
winding is necessary. Track No. 1 is played to the 
end, then the reels are reversed (and turned over) to 
play track No. 2. This procedure rewinds the tape 
into its original condition. 


2. Single track. Two 7-inch reels are required for each 
address. They must be played at 7} inches per sec, 
and, of course, rewound to restore them to their 
original condition. 


Requests for recordings should be made to THomas H. 
Oscoop, Editor, American Journal of Physics, Michigan 
State College, East Lansing, Michigan. To allow some 
flexibility in scheduling, a request should specify a pre- 
ferred date and two alternative acceptable dates, as well 
as the type of recording (single or double track, or either). 

No charge is made for the loan of recordings. The bor- 
rower is expected, however, to pay the return postage to 
the Editor’s office or to any other address the Editor 
specifies. Transportation is to be by first class mail, special 
delivery, which will cost less than one dollar per reel at 
present rates. Reels will be mailed in convenient reusable 
boxes. 

Prospective borrowers may consult the published papers 
in Physics Today, November, 1951, January, 1952, and 
March, 1952. 
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Physics and Physiology in Diving Decompression 


HILBERT SCHENCK, JR. 
Pratt and Whitney Aircraft, East Hartford, Connecticut 
(Received September 12, 1952) 


The fundamental problem of not exceeding the critical saturation of nitrogen in a diver’s 
body requires both physics and physiology in its solution. The basic differential equations of 
nitrogen saturation and desaturation in a specific tissue are integrated over an idealized dive. 
Constants determined through physiological research can be substituted and a time-depth 
ascent curve for any tissue established. A stepwise method is outlined whereby critical tissues 
at each portion of the ascent can fix the proper ascent curve for the entire body. The advantage 
of pure oxygen breathing is demonstrated numerically and the failure of helium-oxygen mixtures 
to shorten ascents is explained by use of the derived equations. 


INTRODUCTION 


LTHOUGH the general subject of deep 
water diving has already received attention 
in these pages,’ a fuller analysis of nitrogen satu- 
ration in organic tissue may prove of interest to 
physicists. This appears to be primarily a physio- 
logical problem, but the solution given here is so 
closely related in method of attack to certain 
physical situations that the boundaries between 
the two fields become difficult to define. It should 
be noted, however, that whereas a physical theory 
can be proved a hundred times over, a physio- 
logical theory is never rigorously applicable to all 
individuals, but only to the few “‘average’’ indi- 
viduals in the group. 

When a man breathes air under pressures 
greater than atmospheric, his body tissues 
gradually take up additional nitrogen in simple 
solution, the amount depending on the partial 
pressure of the inhaled nitrogen and the time 
spent in the high pressure region. The oxygen 
component of air is not dissolved in tissue to any 
extent since it combines chemically within the 
body. Following prolonged exposure to high 
pressures, a sudden return to atmospheric condi- 
tions results in supercritical nitrogen saturation 
in the diver’s body and the formation of nitrogen 
bubbles. The bubbles produce a variety of 
symptoms including itching, pain in the joints, 
and even paralysis and death. 

Although the cause of this disorder (called 
the bends or caisson disease) was known as 
early as 1854, and an investigation accorded it 
by Paul Bert in 1878, it was not until 1908 that 


1L. E. Dodd, Am. J. Phys. 8, 181 (1940). 


Dr. J. S. Haldane produced the first Decom- 
pression Tables, regulating the ascent of a diver 
from the sea bottom.? So excellent was Haldane’s 
work that the tables currently in use by the 
United States Navy? are almost identical to those 
original ones worked out almost fifty years ago. 


THE DESCENT 


The manner in which a specific body tissue 
becomes saturated with nitrogen can best be 
illustrated by the simple model in Fig. 1. Blood 
entering the lung alveoli is almost instantly 
saturated with nitrogen at the diver’s external 
partial nitrogen pressure. As the blood flows 
through a tissue saturated with nitrogen at some 
lesser pressure, some of the gas will diffuse into 
the tissue. The nitrogen pressure rise in the tissue 
will be dependant on the pressure difference 
(P.—P), the area of contact between blood and 
tissue A, and the tissue resistance to nitrogen 
transfer R. The process can be described by the 


Fic. 1. Model of nitrogen satu- 
ration. Blood flows between 
blocks of tissue containing dis- 
solved nitrogen at pressure P,.. 
Nitrogen seeps through mem- 
branes X—X and Y-—Y into 
tissue with nitrogen pressure P. 
Rate of seepage depends on area 
of membrane and resistance to 
nitrogen movement as well as 
pressure differential. 

a 


2 See J. S. Haldane, Respiration (Yale University Press, 
New Haven, 1922), pp. 334-357. 

?Bureau of Ships Manual, Chapter 94, Section II, 
“Diving” (1951). 
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Fic, 2. Safe ascent curves for four representative tissues 
with indicated Ty values for one-hour dive at 200 feet. 
Dotted line is sea-level pressure and stepped curve is 
standard decompression schedule. 


relation, 
(1) 


where T is the time in minutes. (The similarity 
between Eq. (1) and the conduction rate equa- 
tion for heat transfer need not be belabored.) 

If we postulate the simple case of a diver 
whose tissues are fully saturated with nitrogen 
at one atmosphere (P,:), suddenly being sub- 
jected to a higher, ‘‘sea bottom”’ pressure (P.), 
Eq. (1) may be integrated to give 


oe] 


eATIR 


(2) 


P=P.-| 


which expresses the tissue nitrogen pressure as 
a function of T, the time spent ‘‘on the bottom.” 

It is convenient to express A/R in terms of the 
following equation 


[: aS & . ™ _ 


which states that when a given tissue, com- 
pletely free of nitrogen, is subjected to a sudden 
pressure increase (P,), it will require 7; minutes 
to become half-saturated with nitrogen. Inte- 
grating and rearranging Eq. (3) results in 


A/R=log.2 /T;, (4) 


so that the half-saturation time T, of a tissue com- 


JR. 


pletely specifies its behavior providing 7; has 
been determined in accordance with Eq. (3). 
(Note the similarity between 7; as regards satu- 
ration and the half-life of an element as regards 
radioactive disintegration.) 


THE ASCENT 


Having remained a certain length of time on 
the ocean bottom, the diver must be brought to 
the surface in such a manner that the excess 
nitrogen held in the tissues will not form bubbles 
but will diffuse smoothly into the blood and 
be exhaled from the lungs. Following the reason- 
ing in the previous section it is easily seen that 
the desaturation equation is 

—dP'/dT'=A/R(P’ -P.’), (5) 
where P’ is the nitrogen partial pressure in a 
specific tissue, dT’ a time increment during the 
ascent, and P,’ the external nitrogen partial 
pressure at the diver’s depth. (In all formulas, the 
prime symbol refers to ascent conditions.) We 
desire to regulate the ascent such that (P’ —P,’) 
is a safe maximum, and the question arises; how 
should P’ and P,’ be related? 

Through various experiments on men and 
animals Haldane? concluded that a tissue exposed 
to air at a given total pressure could immediately 
be placed in air at slightly less than half the 
original pressure without bubble formation. 
When air of atmospheric composition is used 
throughout the dive this results in the relation, 


P’=nP,/, (6) 


where n for air is slightly greater than 2. Com- 
bining Eqs. (5) and (6) and integrating over the 
limits P to P’ and 0 to T’ yields the desaturation 
curve for any specific tissue, 


- n a, nP,’ , 
--() glee(S I @) 


where P is found from Eq. (2). 

Although the previous equations involve 
nitrogen partial pressures, all P’s can equally 
well represent total air pressures, since the 
nitrogen component remains a constant per- 
centage of the total throughout the dive. Thus 
P.’ will be the absolute pressure surrounding the 
diver, and its variation with time is the decom- 
pression curve for the tissue considered. Now 
Eqs. (2) and (4) can be substituted in Eq. (7) to 
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give the final result, 


rat = 
| tos.( P. a —log.(nP.) | (8) 


Equation (8) will give the safe desaturation time 
T’ from a bottom total pressure P, to any lesser 
pressure P,’ for a tissue characterised by its half- 
saturation time T,; providing, of course, that the 
T’, P.’ schedule specified by Eq. (8) is followed 
throughout. Equation (8) assumes that descent 
is instantaneous as is the initial stage of the 
ascent when T” is negative. 

Haldane? found that the “slowest”’ tissue in the 
body has a 7; of seventy-five minutes and all 
other tissues had varying half-saturation times 
between this maximum value and a few minutes. 

As an example, consider a sixty-minute dive in 
200 feet of water. All pressures can be substi- 
tuted as depths in sea water (33 feet equals one 
atmosphere) and all pressures, as noted pre- 
viously, can be total pressures. Then; P,= 233 ft, 
Pa=33 ft, T=sixty minutes, and n=2.00 
(normal air). Time-depth ascent curves are 
plotted in Fig. 2 for tissues having seventy-five, 
forty, twenty, and ten-minute 7; values. Each 
curve is the safe decompression schedule for the 
individual tissue, and the stepped curve is the 
stage decompression schedule recommended by 
the United States Navy? for this dive. 

Obviously, one tissue cannot specify the 
entire ascent. Obeying the seventy-five minute 
“‘slow’”’ tissue will result in an initial decrease 
that is dangerously large for the ‘‘faster’’ tissues 
which are more highly saturated. Following the 
‘faster’ tissue curves will provide insufficient 
time for the ‘‘slow’’ tissues to desaturate and 
trouble will occur at the surface. The over-all 
decompression curve must be calculated basing 
the early part of the ascent on the “‘fast’’ tissues 
and the portions nearer the surface on the 
“‘slow’’ tissues. Since the ‘‘slow’’ tissues are not 
under a large pressure differential at the begin- 
ning of the ascent, they lose little nitrogen until 
the ‘‘fast’”’ tissues have been partially cleared 
and the diver has ascended to a lower pressure. 
The over-all time from bottom to surface is thus 
longer than that predicted for a specific tissue 
and Eq. (8) is not suitable for a rigorous analysis 
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of the ascent. The ascent differential equation (5) 
must be solved in a tedious series of stepwise 
calculations in which several representative 
tissues are carried along. Equation (8) will yield 
the full curve in one special case, namely the 
condition of complete saturation of all tissues. 
For this situation, To, T,=seventy-five 
minutes so that P=P,. (Needless to say the 
actual, full curve is smooth and not in “‘stages.”’ 
The problem comes in training the diver’s tender 
to pull up his man at a logarithmic rate.) 


PURE OXYGEN BREATHING 


Long periods of time spent in deep water 
require prohibitively long ascents. To alleviate 
this condition, Sir Robert Davis devised a 
chamber to be held in sixty feet of water into 
which an ascending diver climbed, undressed, and 
relaxed while breathing pure oxygen. Bubble 
formation is dependent on gross outside pressure 
regardless of gas composition while nitrogen 
escape depends on external nitrogen pressure 
alone. Thus breathing pure oxygen has the effect 
of reducing P,’ in Eq. (5) to zero, but has no 
effect whatsoever on bubble formation. Then 
Eq. (5) can be integrated for this case (defining 


A/R as in Eq. ve to 
p” 


where P’ is the tissue pressure when pure oxygen 
is first admitted, P”’ is the tissue nitrogen pres- 
sure at any later time T’’. (T’’=0 when P” =P’, 
the instant that normal air is changed to oxygen.) 
Unfortunately, pure oxygen is a poison at pres- 


T= (9) 
= 


"PURE O2 ASCENT 


FEET OF SEAWATER 


oO 200 
MINUTE S 


Fic. 3. Safe ascent curve for fully saturated seventy-five 
minute tissue from 300 ft. 33-ft line is sea level, and pure 
O: line is alternative ascent with pure oxygen administered, 
starting at 60 ft underwater. 
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sures of greater than two-and-a-half to three 
atmospheres so that this technique is limited to 
sixty feet. 

Taking the case of an “‘infinite’’ stay (over four 
to five hours) in 300 feet of sea water, the ascent 
curve for the seventy-five minute tissue as calcu- 
lated from Eq. (8) is shown in Fig. 3. At sixty 
feet we know the tissue is saturated at twice 
(n=2) absolute 60-ft pressure or 2X (60+33) 
= 186 ft. The diver may safely come out of the 
water when his tissues have a nitrogen saturation 
corresponding to twice atmospheric or 66 ft. 
Thus P’ = 188, P’’ = 66 and from Eq. (9), T’’, the 
time spent in the Davis chamber, =113 min. 
Comparing this modified ascent curve with the 
one calculated for air throughout (Fig. 3), it is 
apparent that this technique is very worth while 
for very long dives. 


HELIUM-OXYGEN MIXTURES 


In 1921 it was suggested that helium replace 
the nitrogen component in normal air for respira- 
tion at high pressures.‘ It was thought at this 
time that the lighter molecules would leave the 
tissues rapidly and cut down decompression time. 
Such is not necessarily the case* as a considera- 
tion of Eq. (8) will readily reveal. 

In Eq. (8) the only factors that depend on the 
kind of gas used are m and 7;. Helium, with a 
diffusion rate five-halves that of nitrogen and a 
water solubility of only three-fifths, has a maxi- 
mum 7, of twenty to thirty minutes.’ The 
normal duration of deep (200 ft or more) dives 
is generally not more than twenty to thirty 
minutes according to Navy regulations. For such 
short dives, the ‘‘slow”’ tissues do not play a 
major role in regulating the ascent so that 
helium does not give much assistance. Also, 
helium being more mobile appears to form 
bubbles more readily than nitrogen and requires 
a safe m value of 1.7 rather than 2.0. Taking as an 
example a twenty-minute stay at 300 ft, the 
tables‘ for helium give a total ascent time of one- 
hundred-and-thirteen minutes with pure oxygen 
breathing starting at the sixty-foot level. A 
normal air dive under identical conditions will 
take eighty-four minutes from the sea bottom to 

‘Navy Department, Bureau of Ships, Diving Manual 
(1944). 


5 J. F. Fulton, Decompression Sickness (W. B. Saunders 
Company, Philadelphia, 1951), Chap. III. 
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surface. Should the diver be trapped at a great 
depth for any length of time, however, his ascent 
would be much more rapid with helium than 
with air since his ‘‘slow” tissues would now 
determine the ascent curve. Helium is mainly 
advantageous in that it does not produce the 
narcotic effect under high pressures that is 
always associated with nitrogen-air. 


DISCUSSION 


Although this “‘classical’’ theory of decom- 
pression explains qualitatively a large body of 
phenomena associated with prevention of the 
bends, it is no longer trusted as a complete 
theory of nitrogen desaturation. There is doubt 
as to the physical reality of a variety of tissues 
having fixed half-saturation periods.® Actually it 
appears that when the body is suddenly moved 
from high to low pressure, nitrogen does not 
begin seeping from all tissues into the blood. 
Rather, the “slow” tissues rapidly accept nitro- 
gen from the quickly desaturating tissues and 
thereby act as ‘‘buffers’’ against sudden bubble 
formation in these “‘fast’’ tissues. Wartime tests 
on high altitude bends indicate that a maximum 
T,; of seventy-five minutes is dangerously con- 
servative and that, in fact, this quantitity can 
vary in the same individual depending on his 
work and surrounding conditions. 

While the value of around 2 was fairly satis- 
factory for the conditions that Haldane con- 
sidered (200-ft maximum depth), wartime re- 
search showed that m was variable for different 
tissues, being as high as 5 in ‘‘fast”’ tissues and 
as low as 1.7 in “‘slow”’ tissues. Also, appears 
to be a function of total pressure as well, being 
around 1.75 at 200 feet underwater and 2.7 at 
high altitudes for the same tissue. 

In spite of these added complexities, it must 
be stated that Haldane’s contribution to the 
field was immensely important and showed a 
sure grasp of physical thought that was lacking 
in previous workers. His tables have saved many 
lives and permitted much useful work to be 
accomplished underwater and at high altitudes. 
Whether a completely accurate model of the 
physical processes occurring during bubble forma- 
tion is possible remains to be seen, but the problem 
is an interesting one and its solution becomes 
more urgent as men push deeper under the sea. 
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When changing units and in dimensional analysis it is usual to employ the relation, physical 
quantity = measure X unit. Neither the “physical quantity” nor the “unit’’ is supposed to be a 
real number and yet the formalism is normally effective. This fact is given a rational explana- 
tion and the ideas involved are used to discuss certain questions connected with the rationaliza- 


tion of the electromagnetic equations. 


INTRODUCTION 


LTHOUGH the Giorgi system of electrical 
units (now sponsored by the International 
Electrotechnical Commission) is gradually com- 
ing into widespread use, its relationship with the 
older unrationalized systems is still confused (see, 
e.g., Landolt and de Boer).'! We believe that 
the trouble is caused by the use of the concept 
“‘physical quantity’’ which has objectionable 
metaphysical associations. Two variants of the 
rationalization process have been recognized, and 
it is necessary to refer to fundamental principles 
(on lines inspired. by O’Rahilly)? in order fully to 
appreciate their relative merits. 


PHYSICAL QUANTITIES, MEASURES, AND UNITS 


There is a widespread notion (for which Max- 
well was largely responsible) that ‘physical 
magnitudes” or “physical quantities’’ can be 
regarded as the product of ‘‘measures’’ (pure 
numbers) and “‘units’’ and that equations in 
physics are relations between “physical quan- 
tities.” See, for example, Guggenheim? or Page.* 
Thus it has become fairly common to write 


quantity = measure X unit, (1) 
{x} =xxt, 
where neither {x} nor xt are supposed to be 
ordinary numbers. The “‘units’’ are closely con- 
nected with the ‘“‘dimensions”’ of the correspond- 
ing quantity and are indeed separated out from 
equations of physical quantities and dealt with 


1M. Landolt and J. de Boer, Rev. Gén. de I’Electr. 60, 
499-508 (1951). 


2A. O’Rahilly, Electromagnetics (Longmans Green and 

Company, London, 1938), first edition, Chap. XIV, XV. 
3E. A. Guggenheim, Phil. Mag. 33, 479-496 (1942). 
*C. H. Page, Am. J. Phys. 20, 1-4 (1952). 


by methods similar to those of dimensional 
analysis. Such methods are normally foolproof 
and are, in fact, the quickest way of changing 
from one system of units to another. Unforiu- 
nately, the rational basis for these manipulations 
is rarely understood and this can lead to errors 
when discussing rationalization. 


FUNDAMENTAL PRINCIPLES 


The “measure” of a basic physical quantity 
(usually mass, length or time) is the ratio—found 
by a well-defined physical operation—between 
the quantity in question and an arbitrarily 
chosen conspecific unit. The measure of a derived 
quantity (perhaps already recognized in a quali- 
tative sense) is found by means of an algebraic 
defining equation relating measures of basic 
quantities. This equation may also contain nu- 
merical constants—such as s—chosen for con- 
venience. Measures are thus ordinary numbers, 
and here we shall denote them by plain letters. 

We shall use the symbol {x} for the verbal 
phrase ‘‘some particular example of the phe- 
nomenon whose measure is denoted by x,” and 
the symbol xt for the phrase ‘“‘any example of 
the phenomenon whose measure is 1 on the 
chosen system of measurement,” or, more 
shortly, “the unit of x.’’ Different systems of 
measurement will be distinguished by subscripts. 

For simplicity, the case of velocity, whose 
measure equation is v=kl/t, is considered. The 
“definitional constant’”’ k is usually taken as 1, 
but its explicit inclusion here enables the argu- 
ments to be generalized without any difficulty. 
The measure of {v} depends on the choice of k 
and of the units of length and time. For one 
particular choice we use the label 0 and regard 
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this system as a reference system. We use labels 
1, 2, --- for other choices and then for any of the 


systems labeled 0,1,2,---,7,---, ‘‘measure- 
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ratios” are defined as [x,]=xo/x,, where x can 
be 1, t, v, etc. In particular [xo ]=1. The follow- 
ing statements are then obvious: 


(A) The measure equation v,=k,J,/t, is valid (by definition) for each r. 


(B) The measure-ratio equation [v, ]=[k, ][/,]/[t,] is valid for each r. 
(C) [1-] t/l, measure of /,t in terms of Jot 


[l.] Io/le measure of /,t in terms of Jot 


l, measure of {/} in terms of /,t 


l, measure of (the same) {/} in terms of /,t ; 


[1,]/(l.] is independent of the choice made for system 0. 


[v,] measure of v,t in terms of vt 


(D) - 


[v.] measure of v,t in terms of vot 


(E) Vo = 0401 ]=V2[ V2] = tee =9,[ 0, ]= wioie 


The same results hold for any derived quantity. 

So far, we have only used the symbols {v} and 
vt as shorthand forms; we have not manipulated 
them in an algebraic manner. On searching 
among the purely algebraic numbers occurring 
in (A)-(E) above to find those whose manipula- 
tive properties are the same as those commonly 
ascribed to the ‘“‘unit,’’ the ‘‘measure’”’ and the 
“physical quantity,” the following conclusions 
can be drawn: 

(a) The numbers [1,], [/2], [#1], [te], etc., 
have properties formally identical with those of 
words like “‘foot’”’ and ‘‘second’’ when these are 
manipulated in the traditional way. For example, 
if systems 1 and 2 have the foot and the yard as 
their respective units of length, (C) abbve is 
311 ]= [le] (compare “3 feet = 1 yard’’). It should 
be noticed that although 3[/,]=[l2] is an ordi- 
nary arithmetical equation, the numerical values 
of [/,] and [72] depend on the choice of system 0; 
this fact is only important if two different refer- 
ence systems are used (see later). 

(b) The numbers /;, v:, etc., have their usual 
interpretation; that is, they are measures. 

(c) Numbers such as vp have the same manipu- 
lative properties with respect to [v,] and », as 
‘‘physical quantities” do with respect to “units” 
and ‘‘measures” on the orthodox treatment. For 
example, v)=2[01] corresponds to an equation 
like (1), {v}=v.0:t, which is not an ordinary 
algebraic equality. Since the 0-system is entirely 


arbitrary, equations equivalent to those involving 
‘physical quantities’ have not the special signifi- 
cance that the words ‘‘physical quantity” appear 
to imply. 

The above discussion is illustrated by the solu- 
tion of two problems: 

(1) What is the measure on system 2 (yard, 
minute, k=10) of a velocity whose measure is 5 
on system 1 (foot, second, k=1)? 

In this case, [l2]=3[/1], [te]=60[¢,], and 
[ki ]=10[k- ], so that 


Sf: J=S5Lki J /T1) 
= 1000[ 2 I[/21/[te] 
= 1000[ v2). 


The answer is 1000 and the novel features are the 
use of purely arithmetical symbols and the 
manipulation of k. 

(2) The power required to activate a solenoid 
of cross section A with current J through N turns 
of wire of density p, resistivity o, and total 
winding mass M is 


P=4rpcA(NI)?*M—. 


What is the power in watts when p is given in 
Ib-ft-*, o in ohm cm, A in in.?, J in amp, and M 
in lb? (Taken from Guggenheim.)* 

Consider the following four numbered sys- 
tems: (1) m, kg, sec, ohm; (2) in., kg, sec, ohm; 
(3) cm, kg, sec, ohm; (4) ft, lb, sec, ohm. The 
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problem is to find P; in terms of pu, 03, Ao, Ih, 
and M,. The arbitrary reference system is, as 
before, denoted by 0. Then, 


P,[P1] = Po=4mpoo A o( NI)? Mo 
=4mpaLps loslo3 JA oLA 2 | 
XN TY PMoLMay 
=41pi03A o(NIy)?M 
XLes Cos L421 PLM" 
=4mp.o3A (NY MolM, Lat 
XCRsJLLs Le PL PLM, 


where [R] and [L] are measure-ratios of re- 
sistance and length, respectively. But 


Ci PCRs J=(h PCR =(Pil, 


so that 


PLP: ]=42pyo3A2(NI;)*?M 
X(Ls)[L2PLLat [Pi] 

4x pso3A2(NI)? 

K = 


ineecienttel [Pi]. 
30.48 X 144 M, 


The multiplier of [P,] in the last formula gives 
the measure in watts. In the conventional and 
simpler solution, given by Guggenheim, {P} 
would be substituted for Py and “1 ohm” for 
[Rs], etc. 

These two examples are perhaps sufficient to 
show how the traditional method of dealing with 
changes of units can be justified. It is worth men- 
tioning in passing that, by the use of measures 
and measure-ratios, the whole of dimensional 
analysis can be stated in terms of ordinary num- 
bers. Such a treatment (as O’Rahilly? has shown) 
avoids many of the troublesome metaphysical 
questions associated with the orthodox proced- 
ure. Birkhoff> has also sketched a treatment 
of dimensional analysis in terms of ordinary 
numbers. 


RATIONALIZATION 


The form taken by the rationalized electro- 
magnetic equations is already agreed upon, but 
some physicists treat the symbols as measures 
and others as ‘“‘quantities.”” It is now necessary 
to examine the two variants of rationalization 
(discussed by Landolt and de Boer)! known as 


5G. Birkhoff, Hydrodynamics (Princeton University 
a Princeton, New Jersey, 1950), first edition, Chap. 
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“rationalization of units” and ‘‘rationalization of 
quantities,” to see the consequences of the two 
ways of interpreting the symbols. Only the case 
of the magnetic field {H} is considered because 
the treatment of the other quantities affected by 
rationalization is fully analogous. Following -the 
qualitative concept of field, the operation for its 
measurement is decided upon and a measure 
equation chosen so that numerical values can be 
assigned. One possible process defines the field 
in terms of the current {J} flowing in a ring of 
radius {R}, the measure of {H} at the center 
being given by the measure equation 


H=AI/R, (2) 


which contains a definitional constant A. In 
what follows, current and length are taken as 
basic quantities. The measure H then depends 
on the choices made for A, for the unit of current, 
and for the unit of length. 

Two choices of A are of interest: for the un- 
rationalized system (denoted by u) A=2r=Au, 
for the rationalized system (denoted by 1) 
A=}3=4A,; and also two different measuring sys- 
tems for {I} and {R}: for the e.m. system (de- 
noted by m), {I} measured in emu, {R} in cm, 
for the practical system (denoted by )), {J} 
measured in amp, {R} in meters. With these 
choices, there are four possible measures for the 
magnetic field: 


To convert from one measure of {H} to another 
we already have A,,/A,=42 and need in addition 
the measure-ratios I,/Imn=10 and R,»/R,=100. 
(The latter ratios are the inverses of the ratios of 
the measures of the units—the measures being 
expressed in terms of an arbitrary conspecific 
unit in accordance with (C).) Either directly 
from Eqs. (3) or else by using the measure-ratio 
equation [H]=[A ][J]/[R], it is found that 


1 _ 4a 
— Fy» = 42H m= H,». (4) 
00 100 
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FIRST VARIANT—RATIONALIZATION OF UNITS 


A reference system 0 is chosen in which the 
measure equation 


Hy=Aolo/Ro, (S) 
corresponds to an “equation of quantities’’: 
Ay=2x, A,=}, 
so that 


[Au]=Ao/2e and [A,]=2Ao; 


also 


(MuJ=(A.JUI/CR), 4 J=(4,JU7/C)- 


Hence 


1 
[He J=7 1H, ©) 


where the labels m or p are omitted since either 
could be used. 
Also 


[Hm ]=[A JU mJ/LRm J 
and 

[H,]=(AILI>]/(Re], 
but 


CImVCIp]=10 and [R,]/[Rm]=100. 


Hence 
LH» ]=1000(H,], (8) 


where the labels u and 7 are omitted since either 


could be used. From Eqs. (7) and (8) it follows 
that 


1 1000 
[Hum ]=1000(H.»]=—LHim]=——[Hp], (9) 
An Ar 


which correspond to ‘‘unit”’ equations consistent 
with the measure equations (4). Conventional 
methods can be used to recover Eqs. (9) on con- 
dition that the following device is adopted. 
Introduce the words “‘rat’”’ and “‘unrat’’ related 
to [A,] and [A.,], respectively, in the same way 
as the word “yard” is related to the number 
[/2] in example 1 above. Just as ‘‘3 feet =1 yard” 
so “1 rat=4m unrat.’”’ Then Eqs. (9) corre- 
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spond to 


emu 
1 oersted* =1 unrat —— 
cm 


amp amp * 
= 1000 unrat —— 
m 


=1 unrat 


100 


1 emu 1000 amp * 


=— rat -—— = io , 
4 cm an m 


(10). 


where the terms marked * correspond, in order, 
to those in Eqs. (9). 
Noting that 


Ao [J] 


(IJ 
LH. j= [R] 


[R] 


where the value of Ao is arbitrary, other per- 
missible sets of names for the four units of {H} 
(some of which have been used in the literature) 
can be given. Thus 


1 
— emu/cm, 
an 


and [H,]=2Ao 


Himt: 1 emu/cm 


[Pn 


1 
|r - amp/m 


T 


H,,': 1 amp/m 


is one possibility, 
and 


magts 
be; 


It is not permissible to take the names from 
different sets (e.g., H,pt: 1 amp/m and Hu»pt: 
1 emu/cm) and then to expect the arithmetical 
properties of the names to be preserved. This was 
indeed the difficulty encountered by Landolt and 
de Boer.! (We suggest that the best way of naming 
the units is to adopt the “rat’’ and “‘unrat’’ as 
indicated above, since this preserves the arbi- 
trary nature of Ao.) 

As there is only one reference system 0, there 
is only one ‘‘equation of quantities’ correspond- 
ing to Eq. (5). (The rationalization process has 
been one of changing units and measures in in- 


1emu/cm, HAynt: 


An ae 
Buy: 


1 amp/m, 4x amp/m 


is another. 





A CALCULUS OF MEASURES 


verse ratios.) By suitable choice of the arbitrary 
Ao, the ‘equation of quantities’’ can be given 
the same form as the rationalized or unrational- 
ized equation of measures but not both simul- 
taneously. The process of rationalization of units 
is entirely logical and satisfactory to those who 
use measures exclusively, but it is extremely in- 
convenient for those who wish to interpret both 


H=2rlI/R and H=}I/R 


as equations of quantities. The desire to use the 
idea of ‘‘quantity’”” has apparently led to the 
invention of the second variant of the rational- 
ization process, which will now be discussed. 


SECOND VARIANT—RATIONALIZATION 
OF QUANTITIES 


In this variant two measures of {H} are de- 
fined by equations in which the definitional con- 
stants are chosen and kept fixed. What we call 
the ‘‘unrationalized’”’ and 
ures” of {H} are defined by 


“rationalized meas- 


H=2rI/R_ (unrationalized) (11) 


and 
H'=3I1'/R' (rationalized). 


(12) 


To each of these there corresponds a separate 
reference system in which the measures are 


He= 2rIo/Ro 
H,! = 479 /Ro 


(system 0, unrationalized), 


(13) 
(system 0’, rationalized). (14) 


These analogs of equations of quantities have 
the same form as measure equations (11) and 
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(12), and the quantities themselves have been 
called the ‘unrationalized” and “rationalized 
magnetic field.” 

Changes of units (and measures) between 
either unrationalized systems or rationalized ones 
are carried out in the normal way, e.g., [Hn] 
= (In ]/[Rn]=1000(T,]/[R,]=1000LH,] (un- 


rationalized case), and similarly 
(H.,.’ / =1000CH,’] (rationalized case). 


Because [Hm ]=(Im]/[Rm] and (Hn!) =(Im’]/ 
[Rm’], these measure-ratios for {H} can be given 
the same name, emu/cm. Similarly the two prac- 
tical units can both be called the amp/m. 
Clearly the unit calculus cannot be used to trans- 
form unrationalized measures to rationalized 
ones or vice versa. 

In order to transform between the four chief 
measures of {H} (to recover Eqs. (4)) the follow- 
ing must be noted: 

(a) The unrationalized e.m. and _ practical 
measures of {H} are related inversely as their 
units, and similarly for the rationalized measures. 

(b) Unrationalized and rationalized measures 
of {H}, if {7} and {R} are measured in the same 
units in both cases, are related by 


H=4rH'. (15a) 


(c) Users of quantities (instead of rationaliz- 
ing measures as in (b)) will rationalize quantities 
using the equation of quantities analogous to 


Hy =4nrHy’, (15b) 


but the result is the same as in (b). 


RATIONALIZATION—SUMMARY 


Rationalization of units 


The magnetic field has four 

measures and four unit names: 
Hum unrat emu/cm (oersted) 
Hy» unrat amp/m 


Hym rat emu/cm 
H,, rat amp/m 
Also, 1 rat =4m unrat 


Rationalization of quantities 


The field has two unrationalized 
measures and two unit names: 

H,, emu/cm (oersted) 

H, amp/m 

The field has two rationalized 
measures and two unit names: 

H,,’ emu/cm 

H,' amp/m 

Also, Hy, =4H,,' and H,=41rH,' 


Common to both variants are 
1emu=10amp, 1m=100cm. 
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Therefore 


1 
Hum( = Hn) =——. 


1 

H, (-_#,) 

1000 \ 1000 ~ 
; 4 

= 4rH,m(=4rH,') =—— 


Ps ( ~~ 11’) 
*\ 1000 ” 


1000 


in agreement with Eqs. (4). 

As an example, the ways of specifying the field 
at the center of a single circular turn of radius 
1 cm carrying a current of 1 amp are given: 


(a) rationalization of units, 


H=6.28X10- unrat emu/cm 
= 6.28 X10? unrat amp/m 
=5X10- rat emu/cm=5 X10 rat amp/m. 


(b) rationalization of quantities, 


magnetic field (unrationalized) 
= 6.28 X10-! emu/cm = 6.28 X10? amp/m, 


magnetic field (rationalized) 
=5 X10 emu/cm=5 X10 amp/m. 


The same numbers (measures) appear whichever 
variant is used. 


R. BARKER AND R. O. DAVIES 


CONCLUSIONS 


A calculus of measures has been presented 
which is complete in itself. It has been applied 
to explain the successes of (a) the interpretation 
of the symbols in equations as ‘‘physical quan- 
tities” or ‘‘magnitudes”’ and (b) the manipulation 
of unit names when changing units. It is quite 
unnecessary to use ‘‘quantities’’ and, in view of 
the metaphysical associations of the word, 
dangerous when investigating unsettled ground. 
The manipulation of unit names does, however, 
provide the easiest way of changing units and 
can be regarded as a mnemonic for the more 
logical manipulation of ordinary numbers (i.e., 
the measure-ratios) in accordance with the rules 
of ordinary algebra. 

We have presented two variants of the ra- 
tionalization process. Those who insist on using 
“‘equations of quantities’’ must employ the sec- 
ond variant because only this one provides two 
equations of quantities (13) and (14) of exactly 
the same form as the measure equations (11) and 
(12), respectively. Those who use measures ex- 
clusively can use either variant with success. In 
practice the choice is immaterial. It must, how- 
ever, be clearly stated either in the name of the 
unit (first variant) or in the name of the measure 
or quantity (second variant) whether or not 
rationalization has been employed. 


To Secretaries of Regional Sections of AAPT 


The editor is very anxious to publish reports of Section 
meetings as promptly as possible after they are held. There 
is an unavoidable delay of about three months from the 
time a manuscript leaves my office till the time of its 
appearance in print. These three months are consumed by 
work in the Publications Office of the Institute, in the 
printing house, and by the reading of proof, and there is 
no simple way of shortening the period. 

You can, however, help to assure the most rapid pub- 
lication of your Section news by sending me a report of 
each meeting immediately after the close of the meeting. 
To save time, the report should be written following 
meticulously the standard pattern that we use, as for 
example, on pages 73, 74 of the January, 1953 issue or 
pages 461-463 of the October, 1952 issue. The material 


must be typewritten and everything must be double- 
spaced to allow for editorial markings. Please note that 
this means everything. The contributed papers need not be 
numbered serially if there are only a few. Numbering is 
helpful if there are more than six or eight. 

Abstracts of papers should be included whenever they 
are available and I would urge secretaries of regional 
Sections to collect such abstracts either before or at the 
actual meeting. Abstracts which deal with abstruse re- 
search problems may well be omitted in order that those 
which deal with matters of teaching may be made a little 
longer. Any work that you do in shaping, clarifying, and 
improving abstracts will be much appreciated. 

Tuomas H. OsGoop 
Editor 





Longitudinal Vibrations of a Vertical Column by the Method of Laplace Transform 


L. H.. Har 
Santa Barbara College, University of California, Santa Barbara, California 
(Received October 15, 1952) 


An interesting illustrative problem in mechanics is the motion executed by a vertical rod or 
column when it suffers sudden loss of support at the bottom. A solution is presented by the 
method of Laplace transform in two alternative forms: one represents the motion as a super- 
position of free fall and two oppositely directed progressive waves equal in amplitude and form; 
the other, as a free fall and a pulse, initiated at the bottom and propagated up and down the 


column upon reflection at the ends. 


HYSICALLY a simple type of problem, 

whose solution may deserve illustration in 
a course in mechanics or applied mathematics, is 
the motion executed by a vertical, elastic column 
when it suffers sudden loss of support. Since the 
Laplace transformation applied to this case leads 
to a solution intelligible and natural from a 
physical point of view, the problem serves well 
also to illustrate this important method which 
today increasingly appears in the literature of 
physics and engineering. (In a recent account of 
transform applications,' for instance, are pre- 
sented examples in atomic and nuclear physics, 
hydrodynamics, elasticity, vibrations, and heat 
conduction.) 


DERIVATION OF THE EQUATION OF MOTION 


Consider a vertical, elastic rod or column of 
length L cross section A. Let it be subject to 
atmospheric pressure Po; initially supported at 
the bottom, and in static equilibrium under 
gravity forces. Suppose the support to be re- 
moved suddenly. This action will be represented 
as an instantaneous drop in pressure at the bot- 
tom of the column from the static value, Po+pgx, 
to the atmospheric value Po. We take the case 
of negligible internal and external friction. Let 
the top of the column be initially at the origin; 
x =the distance from the origin to any point or 
particle on the column under consideration; 
u=particle displacement; v=total particle ve- 
locity; P =total pressure; P,=P —Po, the excess 
pressure; p=mass density; w=weight density; 
and E=Young’s modulus. With reference to 


Fig. 1, a small slice of medium lying between x 


17, N. Sneddon, Fourier Transforms (McGraw-Hill Book 
Company, Inc., New York, 1951), first edition. 


and x+-Ax experiences a net force, (—dP/dx)AxA 
+wAxA. By Newton’s law we have 


au oP 

(AxA p)— = ——AxA+wAxA, 
dt? ox 
d*u oP 


p—= ——+w. 
di? Ox 


Since the pressure at any point is P»)+P., 
dP /dx=dP,./dx. The fractional change in volume 
of the slice, AV/V, is simply du/dx. Taking the 
usual linear relation between the excess pressure 
P. and the relative compression, we have 
P.= —EAV/V = —Edu/dx, and therefore 0P/dx 
= —Eé*u/dx*. Letting a? = E/pand limiting study 
to cases where the motion is small enough so 
that total time derivatives can be approximated 
by partial time derivatives, which is true in the 


usual applications, we arrive at the equation of 
motion, 


au 1 du g 
—-——=--. (2) 
os a oF a? 


The initial and boundary conditions are as fol- 
lows: Initially the column is at rest, so the total 
velocity v= du/dt=0. The pressure has the simple 


static distribution, Po-+wx or P.=wx; associated 
with this excess pressure there is an elastic dis- 


Fic. 1. Volume 
element of medium 
under compression. 
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placement of the fluid particles u. It is conveni- 
ent to let «=0 when x =0; then, for compression, 
u will be negative when x is positive. The initial 
displacement is —x?g/2a?, which can be verified 
by noting that it satisfies the stress-strain rela- 
tion, P.= —Edu/dx, and gives the required 
value zero, when x is zero. When the pressure at 
the bottom is instantaneously dropped to Po, we 
have P,=0, du/dx =0, when x = Land when x =0. 
We have therefore the system of equations 

Ou 1 du 


Ox? 


mi 


a? of ~ =a? : 


—x*g 


2a? 


when t=0; 


P.=0; —=0, 


at x equal to 0 or to L 
ox 


when ¢>0. 
THE METHOD OF LAPLACE TRANSFORMATION 


It may not be amiss to review the idea and 
some of the definitions of transform methods. 
Sneddon! gives a nice orientation in his opening 
chapter, and we follow his train of thought. 
Given a set of equations which define implicitly 
a desired function, it is a common technique to 
introduce a transformation of some kind to 
facilitate the solution. An integral transforma- 
tion defined as follows is one of the useful types: 
Let f(t) be a function of the real variable /, 
which ranges over the interval (0, ©), and let 
K(p,t) be a selected function of ¢ and the 
variable », which may be complex. When the 
integral 


f f()K(p, tat 


is convergent, it defines a function f(p) of the 
variable p, which is termed the integral trans- 
form of the function f(#) by the kernal K(f, 2). 
Simple and serviceable kernals are e~?', ¢?-', 
cos(pt), and sin(pt). The first of these is the 
Laplace kernal giving the Laplace transforma- 
tion; the second, the Mellin kernal; the others, 
Fourier cos and sin kernals. In a problem we seek 
to transform the determining equation for the 
desired f(t) into an equation which the transform 
function, f(), must satisfy. It is readily verified 
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that a differentiation of the original function, 
f(t), corresponds merely to a multiplication of the 
Laplace transform by its independent variable 
p with the addition of a constant, which for a 
very broad class of functions is —f(0). (Further, 
in many applications we can adjust conditions 
so that f(0)=0.) Upon solving the transform 
equation, we need of course a way back from 
f(p) to f(t). The more elementary and frequent 
transform pairs are cataloged in transform 
tables; for other cases we look to inversion 
theorems such as the Fourier-Mellin inversion 
theorem, which will be employed later. For linear 
differential equations with constant coefficients, 
a simple and customary way to effect the transi- 
tion from the equation for the original function 
to that for the transform or image function, is to 
multiply through the original equation by e-?! 
and to integrate with respect to ¢ from 0 to o. 
In the Eqs. (3) u is a function of two inde- 
pendent variables rather than one; we shall 
transform with respect to only one, ¢. Multiplying 
Eqs. (3) by e~?‘ and then integrating from 0 to «: 


Ps au 1 © Ou -—g f” , 
f e-?\—_dt —— f e~?'—dt = —— f e-P dt. 
0 Ox? a? 0 or? a? 0 


For the first term, 
- 07u 0? - Ou 
f e~?'—_dt = —— f e—?'udt = —— 
0 Ox? Ox? J Ox? 
(by definition). 


For the second term, integrating by parts, 


; Ou -[ “| 
— e-?—_ dj =—j ¢-?*— 
a? 0 0? a? ot 0 


2 oo 
4 epee f e-?'udt. 
a a Jo 


The first number is zero, for du/dt=0, when t=0 
and (e~?‘du/dt)—-0 as t-, because the factor 
e-?* becomes small more rapidly then du/dt 
becomes large with increasing ¢. The second 
member yields pgx?/2a*. The third member is 
pi /a*. Finally, 

=—— 2 


a? 0 


—gi 
edt =. —, 


a p 
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Hence the equation to be satisfied by the Laplace 
transform (the so-called subsidiary equation) is 
an ordinary differential equation in x, 


ou p pg g 
——a=—- (4) 
0x? a? 2a*  a*p 


with boundary conditions d%/dx=0, when x=0 
or L. These boundary conditions follow from 
du/dx=0 at x equal to 0 or L, since 


:. . - Ou 
=— e~?'udt = f e~?'—dt=0. 
Ox Ox 0 Ox 


The homogeneous part of Eq. (4) has the familiar 
solution 


px 
tinom = Ae?*/*+ Be-rz/e = A cosh (= -0) ‘ 
a 


By inspection, a particular solution is @part. 
= —gx*/2a’p. Thus the solution is 


u=A cosh(“=—0) =. (5) 
a 2a*p 


Introducing the boundary conditions, d%/dx =0 
when x=0 or L: 


ou p 
—=-A sinh (= - “) = ; 
Ox a a*p 


When x=0, 


p 
=—A sinh(—#)=0; therefore 
a 


Ox 


When x=L, 


ou p Lg 
—=-A sinh(pL/a) — 
ox a a*p 


therefore 
Lg 
~ ap*sinh(pL/a) | 


Thus the solution of the subsidiary Eq. (4) with 


boundary conditions is 


Lg cosh(px/a) xg 
ae ere (6) 
ap? sinh(pL/a) ~ 2a%p 


Then u, the inverse Laplace transform of @, is 
equal to the sum of the inverse transforms of the 
two right-hand terms of Eq. (6). It is readily 
verified that the inverse transform of 1/p is 1, 
and thus that the inverse of the second term is 
—x’g/2a?. For the inverse transform of the more 
complicated first term, unlikely to be found in a 
table, we turn to the Fourier-Mellin inversion 
theorem,?~ which states that 


i ytico 
u=— f— evtateyap, 


211 y—iw 


where y is any constant greater than the real 
parts of all the singularities of a(p); i.e., integra- 
tion is along any line parallel to the imaginary 
axis and to the right of the singularities. Thus we 
have to evaluate the complex integral 


Lg ytio 
eP 


2ria 


cosh (px/a) dp 
sere Pp 


This will be done in two ways: (1) 
integration; (2) by residue integration. 


(7) 


y—io 


by series 


EXPONENTIAL SERIES EXPANSION OR 
HEAVISIDE WAVE FORM 


The physical action at the removal of support 
(an abrupt drop in pressure) suggests the idea 
of a pressure pulse (subpressure) which is initi- 
ated at the end where x=L, travels toward the 
end where x=0, is reflected, travels back and is 
reflected again at the bottom, etc. If the integral 
(7) is evaluated by expressing the integrand in a 
series of exponentials (the so-called Heaviside 
wave method) we are led toa solution whose 
form expresses these physical features. 

The integrand of the inversion integral (7) 


2R. V. Churchill, Modern Operational Mathematics in 
Engineering (McGraw-Hill Book Company, Inc., New 
York, 1944), first edition. 
3H. S. Carslaw and J. C. Jaeger, Operational Methods in 
Applied Mathematics (Oxford University Press, London, 
1941), first edition. 
*M. F. Gardner and J. L. Barnes, Transients in Linear 


Systems (John Wiley and Sons, Inc., New York, 1942), first 
edtion, Vol. 1. 
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may be expanded as follows: 
cosh(px/a) 
ert = 


(enz/e + e~Pz/a) 
ecpennensseneesnicneeensnnnecee, pt 
sinh(pL/a) 


(erLla —_ e~PLie) 
e~PLia (erz/e+ e~Pzla) 
(1 ini e~*pL/e) 


= eP!—pL/a(epzia_+ e~P2/a) 


=ePt 


x (1 + e-2pLla+ e~tpLiat .. -) 
« 
as > eP(at—(m+1)L+z)/at eP(at—(2m+1)L—z)/a_ 


m=0 
The displacement may, therefore, be written 


Lg ytim 
=—_——_ . {eP(ot—Qmt1)L+2)/a 
dp xg 


+ ep(at—@mt1) L—z)/a\___ s 
Pp? 2a” 


The excess pressure 


ou wa? du 
[--# 
Ox g Ox 


wh pit™ o@ 
P,.=wx —— DY {ertet-@mt1L+2)/a 
mt —ia m=0 


is then 


— ep (at—(m+ DL aie) - 


~ 


(8) 


The integral consists of a summation of terms of 


the form 
1 y+io 


dp 


Pr, 


y—tw p 


which has the value 1, if k>0 and the value 0, 
if k<0O (refeference 3, p. 165), i.e., the unit step 
function, U(k). Replacing k, we have U(t 
—[(2m+1)L+x]/a), which equals 0 when 
t<[(2m+1)L+x]/a and equals 1 when 
t>((2m+1)L+<x]/a. The excess pressure dis- 
tribution accordingly may be written 


0 2 1)L—- 
P,=xw—wL Z U ae) 


m=0 a 


~_ u(: oo) . (9) 


The total pressure is simply Po+P.. 


Qari 
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The fluid velocity distribution, vf =v, = du/dt], 
is 


vet Eu u(:- (2m+1)L — ==) 


4 “ aa). -" 


DISCUSSION OF HEAVISIDE WAVE SOLUTION 


The displacement distribution will be discussed 
in the next section. 


Pressure Distribution 


To interpret expression (9), consider some 
point along the column, say the mid-point, 
x =L/2, and let the time taken for a disturbance 
to traverse the length of the column be r=L/a; 
then the expression (9) becomes 


wh 00 
oe —wL > U(t—(2mr+7/2)) 


m=0 


— U(t—(2mr+37/2)). (11) 


Consider now a time lapse less than that required 
for an elastic disturbance to travel from the end 
x=L to the mid-point of the column. This is 
a time for which t<7/2, and hence for which 
t<(2mr++17/2) for all values of m. Thus Eq. (11) 
reduces to P,=wL/2, the correct initial pressure 
at the mid-point. When t>7/2 but t<37/2, 
P.=wL/2—wL(1i —0) = —wL/2. When t>3r/2 
butt<57r/2, P.=wL/2 —wL(1 —1) =wL/2. When 
t>57r/2 but t<77/2, P.=wL/2—wL(1+1-1) 
= -—wL/2. The interpretation is clearly that 
a pulse of subpressure wl is initiated at the 
bottom at the withdrawal of support, :=0; and 
that when a lapse of time sufficient for the pulse 
to reach point x=L/2 has occurred (t>7/2), 
this subpressure arrives and is superimposed on 
the original pressure of wLl/2, leaving a net 
pressure of —wL/2. When a time lapse such that 
t>3r/2 but <57/2 has occurred, the pulse has 
been reflected from end x=0; has returned as a 
superpressure, restoring the initial value of the 
pressure at x=L/2. When t>5r/2 but <77/2, 
the pulse has passed on to the end x=Z; has 
been reflected as a subpressure again; and has 
again returned to the mid-point causing the 
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pressure to drop to —wL/2. In this manner the 
pressure pulse travels up and down the column, 
causing on upward passage a pressure decrement 
—wL, which added to the existing pressure 
wL/2, gives a resultant pressure —wL/2; and 
causing on its downward journey a pressure in- 
crement wL, restoring the original positive pres- 
sure wL/2. The pulse of superpressure is reflected 
at the ends as a pulse of subpressure and vice 
versa. Near the bottom of the column, x= Z—AL, 
the excess pressure fluctuates discontinuously 
between 0 and wL—AP. Pressure history dia- 
grams for three points along the rod are shown in 
Fig. 2. 


Velocity Distribution 


When x=L/2, Eq. (10) yields the information: 
the velocity jumps from 0 to gr for t>7/2 and 
remains at this value until ¢>37r/2 when it 
jumps to 2gr. It continues to increase in incre- 
ments of gr every time interval 7, which obvi- 
ously corresponds to the arrival of a velocity 
pulse upon reflection at the ends. When x=L the 
velocity jumps from 0 to gr for ¢>0 and there- 
after increases in increments of gr every time 
interval 2r. 


INVERSION BY RESIDUE INTEGRATION 


The singularities of the integrand in Eq. (7) 
are at the zeros of sinh(pL/a) which give simple 
poles at (pL/a) = +71, n integral, and at p=0. 
The residues for the former singularities are 
given by the formula® 


n+ 5 g 


e?* cosh(px/a) 


n=—0© 2710 


d 
—p* sinh(Lp/a) 
dp 


(Lp/a=+nr i) 


With the help of the relations between hyper- 
bolic, circular, and exponential functions this 
reduces to 

1 2L%g » (-—1)! 


271 wa? n=1 n? 


Xcos(nmrat/L) cos(nmx/L). (12) 


5 R. V. Churchill, Introduction to Complex Variables and 
Applications (McGraw-Hill Book Company, Inc., New 
York, 1948), first edition, p. 122. 


34 7 oe tpeses 


t in units of T/ 


Leterme 
t in units of Th 
(c) 
Fic. 2. Pressure history curves. 


For the residue at p=0, we form the Laurent ex- 
pansion of the integrand about »=0 by expand- 
ing the factors individually and combining, 
Lg ecosh(Lx/a) Lg (1+p+p'?/2+---) 
p? 
X (1+ p*x?/2a?+ ---)(@/Lp —Lp/6a+---), 
which gives for the coefficient of the 1/p term 
1 L’g xg Ff 
5 5 a 
2mit 6a? 2a? 2 


The integral (7) equals the sum of the residues of 
the integrand, expressions (12) and (13), multi- 
plied by 277. We get finally 


L*g 


6a? 


2L*g » (—1)"4 


sae x 


Xcos(nrat/L) cos(nrx/L)+ $gt?. (14) 


Applying a trigonometric identity, the solution 
for the displacement becomes 


L’g 


“= —— 


L’g Ps (-—1)*4! 


wat n=l on? 


nT nt 
x aad ae —x) toes | +4. (15) 





Fic. 3. Pressure wave form when ¢=0. 


The excess pressure [= —(wa?/g)(du/dx) ] is 


. (16) 


TT nT 
sin—(at —x) —sin—(at+~x) 
L i 


The fluid velocity is 


™ nT 
- et ere +gt. (17) 


DISCUSSION OF SOLUTIONS 
Displacement Distribution 
When ¢=0, Eq. (14) or (15) reduces to 


L’g 2L%g «© (—1)"4 
u= —— > —— cos(nrx/L), 


ra? n=1 n? 


which is just the Fourier expansion of the func- 
tion u= —gx?/2a’, the initial displacement dis- 
tribution along the column.® For each n, the 
trigonometric summation term in Eq. (14) gives 
a harmonic standing wave motion; and for all 1, 
a set of superimposed standing waves, constitut- 
ing, a complex standing wave or vibration of the 
column. The last term is the familiar free fall 
motion. Thus the column vibrates as it falls. 


Pressure Distribution 


When ¢=0, Eq. (16) is the Fourier series repre- 
sentation of wx, the initial pressure distribu- 
tion along the column. Terms of the form, 
sinnx(at —x)/L, represent progressive harmonic 
waves traveling in the positive direction, and the 


®W. E. Byerly, Fourier Series and a Harmonics 
(Ginn and Company, New York, 1893), p 
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sum, >> sinnr(at—x)/L, a series of harmonic 
waves—the harmonic components of a wave of 
more complex form. When ¢=0, this wave is 
pictured in Fig. 3. As time passes the wave form 
moves to the right with velocity a. Similarly 
> sinnz(at+x) represents a wave of the same 
form but moving in the negative direction, left. 
The two waves coincide in initial position when 
t=0 and together give }P.+3P.=P,.. They are 
illustrated in Fig. 3. The ensuing distribution of 
pressure for any time / can be obtained graph- 
ically by sliding one of the superimposed wave 
forms to the right a distance x equal to at, and 
the other an equal distance to the left and adding 
ordinates in the new positions to give the new 
pressure distribution. This recipe gives pressure 
history curves identical, of course, with those 
derived from the Heaviside solution. 

While the problem at hand is as readily solved 
by the method of Fourier series as by the Laplace 
transformation, there are in many cases ad- 
vantages to the latter. In a general way this may 
be seen in the fact that differentiation of a func- 
tion has merely the effect of multiplying the 
associated Laplace transform by its independent 
variable and adding to it the constant, —/f(0). 
Thus, in the case of linear differential equations 
with constant coefficients, the transformation 
converts ordinary differential equations into alge- 
braic equations and partial differential equations 
in two independent variables into ordinary differ- 
ential equations. In many instances simplifica- 
tion results also for linear differential equations 
with variable coefficients and for integral equa- 
tions. Nonlinear differential equations, however, 
are less likely to be good subjects. Since the 
integral of a discontinuous function may be a 
continuous or analytic one, the Laplace trans- 
formation also enables the study of discontinuous 
functions (such as occur in impulse or transient 
phenomena) through the medium of the more 
sedate transform functions. To be sure there does 
arise the added problem of evaluating the inver- 
sion integral when the transform is not one 
already cataloged. (It was this inversion task 
which may have made the solution of the present 
problem seem cumbersome.) However, the ex- 
change of a direct integration for an appreciable 
simplification of a differential equation is fre- 
quently a good mathematical bargain. 
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Since the symbols—names, or their abbreviations—for physical units may be incorporated in 
physical equations in exactly the same manner as are the letter symbols for physical quantities 
and the numerical symbols, it follows that the rules for selecting and using these unit-symbols 
should be formulated in the light of established mathematical procedures. A set of such rules 
is proposed. Symbols are suggested for 185 commonly used units and combining forms. 


ARIOUS standardizing and _ publication 
agencies in the physical sciences and en- 
gineering have recommended or adopted lists of 
“abbreviations” for the names of physical units. 
None of the lists, even the most recent ones, seem 
to have been compiled with much regard for the 
view that these ‘‘abbreviations’’ may be usefully 
employed as physical symbols, in the sense that 
they may be incorporated in physical equations 
in exactly the same manner as are the letter sym- 
bols for physical quantities and the numerical 
symbols; moreover, that these symbols (abbrevi- 
ations) are then subject to the same conventions 
and rules of operation as are the other mathe- 
matical symbols of physics. For this view to hold, 
a letter symbol for a physical quantity must be 
interpreted as representing both the magnitude 
and unit of the quantity. It is true that the letter 
symbols may be regarded alternatively as repre- 
senting magnitudes alone. Either one of the 
interpretations may be adopted, so long as one 
uses it exclusively throughout a particular phys- 
ical description. 

No single list of unit-symbols will be generally 
useful, however, unless it is devised mainly in the 
light of the first interpretation and its implica- 
tions; for this interpretation, unlike the second 
one, suggests that considerable economy of 
thought can be effected if the rules for selecting 
and using unit-symbols are formulated, not 
mainly arbitrarily, but in the light of the already 
established mathematical procedures of physics. 
As an analogy, consider an earlier time when 
negative exponents had not yet been defined and 
there were no rules for operating on quantities 
such as x~*. If standardization had been the only 
issue, any one of several definitions agreed upon 
by a majority might conceivably have done. 


However, a definition and rules of operation for 
positive exponents were already available, and 
so the sensible step was to adopt x-"=1/x" as 
the definition, thus making factors involving 
negative exponents subject to the same rules of 
operation as those with positive exponents. Inte- 
grative steps of this kind, that effect economy of 
thought and effort, and simplify learning, are of 
course a main objective of any science. 

In line with these ideas, the writer prepared a 
35-page mimeographed document on unit-sym- 
bols and sent it to a number of physicists in 
various types of educational institutions and 
industrial and Government laboratories. In- 
cluded in the document was a proposed set of 
rules applying to such symbols and a list of ap- 
proximately 500 primary and secondary units, 
together with the various ‘‘abbreviations’’ for 
each of them that are in more or less common use. 
In the light of the various comments and criti- 
cisms received, the present report has been pre- 
pared. Many of the units listed in the original 
document have been eliminated, the present 
Table I being confined to those basic units and 
combining forms that are believed to be in most 
common use. 


PROPOSED RULES FOR FORMING UNIT-SYMBOLS 


A number of the symbols listed in Table I do 
not conform with the following rules, but are so 
firmly established that attempts to change them 
would be futile. Exceptions to rules of language 
are a common occurrence and provide no good 
argument against having rules that facilitate a 
choice between several well-established symbols 
or the selection of a symbol for a new unit. 


(1) The symbol preferably should consist of two 
or three letters, and never more than four. Thus 
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Unit or combining 
form Symbol 
ab- ab 
acre acre 
ampere amp 
a 


angstrom A 

are [=10? m?] are 
atom atom 
atomic mass unit amu 
atomic weight unit awu 
atmosphere, standard atm, As 
atmosphere at 45° atmas, Aas 
bar [ =108 dy/cm?] bar 
barn [=10-% cm?] barn 
barrel bbl 
barye [=dy/cm?] ere 

bel bel 


b 
billion [=109] B 
Bohr magneton, electronic yo 
Bohr magneton, nuclear BI 
Brinell hardness number Bhn 
British thermal unit 
bushel 
calory 
Calory 
candle 
candlepower 
cent (acous.) 
cent (monetary) 


centi- [=10-2] 

centigrade degree 
(temp. difference) 

centigrade temperature 

centimeter-gram-second 
unit 

centimeter-of-mercury 

centimeter-of-oil 

cgs electromagnetic unit 

cgs electrostatic unit 

circular (adj.) 

coulomb 

cubic 

curie 

cycle 


day 

deca- [ =10] 
deci- [=10-!] 
degree (of arc) 


degree absolute 


degree Baumé 
degree centigrade 
degree Fahrenheit 
degree Kelvin 
degree Rankine 
deka- [ =10] 
diopter 

division 

dollar 


dozen 

dyne 
electromagnetic unit 
electron (charge) 
electron volt 
electrostatic unit 


erg 

Fahrenheit degree 
(temp. difference) 

farad 


foot 

foot-candle [ =lu/ft?] 
foot-lambert [ =ca/zft?] 
foot-pound-second unit 
fresnel [ =10!2 cy/sec] 
gal [ =cm/sec*] 

gallon 

gamma [10-5 oer] 


gill 
grad [ #107? rt. angle] 
grain 
gram (mass) 


gram of force 
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TABLE I. Proposed symbols for physical units. 


Examples of use; 
comments 


abamp, abcoul/cm? 
acre ft 

amp turn/m 

With prefixes: ma, ya 


Use sparingly 
atom/gm-awu, atom/mole 
1 amu =931 Mev 

1 awu =1.0002 amu 


Use wbar or dy/cm? 
| prefixes: db 


1 wo =9.27 X1072! erg/gauss 
1 wt =p0/ 1836 erg/gauss 


Btu/lb F° 
bu/acre 


See kilocalory 
ca hr, ca/m? 
See candle 


ct/gal, ct/kw hr 
= crowded tables, etc. 


ceifamn ane oo, 

100C° = =180F° 
See degree centigrade 
Use sparingly 


See also ab- 

See also stat- 

cir-mil [=0.7854 mil?] 

coul/m? 

cm? (never “‘cc’’), ft?/sec 

me, pc 

cy/sec [=hz] 

Often used with prefixes: 
ke, Mc 


Deprecated 


. dm 
deg/sec (ang. velocity) 
90 


See: degree Kelvin; d. 
Rankine 


0°C =32°F 

32. get (ice point) 
273.16 

491. ook 
Deprecated 


div/sec, div/pv 

dol/hr, dol/ton 

In tables, etc. 

doz/hr 

dy/cm?, dy cm (torque) 
Ambiguous; deprecated 


1 ev =1.601 X107"? erg 
Ambiguous; deprecated 
erg sec, erg/C°, erg/deg 
Btu/F° 


fd/m 

With prefixes: abf; uf 
ft/sec? 

lu/ft? is preferable 


Use sparingly 
[=MM cy/sec] 
go =980.665 gal 
gal/min 


Synonym for slug 
gm/cm? (density) 


With prefixes: kg; mg 
gf/cm?, gwt/cm? 


Unit or combining 
orm 


gram atomic weight 
gram calory 

gram molecular weight 
grav [ =32.174 ft/sec?] 
Hartree unit 

hecto-, hect- [ =10?] 
hekto- 

henry 


hertz [ =cy/sec] 
horsepower 
hour 


inch 

inch-of-mercury 

inch-of-oil 

joule 

Kelvin degree 
(temp. difference) 


kilo- [ =103] 
kilogram-calory 
kilocalory 
kilogram of force 
kilogram-mole 
kilomega- [ =10*] 
kilomole 


kcal 

kef, kgwt 
kM 
kmole 
knot knot, kn 
lambert [=ca/x cm?] lam 
large calory see 
light-year It-yr 

line 

liter 


lumen 

lux [=lu/m?] 
magnetic pole 
magneton 

maxwell 

mega-, meg- [ =10] 
megamega- [ =10!2] 
megohm [=10* ohm] 
meter 

meter-candle [ =lux] 
meter-kilogram-second unit mksu 
mho mho 
micro-, micr- [=107§] “ 
micrometer eee 
micromicro- [ =107!2] uu 
micron [=10-* m] um, # 
mil [=1073 in.] mil 
mile mi 
mil-foot mil-ft 
milli- [=1073] m 
millimicro- [510-9] mu 
millimicron [=10~ m] mu 
million | =108) M 
minute (of arc) 


meg, Mohm 
m 


m-ca 


min 
, 


minute (of time) min 
mn 
mks electromagnetic unit 
mole mole 
molecule molecule 
month mo 
myria- [ =10*] myria 
neper [ =8.686 db] nep 
newton [ =kg m/sec?] new 
normal atmosphere eee 
number 

ohm 


mksm 


oersted 
ounce (force) 


phon 
phot [ =lu/cm?] 
pint 
pole, unit magnetic 
poise [ =dy sec/cm?] 
pound (force) 
pound of mass 
omadel 
pulse 
quart 
radian 
Rankine degree 
(temp. difference) 
revolution 


rhe [ =1/poise] 


Examples of use; 
comments 


See calory 
See mole 


Use sparingly: hm 
Variant of hecto- 


y/m 

With prefixes: mh 

hz sec [=cy] 

hp hr 

hr/day 

Astron. text and tables: 35 
in./sec 


j/mole deg 


Use when not ambiguous: 
erg/deg molecule 

kcal, kev, kmole 

See kilocalory 


kgf/m? (weight-density) 
See kilomole 
1 kMev=1 Bev 


knot hr, kn hr 
See kilocalory 
line/cm? 


With prefixes: ml, k 
lu hr, lu/watt 

lux sec 

See pole 

See Bohr magneton 


Mev, Mm 
meg/m, Mohm/m 


lux is preferable 
Use sparingly 
mho/cm 

usec, HV, HE 

See micron 

wut 


mi/gal 
ohm/mil-ft 
ma, ml (not ‘‘cc’’) 


Mgal/day; see mega- 
min/sec (ang. velocity) 
In tables, etc. 


Astron. text and tables: 5™ 


Use sparingly 

new/m?, new/m? amp turn 
See atmosphere, standard 
no./hr 

ohm cm 

In crowded diagrams, etc. 


dy/pole 


lb sec or Ibf sec (impulse) 
lbm/ft? (density) 
ft pdl 


pulse/sec 


rad/sec 


_Btu/Ib R° 


rev/min, rev/sec (“rpm" 
and “rps” only in 
crowded tables, etc.) 





SYMBOLS FOR PHYSICAL UNITS 


TABLE I.—Continued. 


Unit or combining Examples of use; 


comments 


rod 

roentgen mr 
rowland 

r-unit 
rutherford 
second (of arc) 


Synonym for roentgen 


d( In tables, etc. 

second (of time) 

Astron. text and tables: 10* 
slug/ft®? (density) 
Synonym for x-unit 

See calory 

in. 

statamp, statcoul 


slug [ =32.174 Ibm] 
Siegbahn unit 
small calory 


steradian 


“lu” is preferable to “Il” for ‘‘lumen.”’ Well- 
established exceptions to this rule include 
“ogamma’”’ for “gamma,” ‘j” for ‘‘joule,’”’ and 
“vy” for “‘volt.”” If the name of a unit contains 
only two or three letters, this name, rather than 
an abbreviation, should be used; examples are 
“bar,” “day,” “erg,”’ “lux,” and ‘‘ohm.” 

(2) The symbol preferably should consist of the 
first two or three letters of the name of the unit. This 
facilitates recognition and pronunciation, as in 
“dy” for “‘dyne,” “lu’’ for “lumen,” and “‘oer’’ 
for ‘‘oersted.”” Among the firmly established ex- 
ceptions to this rule are: “ct,” “ft,” “hr,” “hy,” 
“ib,” aad “ez.” 

(3) Preferably there should be only one symbol 
for a particular unit. However, it seems desirable 
to recommend certain exceptions; for example, 
both ‘‘a” and ‘‘amp”’ for “ampere,” and ‘“‘b” 
and ‘‘bel”’ for “‘bel,”’ the ‘‘a” and “‘b” to be used 
only with combining forms, as in ‘‘ma”’ and “‘db.”’ 

As the rule implies, the same symbol should be 
used for both singular and plural forms; thus, 
10 amp, not 10 amps. 

One of several objections to such symbols for 
secondary (derived) units as ‘‘fps” and ‘‘rpm”’ is 
that the “‘f,’’ ‘‘s,”’ “‘r,”” and “‘m” are used in place 
of the well-established ‘‘ft,”’ “‘sec,” “rev,” and 


(4) Periods should be omitted from symbols. An 
exception is ‘‘in.,”’ for ‘‘inch,” since omission of 
the period would often result in confusion. 

(5) The symbol for a combining form denoting a 


multiple or submultiple should be a single letter. 
Thus ‘“p”’ 


and so on. 


for ‘‘micro-,” ‘‘M” for ‘‘mega-,” 


Unit or combining Examples of use; 


Symbol comments 


stilb [ =ca/cm?] eee 
thousand [ =10*] k 
ton (force) ton, tont 
ton of mass tonm 
turn turn 
vibration yi vib/sec 
volt rv v coul, volt coul 
watt watt hr 

With prefixes: kw 
weber web/m? 
x-unit 1 xu =1.0020 mA 
week 
yard 
year 


Deprecated; use ca/cm? 
kBtu, kft 

ton mi or tonf mi (work) 
tonm mi (transportation) 


USE OF UNIT SYMBOLS IN PRINTED TEXT 


(i) Set the symbol for a unit in roman type. 

(ii) Use an abbreviated unit-symbol only when it: 
(a) is preceded by a numerical value, or (b) appears 
in headings of tables or in crowded text, in which 
cases the symbol is enclosed in parentheses. Thus: 
“25 cm”; “several centimeters’’; ‘‘volumes, in 
cubic centimeters”’; ‘‘volumes (cm')’’; ‘‘v (cm?).”’ 

(iii) Use standard signs to indicate all mathe- 
matical operations with unit symbols. Thus indi- 
cate: (a) multiplication by a space, a center dot, 
or even a X;; (bd) division by a solidus, a negative 
exponent, or the ordinary fractional form; (c) a 
power by a positive exponent; (d) a root by a 
fractional exponent or \/. Thus dy cm, dy-cm, 
or even dy Xcm, but not dy-cm; ft/sec, ft sec, 


ft : ; 
or —, not fps, fs, or (in equations) ft per sec; 
sec 


cm?, not cu cm or cc. 

A hyphen should not be used to indicate multi- 
plication, or a ‘‘p’’ to indicate division, since 
these devices are not used elsewhere in mathe- 
matics for these purposes. A hyphen should be 
used only when it is desired to separate parts of a 
single unit symbol; thus a hyphen is appropri- 
ately used in “‘ft-ca,” since the “foot-candle”’ 
does not represent the product of “foot’’ and 
“candle.” 

One report on “Abbreviations for Scientific 
and Engineering Terms’ advocates abbrevia- 
tions such as the following for secondary units: 
cfm (for ft?/min); kgps (for kg/sec); mphps (for 
mi/hr sec); psf (for lb/ft?); rpm (for rev/min). 


1 American Standards Association, document Z10.1-1941 
(American Society of Mechanical Engineers, 1941). 





296 


A form such as “rpm’’ may be used on crowded 
diagrams or apparatus labels; but from the stand- 
point of suitability for substitution in physical 
equations and ease of mathematical manipula- 
tion, the aforementioned forms are comparable 
to such unconventional and ambiguous algebraic 
expressions as: cab (for a*/b); apb (for a/b); 
apbpe (for a/bc); asb (for a/b?); and so on. Of 
course, forms such as “‘cfm’’ or ‘‘mphps’’ may 
be regarded merely as special names for second- 
ary units, in the sense that ‘‘dyne”’ and “‘erg”’ are 
names for such units. But the suggestion that 
special names, many of them unpronounceable, 
be given to all commonly used secondary units, 
metric and English, multiple and submultiple, is 
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certainly not a move in the direction of greater 
economy of thought and learning. 

This proposal has been submitted to the re- 
cently reorganized Committee on Terminology, 
Symbols, and Abbreviations of the American 
Association of Physics Teachers. Interested 
readers are urged to send criticisms and sugges- 
tions concerning the present recommendations 
to the committee. The committee consists of: 
J. D. Elder, H. K. Hughes, T. H. Osgood, D. 
Roller (chairman), L. D. Weld, and M. W. 
Zemansky. It is contemplated that this com- 
mittee will be assisted by a group of consultants, 
to be selected from among those who have an 
evident interest in contributing to its work. 


Teaching a Philosophy of Experimentation in a Course in Electrical Measurements* 


Francis T. WORRELL 
Rensselaer Polytechnic Institute, Troy, New York 


(Received September 12, 1952) 


A course in electrical measurements is described wherein the emphasis is shifted from the 
learning of various electrical measuring techniques to the development of an appreciation of 
experimental problems in general. This is accomplished by dropping the requirement of written 
reports, by allowing sufficient time for the performance of each experiment for thoughtful work, 
and by selecting the experiments for certain features of general experimental significance rather 
than for specific electrical measuring techniques. Examples are given of typical experiments. 


HE course in electrical measurements for 

physics students contains a number of 
pitfalls for the unwary instructor who is called 
upon to organize it. Not the least of these pitfalls 
is the temptation to start by gathering a large 
number of resistance boxes, galvanometers, po- 
tentiometers, bridges, and other standard elec- 
trical measuring instruments; then to select from 
the numerous textbooks and manuals available a 
number of standard experiments, and have the 
student run through them—the more experi- 
ments the better. At the end of a course of this 
sort the student has probably acquired a certain, 
possibly considerable facility with electrical 
equipment, and has learned a number of specific 
techniques. But if he has acquired no more than 
this, have we accomplished a desirable result? 
Have we done any more than a good trade 


* Presented before the American Association of Physics 
Teachers, Iowa City, June 12, 1952. 


school could have done? In other words, have 
we contributed to the “liberal” scientific educa- 
tion of the student? If the answer to these 
questions is in the negative—and that would be 
the author’s answer—then we should try to see 
what more could or should be accomplished in 
such a course. 

The subject material in a course in electrical 
measurements is relatively easy to grasp. The 
experimental conditions in most of the standard 
experiments are generally rather good (compared, 
for instance, with conditions prevalent in acous- 
tical or thermal measurements), and the tech- 
niques and instruments are well developed. A 
facility with the instruments is not difficult to 
acquire; further, electrical measurements are 
made so often in actual work in the research 
laboratory that the student, or graduate, will 
have ample opportunity to get the facility later. 
An appreciation of the potentialities and limita- 





TEACHING A PHILOSOPHY OF EXPERIMENTATION 


tions of different experimental techniques is less 
easily acquired, however, and is therefore worth 
more emphasis. If a student performs a number 
of well-defined experiments with almost sure 
results he may learn little about these matters, 
or about the philosophy of experimentation in 
general; and, after all, this really should be one 
of the principal objectives of our laboratory 
courses. 

The student taking a laboratory course in a 
subject which is somewhat difficult to under- 
stand exerts most of his effort in the comprehen- 
sion of the subject; but when the subject ma- 
terial is not difficult the instructor has the 
opportunity to acquaint the student with some 
of the mechanics of performing the experiments 
themselves, and to get him to appreciate more 
fully the reasons for the methods which are used, 
and some of the limitations of experimental 
work. The way to do this, in the opinion of the 
writer, is, in the first place, to plan to do only a 
small number of experiments, allowing plenty of 
time for each one, and, in the second place, to 
select as far as possible experiments which illus- 
trate general experimental procedures rather than 
specific electrical measuring techniques. A start 
has been made in this direction in the course 
which will now be described. 

The course in Electrical Measurements, a 
first-term junior course for students majoring in 
physics at Rensselaer Polytechnic Institute, is a 
two-credit course consisting of two three-hour 
laboratory periods a week.’ The first hour of 
the first period of each week has been used 
generally for a discussion of some basic ideas, 
and of some experimental methods not covered 
in the laboratory which are considered worthy of 
special mention. No attempt is made to cover 
the.field in the detail in which a textbook would, 
nor do the discussions follow the textbook used 
in the course. Rather, it is hoped that the student 
will get the general picture from his reference 
reading for the experiments which he performs. 

Reports on the experiments are not required. 
It is felt that the writing of reports, important 
as it is—and poorly done as it is, after all the 
practice which students get—would be done at 


‘This has been preceded by a four-semester course in 
general physics, and is concurrent with a course in circuit 
theory and a course in electric and magnetic fields. 
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the expense of the thorough experimentation 
which this course emphasizes. Instead, the stu- 
dent is asked to keep what the instructor hopes 
may come to bear a reasonable resemblance to a 
research notebook. The student is asked to use a 
bound notebook, since the looseleaf notebook, 
with all its convenience, encourages the discard 
of sections containing information thought to be 
incorrect, but sometimes later found to be of 
interest. He is encouraged to use his book not 
only as a place in which to record the usual data, 
but as a place in which to record passing observa- 
tions made during the course of the experi- 
ment, notes on changes in the originally proposed 
procedure, and such things; in other words, 
comments which indicate that he was thinking 
while he did the experiment. The student is 
asked to record data directly in the book as he 
proceeds, even if it means crossing out some 
mistakes and starting a tabulation over, rather 
than to take the data on loose pieces of paper 
with the idea of sifting these data before putting 
them in the book. The instructor freely admits 
to them that the book that results may not be as 
pretty in the first case as in the second, but 
points out that there may be some valuable 
information lurking in the unwanted data (the 
old story of the Ugly Duckling) which would be 
lost in the selection process. The student who 
wants a neat book will probably make more of 
an effort to use foresight rather than hindsight. 

In preparation for a given experiment before 
performing it in the laboratory the student is 
expected to read the instruction sheet? for the 
particular experiment, also the reference ma- 
terial listed there, and then write in the notebook 
a brief prospectus including diagrams, procedure, 
and equations to be used. (This material’ is not 
referred to as ‘‘Purpose’’ and ‘Procedure,’ be- 
cause this induces a conditioned reponse in the 
student, and he sees in his mind’s eye a long 
series of reports reaching back into the past. We 
want to avoid an association with what may 
seem to be a routine report which, if filled out 

?In the past the instruction sheets have been quite 
complete in themselves, as far as instructions were con- 
cerned, but a recent conversation with one of my colleagues 
has led me to believe that fewer instructions in advance 
would make for more careful planning on the part of the 
student. The detailed instructions could be available at 


the site of the experiment. In case there is no good refer- 
ence, the sheets must still be detailed. 





298 


in the proper way, without need for thought, 
will pass.) It is emphasized at the beginning of 
the term that he is merely doing what anyone 
should do in preparation for tackling an experi- 
mental problem 

When the student arrives at the laboratory 
his preparation is evaluated by an inspection of 
the prospectus which he has written in his 
notebook. Suggestions are made at that time for 
improvements, and a grade is given. This grade 
will be considered in giving a grade on the report. 

Also considered in the over-all grade is the 
performance of the student in the laboratory. 
The instructor observes his general conduct, the 
efficiency with which he performs the experiment, 
his ability to observe experimental difficulties 
along the way, etc. Occasionally a student will 
be asked a question during the experiment, 
usually on some point which might be overlooked 
by the thoughtless student. 

The student is asked to make some semi- 
quantitative statement as to the precision of 
the result. In his thinking he must consider (1) 
instrumental and, occasionally, statistical errors; 
(2) observational errors resulting from limita- 
tions in the readability of the instruments; (3) 
the sensitivity of the method. In the first item 
the students have already been drilled in their 
general laboratory, but in the others they have 
had no practice, and their evaluation requires 
more mature thought. 

.When the student has made the measure- 
ments, he of course calculates the results and 
evaluates them, including in his evaluation con- 
sideration of the precision, as mentioned above. 
He is encouraged to respect his results, and not 
to worry about the “correct’’ value. In cases 
where he measured an “unknown,” it remains 
“unknown,” except in so far as he believes his 
results. He is not told what he should get. In one 
experiment, where the dielectric constant of air 
is measured, he does have a handbook value 
with which to compare his results. Since the 
results he gets are generally in disagreement with 
the accepted value he is asked to see if he can 
explain the deviation in terms of certain experi- 
mental conditions. 

All this is intended to encourage him to take 
the attitude that his meters have told him 
something, and that it is his responsibility to 
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decide what they have told him. Have they 
told him what he thinks they did, or something 
else; that is, is the experiment valid? If it dis- 
agrees in its results with some other determina- 
tion this does not mean that the student’s work 
is ‘“‘wrong,” but that the underlying experi- 
mental conditions were not suitable. The stu- 
dent is allowed one week for each experiment, 
which means five or six hours in the laboratory, 
depending on whether some time is used for 
general class discussion. The experiments them- 
selves are mostly quite standard, but it is the 
way in which they are used which the author 
wishes to emphasize. 

In illustration of some of the points which 
have been discussed so far, let us look at some 
of the experiments which are performed in the 
laboratory, confining our discussion only to the 
particular features which are illustrative. The 
measurement of resistivity using the Kelvin 
double bridge provides two points of interest. 
First, the student, if he evaluates his precision 
correctly, will find that the principal source of 
error is not the measurement of the length of the 
test bar (about 1 meter) with a meter stick but 
rather the measurement of the diameter (about 
+ inch) with a micrometer. Second, the good 
student may have some time to check his 
measurement using a type-K potentiometer. He 
soon comes to realize that although the type-K 
is potentially more precise than his bridge, his 
results obtained with the potentiometer are 
poorer because of unsteady current through the 
unknown resistance. 

In one experiment a Weston Model 45 am- 
meter is calibrated against a Leeds & Northrup 
Student potentiometer. The student is asked 
whether it would have been better to calibrate 
the meter using (a) another Model 45 meter in 
series with it, or, (b) a type-K potentiometer. 

In another experiment a 1-yf condenser is 
measured in a Wien bridge, then in a Maxwell 
bridge. While the use of the latter to measure 
capacitance in terms of inductance is admittedly 
the reverse of the usual procedure, no harm 
seems to have resulted. The use of the two 
methods provides the student with the oppor- 
tunity of comparing different methods of making 
the same measurement. He is asked to compare 
the two methods with regard to the standard 
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used in each case (a stepwise-variable capacitor 
vs a continuously variable inductor), the ease of 
balancing, and the precision. In connection with 
the last item, he is asked to determine the 
sensitivity of the bridge, i.e., the minimum 
detectable unbalance. He is also shown how to 
improve his balancing procedure by interpolating 
between two points of equal unbalance, rather 
than by trying to find the null directly.* For 
detectors we firmly stick to the headphone, 
praising the detecting qualities of the ear. This 
is sometimes over the loud protests of the 
gadgeteers who demand nothing less than an 
oscilloscope. 

In another bridge experiment the frequency 
bridge is used to measure the variation of 
inductance of a variable inductor with frequency. 
The student is to plot his results in the form of 
curves of apparent inductance vs frequency for 
three chosen scale readings of the inductor. 
Since the capacitor used as a reference in this 
bridge has a stepwise adjustment it is not 
possible, in general, to achieve a balance at the 
chosen scale setting. Instead, one must plot 
inductance as a function of scale setting for 
several different frequencies, then from these 
curves obtain data from which to plot the desired 
curves. This is an important technique for them 
to know. 

In the measurement of the volt-ampere charac- 
teristics of rectifiers several are used, including 
a vacuum diode. A test of the student’s alertness 
is to see whether in the case of the diode he 
records unthinkingly that zero applied voltage 

*In this connection we can point out the applicability 


of this technique to such things as photometric measure- 
ments. 
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gives zero current. You'd be surprised how many 
do! 

In the calibration of a Hewlett Packard audio 
oscillator the student is presented with the diffi- 
culty that the only available laboratory standard 
is another Hewlett Packard oscillator no more 
precise than the one to be calibrated. This 
difficulty is resolved by setting the reference 
oscillator to zero beat with the 440-cps note 
from radio station WWYV, then calibrating the 
unknown in the range 44 cps to 4400 cps by 
using Lissajous figures on the screen of an 
oscilloscope. To extend the range of calibration 
upward the unknown is set with care at 4400 cps, 
the reference brought up that frequency, and 
the calibration continued. This equipment thus 
exhibits two techniques of value. The latter 
technique, for instance, has been used in the 
calibration of an optical pyrometer. 

The six experiments which have just been 
discussed represent, of course, only part of the 
list. The intent has been not to present a detailed 
discussion, but to point out the aspects of some 
typical experiments which are considered of 
importance in the course as we teach it. 

In summary, the laboratory in electrical 
measurements which has been described is in- 
tended to give the student an appreciation of 
the general problems of experimentation, rather 
than a knowledge of specific techniques. To this 
end ample time is given for the performance of 
each experiment, and the general nature of the 
experimental method is emphasized. Experi- 
ments are chosen on the basis of their general 
experimental value rather than on the basis of 
any specific technique. 


About 1904 I was asked by Rutherford to make a sensitive small-capacity electroscope the gold 
leaf of which would remain charged for two or three days. This I failed to do. So Rutherford said: 
“Lester Cooke used to make them; why can’t you? Get Jost the mechanic to make you onel"’ So I 
went to Jost and repeated this. He said: ‘If I could not make a better electroscope than Cooke, I'd 
shoot myself!" So he made a beauty to look at, but a bad one to go. Its leaf collapsed in twelve hours. 
This puzzled me. One night I could not sleep and got up in my diggings and made an electroscope 
of a tobacco tin, an amber mouthpiece of a tobacco pipe, and some Dutch metal foil; charged it with 
sealing wax and went to sleep. The leaf of this home-made freak electroscope remained open and 
charged for three days, and solved the problem. An electroscope made of material outside the labora- 
tory would remain charged for 48 hours inside the laboratory, but all the material inside the building 
was contaminated and coated with active deposit including the slow period transformations of 
radium. . . . Precautions were then taken to prevent the escape of radon which had been at the 
root of the mischief, hitherto unsuspected.—A. S. Eve. 
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The curriculum of the Oak Ridge School of Reactor Technology includes a laboratory course 
in reactor physics. The first experiments are concerned with the fundamentals of nuclear 
physics and the later ones are of direct relevance to reactor engineering. The laboratory manual 


of the course is available as AECU-2164. 


HE Oak Ridge School of Reactor Tech- 

nology (ORSORT) presents a year’s basic 
training in reactor engineering. The curriculum 
consists of about 38 weeks of classroom work 
including 26 weeks of laboratory, followed by 12 
weeks of either in-service training or participa- 
tion in a reactor design study. Lectures are 
presented in reactor theory, mathematics, reactor 
materials, reactor physics, reactor chemistry, 
health physics, reactor technology, and reactor 
engineering. Students are either on leave from 
industries associated with The Atomic Energy 
Commission’s work or are recent graduates in 
physics, chemistry, or any of the various branches 
of engineering. 

The large variation in the scientific back- 
grounds of the students presents peculiar educa- 
tional problems. In particular, the usual engineer 
has little or no knowledge of nuclear phenomena. 
Thus the first few experiments in the laboratory 
and some of the lecture material in the reactor 
physics course are not essentially new to most 
of the physicists and chemists. However, since 
the majority of students are engineers, it is 
necessary at present to include this material in 
lectures and laboratory. It is hoped that in the 
future much of the introductory laboratory ma- 
terial can be included in the students’ under- 
graduate training. 

The laboratory course at ORSORT is the 
product of the work of a number of people at 
Oak Ridge National Laboratory who have de- 
signed the experiments and instructed in the 
course. In particular, we should like to acknowl- 
edge the work of Dr. E. C. Campbell who 
initiated this laboratory course in 1950. At the 
beginning of the training period in September 


* On leave from The College of Wooster, Wooster, Ohio. 


the students are given a set of instruction sheets 
covering the experiments for the year. These 
contain references, the pertinent theory and 
procedure of the experiment, and questions for 
discussion. It is neither desirable nor possible 
to give anything but a brief account of the 
experiments in this article. However, anyone 
who is further interested in these experiments 
may purchase a complete set of instruction 
sheets for the 1951-1952 course.! 

During the year the student performs some 
26 experiments. The first group of these deals 
with alpha- and beta-particles, gamma-rays, 
Geiger counters, simple and complex radioactive 
decay, and health physics measurements. Most 
of these experiments are or could be performed 
in any modern physics course. Following these 
are five experiments in chemistry and the re- 
mainder are largely group experiments performed 
at the Oak Ridge graphite reactor. The student 
body is broken into four laboratory groups of 
approximately 18 each. Wherever possible, the 
students work in pairs on individual pieces of 
apparatus, though for many experiments as- 
sociated with the reactor this is not possible. 

The instructor holds a one-hour discussion on 
each experiment before it is performed. Questions 
on the previous experiment may also be discussed 
at this time. After the experiment is performed 
the student hands in a report consisting of the 
instruction sheets, the original data taken in the 
experiment together with appropriate calcula- 
tions and graphs, a discussion of the errors 
involved in the experiment, and answers to the 
questions. After the first experiment has been 


1 AECU-2164, “‘A Laboratory Course in Reactor Phy- 
sics,’’ Office of Technical Services, Department of Com- 
merce, Washington 25, D.C.; price $1.05. 
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completed, including the report, each student is 
given a copy of a sample report of this experi- 
ment. made by an instructor. This is done to 
indicate to the student the degree of complete- 
ness and detail required in the reports. In general, 
this has proved a very useful device both for the 
student and for the instructor. 

The first experiment is given largely to ac- 
quaint the student with the build-up and decay 
of radioactive nuclei produced by neutron 
bombardment of a small piece of silver or 
indium. The activation is done in a “neutron 
howitzer” consisting of a Po-Be neutron source, 
emitting about 107 neutrons per second, placed at 
the center of a cubical paraffin block two feet on 
edge. By placing the sample about one inch from 
the source an irradiation time of ten minutes 
sufficiently activates the sample so that it can 
be counted with the usual Geiger counter and 
scaling equipment. 

The subjects of the next three experiments are 
the specific ionization due to alpha-particles, the 
absorption and range of beta-particles, and the 
absorption of gamma-rays. In the alpha-particle 
experiment the ionization produced in a shallow 
ionization chamber is measured by a vacuum 
tube electrometer as a function of the distance 
from the alpha-source. The results are then 
plotted, giving a Bragg curve. The beta-particle 
experiment is concerned with measuring the 
intensity of beta-particles from P® after passing 
through various thicknesses of aluminum. An 
end-window Geiger counter is used to measure 
the intensity of the beta-rays. By plotting the 
logarithm of the intensity against the absorber 
thickness the maximum beta-particle range is 
found from the curve. From the empirical rela- 
tionship between the maximum energy and range 
given by Glendenin, the maximum energy of the 
beta-particles emitted by P® is determined. In 
the next experiment the absorption coefficients 
of gamma-rays from a Co® source are measured 
for lead and steel. The gamma-ray intensity 
transmitted through various thicknesses of lead 
and steel is measured by means of a Lauritsen 
electroscope. Absorption coefficients are obtained 
from the slopes of the semilog plots of intensity 
against thickness of absorber. 

The next experiments are concerned with 
typical health-physics measurements. In the 
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first of this group of three experiments, the 
@udents calibrate a beta-gamma survey meter, 
using a calibrated radium source. The dose rate 
in milliroentgens per hour is calculated for 
various distances between the source and survey 
meter and the corresponding meter readings are 
taken. A curve is then plotted of meter readings 
against the dose rate for each of the three scales 
on the instrument. The beta-gamma survey 
meter consists of a portable Geiger counter and 
appropriate electronic equipment using a count- 
rate circuit having three scales corresponding to 
maximum doses of about 20 mr per hr, 2 mr per 
hr, and 0.2 mr per hr. The second experiment 
consists in determining the amount of radium in 
an unknown source by comparing it with a 
known source. A Lauritsen electroscope is used 
as the measuring device. In the final experiment 
in this group the students use a variety of in- 
struments to survey the gamma-ray flux and 
the fast and slow neutron fluxes around the Oak 
Ridge graphite reactor. An ionization chamber 
type of instrument called a ‘‘cutie-pie”’ is used to 
measure gamma-ray intensities. For the thermal 
neutron flux a pair of Lauritsen electroscopes is 
used, one of which has its chamber coated with 
boron. In general, measurements are made of 
thermal neutron fluxes in the presence of a 
fairly large gamma-flux. The difference in the 
readings of the coated and uncoated instruments 
gives a measure of the thermal neutron flux. A 
calibration curve is supplied with the instruments 
so that it is possible to measure thermal neutron 
fluxes from zero to twice the maximum permissi- 
ble dose of 1750 thermal neutrons per cm? per sec. 
The fast neutron flux instrument consists of a 
paraffin-lined proportional counter in which the 
proton recoils from the paraffin are ‘measured. 
This instrument is designed to detect and 
measure approximately the fast neutrons in the 
0.2 to 10 Mev energy range and from 10 to 
10 000 fast neutrons per cm? per sec. 

In the next experiment the students have their 
only experience during the course with elec- 
tronics as such. They wire a Higinbotham 
scale-of-two circuit and test it for correct 
operation, examining operating voltages and 
wave forms of the voltage at various points in 
the circuit. Finally the individual units are 
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placed in series to simulate a scaling unit of 512, 
1024, etc. z 

The experiment on Geiger-Mueller counters is 
designed to familiarize the students with some of 
the characteristics of these counters and also to 
give a brief introduction to the statistics of 
counting. Natural uranium samples are used as 
sources. The resolving time of the counter and 
associated circuit is found by the two-sample 
method. A plot of counting rate vs tube voltage 
is also obtained. Finally, the number of counts 
occurring in a short time interval is observed 
over a large number of time intervals, and a 
statistical analysis is made of the observed 
distribution. 

In the experiment on the analysis of a complex 
radioactive decay curve, the half-lives and initial 
activities of two independent radioactive sub- 
stances are measured. The mixed substances are 
chosen so that one has a relatively short half-life 
of a few minutes while the other has a relatively 
long half-life of the order of an hour. The re- 
sultant activity, measured with a Geiger counter 
and scaler, is plotted on semilog paper against 
the time after removal of the sample from the 
reactor. By subtracting (‘‘peeling off’’) the ac- 
tivity due to the long half-life substance from 
the total activity, the activity of the short 
half-life substance is found. 

The next five experiments are concerned with 
chemical techniques and processes of importance 
to reactor engineers. The first of these involves 
the separation of metallic elements by solvent 
extraction. A single-stage water-benzene extrac- 
tion is used to separate a mixture of aqueous 
zirconium and thorium perchlorates. Radioactive 
tracer techniques are used for following the 
thorium and zirconium through the extraction 
cycle. In the next experiment, an ion-exchange 
separation of cobalt and nickel is performed. An 
ion-exchange column containing a strong base- 
type anion resin is used. The cobalt is followed 
with a radioactive tracer, while the presence of 
nickel is determined by chemical test. The third 
experiment in this group is concerned with 
radioactivation analysis. The manganese content 
of an aluminum alloy sample is determined, 
first, by calculation, using the known neutron 
flux, cross section, etc., and the measured activity 
of the sample, and, second, by comparing the 
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activity of the alloy sample with that of a 
sample of pure Mn. In the first method, a 
calibrated high-pressure gamma ionization cham- 
ber is used for obtaining the sample activity. 
In the latter method the samples are dissolved, 
then evaporated on watch glasses in a standard- 
ized manner, and counted with an end-window 
Geiger counter. The next experiment deals with 
radiation chemistry. The build-up of pressure 
over a distilled water sample is observed as a 
function of time of exposure to gamma-rays 
from a Co® source. The experiment is then 
repeated using a water solution of potassium 
bromide. Next, a crystal of sodium chloride is 
irradiated by gamma-rays at liquid nitrogen 
temperature. The crystal is allowed to warm 
slowly, and the release of stored energy in the 
form of visible light is observed by means of a 
photomultiplier tube and associated circuit. The 
final chemistry experiment consists of deter- 
mining the neutron flux in the Oak Ridge 
graphite reactor by measuring the amount of 
fission product Ba™® in a natural uranium sample 
that has been irradiated in the reactor for a 
known time. The radioactive barium is separated 
from the uranium chemically, using natural Ba 
as a carrier. The amount of Ba"® present is 
determined by counting with a calibrated end- 
window Geiger counter. As the fission yield, 
decay constant of Ba"™®, and the U* fission cross 
section are known, the neutron flux can then be 
calculated. The flux is also measured by irradiat- 
ing an aluminum alloy sample of known Mn 
content. The activity of the sample is deter- 
mined with a calibrated high pressure gamma 
ionization chamber, and the neutron flux com- 
puted. 

For the experiment on neutron scattering and 
absorption, the neutron howitzer is used together 
with boron-trifluoride proportional counters. The 
back scattering of neutrons by paraffin blocks is 
measured as a function of the number of blocks, 
while the total cross section is measured for 
aluminum, iron, and paraffin by the transmission 
method. Relative scattering cross sections of 
carbon, aluminum, and lead are determined by 
scattering a beam of thermal neutrons from the 
graphite reactor through 90° by the different 
substances and measuring the scattered beam 
with a BF; proportional counter. 
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The experiment on the measurement of ther- 
mal neutron capture cross sections gives the 
students an opportunity to operate the ‘‘pile 
oscillator,’” a piece of research equipment regu- 
larly used at the Oak Ridge graphite reactor. 
With this apparatus the students measure the 
thermal neutron capture cross sections of various 
materials relative to that of gold. They also 
measure the absolute ‘‘capture areas’’ of various 
cadmium and boron samples. 

In the next experiment the students obtain 
beams of monochromatic neutrons by crystal 
diffraction. By means of transmission measure- 
ments a set of boron absorbers is calibrated for 
use in a later experiment. Since a beam hole at 
the Oak Ridge graphite reactor is available to 
the school, it is possible to set up a simple 
neutron diffraction apparatus for the exclusive 
use of the students instead of having them 
perform the experiment on apparatus ordinarily 
used for neutron diffraction research. In this 
manner the essential features of the apparatus 
are made clear at a not too serious sacrifice in 
the quality of the data obtained. 

The calibrated boron absorbers are used in 
the next experiment to measure the neutron 
energy at which the main resonance in indium 
appears. This is done by measuring the activa- 
tions of indium foils by a cadmium-filtered neu- 
tron beam from the Oak Ridge graphite reactor 
with and without one or more of the boron 
absorbers. Use is made of the fact that the 
absorption cross section of boron varies as the 
reciprocal of the neutron velocity. 

For the experiment involving measurement of 
the number of neutrons emitted per thermal 
fission and the thermal fission cross section of 
U5 the students are given no detailed instruc- 
tions. Rather they are provided with the equip- 
ment necessary for the performance of the experi- 
ment, and are expected to work out the detailed 
experimental procedure themselves. The equip- 
ment and facilities provided consist of a U2 
sample, a collimated neutron beam from the Oak 
Ridge graphite reactor, a calibrated parallel- 
plate fission chamber and associated equipment, 
a calibrated hydrogen-recoil fast neutron counter 
and associated equipment, cadmium shields, 
indium foils, a ‘‘standard pile’ in which the 
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neutron flux is known, and an end-window 
Geiger counter. 

In the next experiment the students measure 
the half-lives associated with the delayed neu- 
trons emitted following thermal neutron induced 
fission in U*, This is done using the “fast 
pneumatic,” a piece of research apparatus at 
the Oak Ridge graphite reactor which enables 
one to remove a sample from the reactor and to 
start counting it in a fraction of a second. The 
delayed neutrons are detected with a suitably 
shielded BF; counter, and the apparatus plots 
directly the logarithm of the counting rate as a 
function of time. From this record the students 
find the half-lives of the various neutron groups 
using the ‘‘peeling off’ technique learned in an 
earlier experiment. 

Next the students measure the diffusion length 
of thermal neutrons in graphite by exposing a 
series of indium foils along the vertical axis of a 
rectangular parallelepiped of graphite which con- 
tains a Po-Be neutron source. The saturated 
activities of the foils are computed for distances 
from the source at which the neutron flux is 
nearly pure thermal, and the logarithm of the 
saturated activity is plotted against the distance 
from the source. From the slope of the straight 
line thus obtained and the known dimensions of 
the graphite block the diffusion length is cal- 
culated. 

In the next experiment, on the moderation of 
neutrons in graphite, a graphite stack containing 
a Po-Be neutron source is used. Cadmium- 
covered indium foils are inserted in the stack at 
various distances from the source. The foils are 
counted with Geiger counters and scalers, and 
the Fermi age of the neutrons of In resonance 
energy is then obtained from a curve of logarithm 
of saturated foil activity vs the square of the 
distance from the source. 

The next experiment consists of a study of the 
neutron flux distribution in a subcritical reactor 
or ‘‘exponential pile.’’ The pile used is an eight- 
foot cube duplicating the lattice of the Oak 
Ridge graphite reactor. This cube rests on a 
graphite base containing a Po-Be neutron source. 
Traverses are made vertically and horizontally 
using BF; counters. (The horizontal traverse is 
made parallel to the fuel channels.) A Cd- 
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covered rod, simulating a control rod, is then 
inserted vertically in the pile and the horizontal 
traverse repeated. The material buckling of the 
lattice is determined from the data of the vertical 
traverse, and is then used to estimate the size 
that the pile would have to be in order to be 
critical. 

The flux distribution perpendicular to the fuel 
channels in the Oak Ridge graphite reactor is 
studied in the next experiment. One end of a 
long copper wire is attached to a small block or 
“rabbit” which fits snugly into a tube running 
into the pile. The rabbit is drawn into the tube 
by the partial vacuum maintained in the reactor. 
After a suitable exposure time, the wire is reeled 
out. The end of the wire is then soldered to 
another copper wire which passes through a lead 
pig and on to a take-up reel which can be driven 
at constant speed by a motor. The lead pig 
contains a beta-proportional counter whose out- 
put is connected through an amplifier to a 
count-rate meter and speedomax recorder. After 
waiting about an hour for the decay of the short 
lived 5.1-minute Cu®, the motor is started and a 
record of the activity of the wire as a function 
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of its position in the pile is thus recorded by the 
Speedomax. 

The final experiment of the course at present 
is a study of the transient behavior of the Oak 
Ridge. graphite reactor during shutdown and 
start-up. The logarithm of the reactor power and 
the instantaneous reactor period are recorded 
graphically from the neutron-induced ionization 
in a gamma-compensated chamber. From the 
measured period on start-up the multiplication 
factor, k, (the factor by which the neutron flux 
changes from one neutron generation to another) 
is calculated for the supercritical reactor. After 
shutdown, the delayed neutrons constitute a 
neutron source and this enables one to calculate 
the multiplication factor of the subcritical reac- 
tor. The temperature coefficient of reactivity is 
also calculated from a temperature chart of a 
uranium slug. 

This year it was found possible to have 
students actually operate a reactor. The student 
body was divided into groups of six and each 
group spent two hours at the low intensity 
research reactor at Oak Ridge National Labora- 
tory. This formed a suitable culmination to both 
the lecture and laboratory courses. 


The Steady-State Theory 


The expansion of the universe, which can be inferred either 
from thermodynamics or from astronomical observations, 
would seem to lead to a thinning out of material. By the per- 
fect cosmological principle the average density of matter must 
not undergo a secular change. There is only one way in which 
a constant density can be compatible with a motion of expan- 
sion, and that is by the continual creation of matter. Only 
if the diminution of density due to the drift to infinity is 
counteracted by a constant replenishment of newly created 
matter can an expanding universe preserve an unchanging 
aspect. 

There is little doubt that the continual creation of matter 
necessary in this [steady-state] theory is the most revolu- 
tionary change proposed by it. There is, however, no observa- 


tional evidence whatever contradicting continual creation at 
the rate demanded by the perfect cosmological principle. It is 
easily seen that this is, on the average, 


3X (mean density of matter in the universe) 
X Hubble's constant = 10-* g/cm? sec 


approximately. In other words, on an average the mass of a 
hydrogen atom is created in each liter of volume every 10° 
years. . . . It is clear that it is utterly impossible to observe 
directly such a rate of creation. There is therefore no contra- 
diction whatever with the observations, an extreme extrapola- 
tion from which forms the principle of conservation of matter. 
—Cosmology, by H. Bonp1, Cambridge University Press, 
1952, p. 143, by permission. 
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The role that pulverization of the target material may play in absorbing the energy of an 
impacting missile is discussed. The energy that can be absorbed in this way will, in general, 
depend on the area of new surface that is created. Approximate calculations suggest that, for 
impacts against rock-like materials and glass, very considerable amounts of energy could be 
absorbed. The production of rock flour at the Barringer Meteorite Crater in Arizona is cited 


as an example. 


ANY processes can be operative during 
conversion to a stable situation of the 
transitory energies of impacting missiles and of 
exploding charges. It seems not unlikely that, 
in some cases, pulverization of a portion of the 
target material by the impact of the missile or 
the force of the explosion will play an important 
role. Such, of course, would not be the case for a 
perfectly plastic impact, viz., a steel bullet 
striking a lead target. It is well known, however, 
that glass and glass-like materials are broken 
up into very fine flour-like material when sub- 
jected to shock loading.! The crushing of glass 
under explosive attack has been discussed re- 
cently by Pugh and his co-workers?* and by 
Kolsky and his co-workers.4> The Barringer 
Meteorite Crater in Arizona contains millions of 
tons of rock flour. Barringer® has pointed out 
that it is almost certain that this rock flour is 
sandstone that was pulverized at the time the 
meteorite impacted. Considerations that will be 
discussed here suggest that in the foregoing 
cases, storage of energy on the fracture surfaces 
of the comminuted material may be an important 
way by which energy is absorbed by the target. 
The amount of energy E, that can be con- 
verted by this process will depend upon (1) the 
total amount of new surface that is generated by 


1John S. Rinehart, Meteoritical Society, 
Astron. .58, 9 (November, 1950). 

2? Pugh, Heine-Geldern, Foner, and Mutschler, J. Appl. 
Phys. 22, 487 (1951). \ 

3 Pugh, Heine-Geldern, Foner, and Mutschler, J. Appl. 
Phys. 23, 48 (1952). 

* Kolsky, Snow, and Shearman, Research 2, 89 (1949). 

5D. G. Christie and H. Kolsky, Trans. Soc. Glass Tech- 
nology 36, 65 (1952). 

6 BM. Barringer, Paper read before National Academ 
of Sciences at its Autumn Meeting at Princeton Uni- 
versity, November 16, 1909. 
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the impact of the blow, and (2) the energy that 
is required to form each unit area of this surface. 
Analytically, these facts may be expressed by 
the equation 


E,=<Ae, 


where A is the total surface area and a is the 
energy stored per unit surface. In the case of 
liquid droplets, o is equal very nearly to the 
surface tension of the liquid. 

It is convenient physically to think in terms 
of the mass of material that will be pulverized. 
If a is taken as the new area generated per gram 
of material comminuted, the energy stored on the 
surface of a gram of material will be given by 
the equation 


é,=do. 


If now we assume that a projectile of mass m and 
velocity v strikes a target, the mass M of target 
material that will be crushed can be assumed to 
be given roughly by the equation 


u=(™)«, 


where K is a factor that represents the efficiency 
of the conversion of the kinetic energy of the 
projectile to surface energy. Absorption of energy 
by the projectile itself will not be considered 
here although it could play an important role in 
many Cases. 

The primary problem in deciding how import- 
ant comminution may be as an energy absorbing 
process is, from a physical point of view, one of 
evaluating e, and K. First, let us consider the 
factors upon which K depends. Some 80 years 
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ago, Rittinger’'* stated that the useful work 
accomplished in crushing and grinding of ores is 
directly proportional to the new surface area 
produced. Most investigators who have discussed 
comminution problems since that time have 
assumed that most of the energy goes into heat- 
ing up the rock and that only a few percent goes 
into breaking it up. Thus in normal crushing 
operations, the amount of energy stored on the 
surface is small compared with the work done in 
crushing the material. However, it seems not 
unlikely that during high speed impact or during 
impulsive loading of the type produced by the 
detonation of an explosive charge, a much larger 
fraction of the energy will go into surface energy. 
It is apparent that K will depend strongly on 
the physical properties of the target material. 
If the target material flows plastically, then K 
might very well be equal to zero. It is to be 
remembered in this connection that the target 
material may behave very differently under the 
extreme dynamic loading conditions that exist 
during comminution than it would under static 
loading. For materials that seem to behave in 
a brittle fashion under impact, such as glass and 
crystalline quartz, K may be quite large. Since 
in a truly brittle material no plastic flow can 
occur, K might approach a value nearly equal to 
one. In so far as the author knows, no quantitative 
determinations of the value of K have been made 
so that whatever numerical value is assigned, 
K is at best an estimate. Indeed certain investi- 
gators feel that, under the high pressures created 
by the impact, glass will flow plastically so that 
the value of K will be negligibly small. The 
author does not hold this point of view. 

In order to make a realistic estimate of the 
value of e, two factors must be considered: (1) 
the energy stored per unit area of surface, and 
(2) the distribution in shapes and sizes of the 
comminuted material. The exact value to use 
for o in a given case is not clear. In the case of 
quartz, investigators® concerned with comminu- 
tion problems often take a value of about 900 


7Fred C. Bond and Jen-Tung Wang, Trans. Am. Inst. 
Mining Met. Engrs. 187, Mining Engineering (1950). 

8 Fred C. Bond and F. E. Briber, Jr., Pit and Quarry 
44, 173 (1950). 

® A. M. Gaudin, Principles of Mineral Dressing (McGraw- 
Hill Book Company, Inc., New York, 1939). 
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ergs/cm?*. Griffith’? used 520 ergs/cm? in con- 
nection with his studies of fracturing in glass. 
Recent work" in connection with silicosis pro- 
duced by quartz particles indicates the existence 
of a highly soluble surface layer on these par- 
ticles. An experimental determination of sur- 
face energy is extremely difficult because a 
freshly created surface will become contaminated 
almost at once. In brief, the whole problem of 
the amount of energy represented by a square 
centimeter of surface on a quartz or glass particle 
can be only very roughly approximated. It seems, 
however, not unreasonable to assume that the 
energy probably ranges from about 200 ergs/cm? 
to about 900 ergs/cm?. 

Similarly, distribution of shapes and sizes of 
particles is also not known exactly for particles 
whose dimensions are less than a few microns. 
Indeed, accurate size distribution measurements 
below 10 microns are rarely accomplished and are 
usually of doubtful validity.? The rock flour 
found in the bottom of the Barringer Meteorite 
Crater is an excellent example of crushing due to 
impact. Merrill'? has described this rock flour as 
having a sharp gritty feeling when rubbed be- 
tween the fingers. Under the microscope, the 
individual particles of the flour are found to be 
sharp angular bits of quartz that must have come 
from shattering of individual grains of sand. 
There are millions of tons of this rock flour in 
the bottom and along the sides of the crater. 
We found that about 30 percent of one sample 
of flour passed a 325 mesh (43 micron) screen. 
Gaudin® estimates that particles as fine as 0.001 
micron might be produced by comminution. 
Rosenthal," in connection with other studies, has 
recently examined 0.05-micron marble particles 
that were produced by crushing. As regards sur- 
face area per unit mass, Gaudin® suggests that 
the area of an individual particle is likely to be 
about 1.7 times that of a cube or sphere of 
equivalent mass. 

It is possible to form a very rough estimate of 


10 A. A. Griffith, Trans. Roy. Soc. (London) A221, 163- 
198 (1920). 


ul a. Gordon, and Griffin, Nature 169, 538- 


540 (1952). 
( one P. Merrill, Smithsonian Misc. Coll. 50, Part 4 
1 i 
1D. Rosenthal and M. Kaufman, J. Appl. Phys. 23, 
600 (1952). 





MODE FOR ENERGY ABSORPTION DURING IMPACT 


the amount of material that would be crushed by 
a particular missile, such as a large meteorite, if 
we make some very simplifying assumptions. 
As an example, assume (1) that the target ma- 
terial is broken up into cubes all one micron on a 
side, (2) that the energy stored on the surfaces 
is equal to 500 ergs/cm*, and (3) that K, the 
conversion coefficient, is equal to one. We find, 
with these assumptions, that 1.2107 ergs of 
energy could be stored on the surfaces of one 
gram of comminuted target material of density 
2.5 g/cm’, i.e., about that of quartz. It is inter- 
esting that experimentally! it has been found 
that small caliber ammunition will on impact 
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lose energy at the rate of about 2.6X107 ergs 
per gram of rock crushed. 

It is recognized, of course, that the above 
treatment is neither rigorous nor exact. None- 
theless, it is believed that storage of energy on 
surfaces of comminuted material may be of much 
significance as regards the partition of energy 
when large high velocity missiles such as meteor- 
ites form craters and that further pursuit of the 
problem might be profitable. 

In conclusion, appreciation is expressed to 
Dr. H. H. Nininger who kindly supplied a sample 
of rock flour from the Arizona crater and to 
Ensign W. G. Feuerstein who examined the flour, 


An Undergraduate Experiment for the Determination of Alpha-Particle Range 


WaLtTER C. MICHELS AND Eva WIENER* 
Bryn Mawr College, Bryn Mawr, Pennsylvania 
(Received August 8, 1952) 


Simple apparatus, using a scintillation counter with a 931-A photomultiplier, is described. 
The range of alpha-particles from Po is determined with this apparatus by varying the pres- 
sure in a shielded container, keeping the geometry constant. Data are presented on a 0.06 
millicurie source, yielding a range of (3.81-0.03) cm at 760 mm Hg and 15°C. This value is 
obtained after correction for finite source thickness. The experiment offers particular value in 
giving students an opportunity to consider statistical fluctuations and to analyze the effects 


of finite source thickness, straggling, etc. 


ECENT developments in connection with 
the simple and low cost scintillation counter! 
make feasible the increasing use of radioactivity 
and nuclear experiments in the undergraduate 
laboratory. The apparatus and the experiment 
described here offer one such possibility. 

In order that the range of alpha-particles may 
be determined with accuracy, it is desirable 
that the geometry of the source-detector system 
be kept constant. This is accomplished in the 
apparatus shown in Fig. 1 by the use of a rigid 
chamber which may be evacuated to vary the 
“effective path” of the alpha-particles without 
displacing either the source or the detector. 
Since the container has been constructed to pro- 
vide shielding for a variety of experiments, it is 

* Now at Stanford University, Stanford, California. 

1For a discussion of these developments and an ex- 
tensive bibliography, see W. H. Jordan, “Detection of 


Nuclear Particles,’’ Annual Review of Nuclear Science 
_— Reviews, Inc., Stanford, 1952), Vol. 1, pp. 207- 
226. 


considerably heavier than is necessary for alpha- 
particle work. It consists of two coaxial cylin- 
drical shells of Cerrobend,? separated by a one 


pare & manometer 


ora 
a RZ 


Cerrobend N 
Poraffin N 
N 
N 


OLA 


Scaler Amplifier 


Fic. 1. Apparatus for determination of 
alpha-particle range. 


*Cerrobend is a commercially available material pre- 
pared by Cerro de Pasca Company. Its composition is Bi 
50.0 percent, Pb 26.7 percent, Sn 13.3 percent, and Cd 
10.0 percent. The bismuth and lead offer high mass ab- 
sorption for beta- and gamma-radiations, while the cad- 
mium guarantees absorption of slow neutrons. The ma- 
terial melts at 70°C, so that it can be cast with great ease. 
It also machines readily. 
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centimeter thickness of paraffin. The cover is 
similarly constructed. 

Provision is made, as shown, for cooling the 
931-A photomultiplier, which is operated at 
980 v and with 10° ohms across each stage. Both 
batteries and a rectifier have been used to 
supply the potential to the photomultiplier, but 
it has been found extremely difficult to filter the 
rectifier output sufficiently to reduce the back- 
ground count to the level which is achieved with 
a battery supply. The output of a two-stage 
R-C coupled amplifier, operating at a voltage 
gain of approximately 60, is fed into a Berkeley 
scaler. When weak sources are used, the scaler 
could be replaced by a single decade unit and a 
high speed mechanical counter. 

We have found zinc oxide?’ to be a satisfactory 


Equivalent Distance (cm at76cm & 15°) 
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Oo 
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‘ 
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Pressure (cm Hg) 27.3cm Hg 


(3.81 + 0.03)cm 
Fic. 2. Analysis of the data. 


phosphor for alpha-particle detection. The dis- 
tance between the polonium foil source and the 
ZnO is 9.80 cm. The “equivalent distance’’ in 
air, under the standard conditions of pressure 
P, and Kelvin temperature T> is therefore 


y= (9.80 cm)(PT0/PoT), (1) 


where P and T are, respectively, the pressure 
and Kelvin temperature in the chamber. The 
pressure is varied either by the use of a controlled 
leak or by variation of the pumping speed. The 
temperature correction is particularly important 
when the photomultiplier is cooled. With dry ice 
and acetone used as a coolant, the chamber of 
our apparatus runs at —1.5°C. 


3 Fluorescent ZnO supplied by J. R. de Vore of the New 


Jersey Zinc Company. This is an experimental phosphor 
labeled CP-215. 


MICHELS AND E. WIENER 


With our arrangement, the phosphor subtends 
an angle of about 7X10 steradian. The effi- 
ciency of the counter, with the “light piping 
system” shown, is about four percent. The count- 
ing rate is therefore approximately 210-5 of 
the disintegration rate of the sample. The photo- 
multiplier used in this work, a typical 931-A 
tube, shows a background count of about 100/ 
sec at room temperature and about 40/sec at 
dry ice temperature. Taking the minimum prac- 
tical source as one which will give a counting 
rate equal to the background of the cooled sys- 
tem at low pressures, we find that this source 
must exhibit about 10° disintegrations per sec- 
ond, i.e., must have an activity of about 0.05 
millicurie. The upper limit of source strength is, 
of course, determined by the resolution time of 
the amplifier and the scaler. 

The experimental points shown in Fig. 2 indi- 
cate the net counts observed at various pres- 
sures, using a 0.06 millicurie polonium foil 
source, and the uncertainties as determined from 
the statistics. The scale shown at the top of the 
curve represents the reduction of the observed 
pressures to equivalent distances at 760 mm Hg 
and 15°C, as determined by Eq. (1). 

It is clear that the data shown are typical of a 
thick source. In order that these data might be 
reduced to those for a thin source, the following 
reasoning was used. It was assumed that the 
straggling of alpha-particles about their mean 
range R could be represented by a Gaussian 
curve. The rate at which alpha-particles reach 
a distance y from a thin source is then given by 


Ny= Nol 1+@/v) f eax, (2) 


where No is a constant and x =(R—y)/a, a being 
a straggling constant. Let us now suppose that 
there is a uniform distribution of radioactive 
material between the surface of the foil and a 
depth Z, measured in equivalent air distance, 
below the surface. Then the rate at which par- 
ticles leave the layer at depth between z and 
+dz from the surface is 


dNo=(No/Z)ds. 


The rate at which particles starting from: this 
layer reach a point at a distance y from the 
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surface is then 
z—z/a 
an =dNel1+(2/V ~f eters, 
0 


Hence the counting rate at a distance y is 


Z 
N=(N/Z) f dz 


z—z/a 


x 1402/9 of emnatas, (3) 


A first trial range was obtained by extrapolat- 
ing the linear section of a curve drawn through 
the experimental points. This range was found 
to be 3.64 cm. Near the tail of the curve, only 
particles originating near the surface of the foil 
can be counted. Hence in this region, s<]ax|, 
and the curve has very nearly the same shape as 
it would have with a thin source. We may there- 
fore choose a in Eq. (2) to obtain a good fit to 
the tail of the curve. For the data shown, a=0.15 
cm. The dashed curve shows N,, as calculated 
from Eq. (2) with these constants. Extrapolation 
of the tangent of this curve at the point of in- 
flection yields the extrapolated range of (3.81 
+0.03) cm. 

Since Eq. (3) cannot be integrated in closed 
form, the transition from a thin to a thick source 
was accomplished by numerical integration. It 
was found that the summation of 20 thin sources, 


each characterized by a counting rate of 2.3/sec 
and by the dashed curve of Fig. 2, with an air 
equivalent distance of 0.125 cm between suc- 
cessive sources, gives the solid curve of the 
figure, in good agreement with the experimental 
points. 

This apparatus has made possible the attain- 
ment of results of fair precision in undergraduate 
laboratories. The experiment serves as an intro- 
duction to counting techniques and to the study 
of statistical fluctuations. In addition, the analy- 
sis of the data is sufficiently simple to be quite 
clear to undergraduate students, and yet is 
sufficiently complex to indicate that significant 
results are not necessarily revealed directly by ex- 
perimental measurements without some thought 
as to the interpretation of such measurements. 

A number of improvements could undoubtedly 
be made in the apparatus without undue trouble. 
A larger area covered with phosphor would in- 
crease the counting rate without appreciable 
sacrifice of favorable geometry. A still greater 
gain could probably be made if two photo- 
multipliers were used with a simple coincidence 
circuit‘ to reduce the background count. Finally, 
the use of 5819 or similar tubes in place of the 
931-A type would increase the efficiency and 
might well eliminate the necessity of cooling, 
even with weak sources. 


4W. C. Elmore, Rev. Sci. Instr. 22, 649 (1950). 


Summer Meeting of AAPT 


The Summer Meeting of the American Association of Physics Teachers will be held in , 
Pittsburgh, Pennsylvania, on June 25, 26, and 27, 1953. Arrangements have been made with ' 
the Webster Hall Hotel for room reservations for our members. It is important and necessary 
that members write to the hotel individually to receive confirmation of their reservations, men- 


tioning their connection with the Association. 


Most of the meetings will be held in the air-conditioned rooms of the nearby Mellon Institute. 
Titles and abstracts of papers to be contributed to the program should reach the program 
chairman, Professor Marsh W. White, Department of Physics, Pennsylvania State College, 
State College, Pennsylvania, not later than April 30. Several symposiums are planned, and 
active local preparations are being made by the Western Pennsylvania, Central Pennsylvania, 
and Appalachian Sections in collaboration with Professor Chas. Williamson, Carnegie Institute 


of Technology. 





Relativistic Rocket Theory 


W. L. BADE 
University of Nebraska, Lincoln, Nebraska 


(Received August 14, 1952) 


In the relativistic theory of one-dimensional rocket motion in empty space, Ackeret has 
found the differential equation of motion and the mass-ratio formula. One can obtain a formal 
solution for the integrated equation of motion X = X(T) in the general case by determining X 
and T parametrically as functions of the proper time s of the rocket. 


HE subject of the following discussion is a 

rocket! which moves in field-free empty 
space. In case either the velocity of the exhaust 
or that of the rocket itself becomes very large, 
it is necessary to use relativistic in place of 
classical dynamics in order to secure exact re- 
sults. Although such velocities are far out of 
current technical reach, the problem considered 
here is of interest both as a legitimate part of 
the theory of rockets and as an exercise in rela- 
tivity theory. 

This problem may be stated as follows. A 
rocket of specified construction moves in field-free 
empty space, to determine the expected performance. 

The term ‘specified construction” implies 
knowledge at least of 


1. the initial rest mass m, of the rocket; 

2. the exhaust velocity —w, constant rela- 
tive to the rocket; 

3. the rate defined by p(s) =dm/ds at which 
the rest mass of the rocket changes, as 
measured by a clock carried by the 
rocket; s is the proper time of the rocket. 


This rate p(s) may be any sectionally-continuous 
nonpositive function of s such that Jo%p(s)ds 
+m,>0 (m=m, when s=0). The rest mass- 
ratio, given by the equation v=m,/m2= (initial 
rest mass)/(final rest mass), can be computed 
from items 1 and 3 above. The term ‘“‘per- 
formance”’ means 


(a) the terminal velocity U2 of the rocket 
(the initial velocity U; is assumed 
given) ; 

(6) the integrated equation of motion X 
=X(T) of the rocket in an inertial 
system. 


1 That is, any sort of reaction-propelled vehicle. 


Ackeret? has developed the essentials of the 
relativistic theory of rocket motion in empty 
space.’ He has given the fundamental differential 
equation [Eq. (6) below] of such motion, and 
has integrated this equation to obtain a formula, 
Eq. (8), which suffices to determine the final 
velocity U2 in terms of the known mass-ratio 
m,/m:2. However, he has not discussed the prob- 
lem of determining the integrated equation of 
motion. This latter has been solved in the special 
case of ‘‘hyperbolic motion’ by Shepherd.§ A 
general solution is given here. 

To obtain Ackeret’s relativistic differential 
equation of rocket motion, let m be the rest 
mass of the rocket at any point in the motion, 
U the velocity of the rocket in an inertial frame 
K(X YZT), dm’ the rest mass of an element of 
exhaust, U’ the velocity of this element of ex- 
haust in K, and —w the exhaust velocity relative 
to the rocket, a constant. 

The rocket is assumed to move along a recti- 
linear path in the inertial frame K; let the X axis 
of K lie along this path. Assume that all the 
elements of exhaust are ejected straight back 
along the X axis, so that the internal efficiency 
of the rocket drive is 100 percent and the problem 
is essentially one-dimensional.® 

On account of the assumption of “‘field-free 
empty space,” the rocket, together with its ex- 
haust, constitutes an isolated system. Hence the 
total four-momentum of this system is constant; 


2 J. Ackeret, Helv. Phys. Acta 19, 103 (1946). 

3For a readily available review of Ackeret’s work in 
English, see Seifert, Mills, and Summerfield, Am. J. Phys. 
15, 255 (1947), pp. 267-269. 

4 That is, motion in which the acceleration of the rocket 
as measured in the instantaneous proper system is constant. 
See C. Mauguin, Compt. rend. 234, 1004 (1952). 

5 L. R. Shepherd, J. Brit. Interplan. Soc. 11, 149 (1952), 
Appendix. 

6 Shepherd has investigated a case where this assump- 
tion is not satisfied. See L. R. Shepherd, reference 5. 
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that is, 
A & mdx*/ds=0. (1) 


The summation is over the particles of the sys- 
tem, m is the rest mass of a sample particle and 
dx*/ds a component of its four-velocity, while A 
denotes taking the difference of values before 
and after an interaction among the particles. 
The relativistic formula for addition of velocities 
gives in this case 


U-—w 
~ 1—(Uw/e) 


, 


(2) 


When Eq. (1) is applied to the expulsion of an 
element of exhaust from the rocket, there results 


d[{mdX /ds|+ (dX /ds)'dm'=0, (3) 
d{mdT /ds]+(dT/ds)'’dm' =0, 


in which (dX/ds)’ and (dT/ds)’ are components 
of the four-velocity of the ejected material. If 
now dm’ is eliminated from Eq. (3), there results 


d[ mdX /ds]— U'd[mdT/ds]=0, 


where 


U’ = (dX /dT)' = (dX /ds)'/(dT/ds)’. 
Let 
ds? =dT?—c?*(dX?+-d Y?+dZ?’). (4) 


Then, since dY=dZ=0 here, 
dT/ds=(1—(U/c)?}4, 
dX /ds= U(1—(U/c)*}-3*. (5) 


Making use of Eq. (2) and Eq. (5), one obtains, 
after algebraic reduction, Ackeret’s equation 


wl1—(U/c)? |}dm+mdU=0. (6) 


One readily finds as a first integral the mass- 
ratio formula 


c+ U2 c/2w: c— U, c/2w 
Ey « 
Me c— U2 c+ U, 
in which U, is the initial, U2 the final velocity 


of the rocket in K. When J, is set equal to zero 
and Eq. (7) is solved for U2, one finds, with 


Ackeret, 
prwic_ ff 
vid | (8) 
prwlet 4 
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Formulas (7) and (8), determine part (a) of the 
“performance,” since w and v=m,/m, are 
known. 

When it comes to determining the second part 
of the “‘performance,”’ i.e., the integrated equa- 
tion of motion, a slight difficulty appears. If the 
function P(T)=dm/dT were given (T=time 
read by a clock at rest in K), it would be possible 
to integrate Eq. (6) directly for an equation 
of motion. However, the rate at which rest 
mass is lost by the rocket is determined by 
mechanisms or agencies on board the rocket. 
Hence, not P(T)=dm/dT but p(s)=dm/ds is 
specified by the construction of the rocket. Be- 
cause of the Einstein time dilation, these two 
rates are not the same. 

One can overcome this difficulty and obtain 
the integrated equation of motion (in para- 
metric form) as follows. From Eq. (4), with 
dY=dZ=0, it is clear that a function 6(s) 
exists, such that 


cosh@=dT/ds, csinhé@=dX/ds. (9) 


Then 
U=dX/dT =c tanhé. (10) 
Substituting Eq. (10) into Eq. (6), one obtains 
d6/ds = —(w/mc)p(s). 


As the quantities on the right are specified by 
the construction of the rocket, it is possible to 
integrate and obtain 


6(s) = =f as) 5. 


c ™m 


(11) 


m=m,+ f p(s)ds. 
0 


Then, from Eq. (9), 


r= f coshé ds, x=cf sinhé ds. (12) 
0 0 


Equations (12) determine X and T parametri- 
cally as functions of the proper time s of the 
rocket. It is assumed that the system K has 
been so chosen, and the clocks in K and on the 
rocket have been so set, that when s=0, X 
=T=U=0. 

It is of some interest to note another way of 
obtaining this result, Eq. (12). The preceding 
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calculation was carried out entirely in the in- 
ertial system of reference K(X YZT). One might 
instead choose to employ the system k(xyzt) in 
which the rocket is always at rest at the origin. 
As this system is not an inertial one, it is then 
necessary to borrow techniques from the general 
theory of relativity. 

The problem of accelerated frames of reference 
has been discussed by Mller.? He has shown 
that the transformation equations from k to K 
are 


t 
xX =~ cosho-+< f sinhé dt 
: (13) 


t 
T= (x/c)sinho+ f coshé dt, 
0 


and that the metric in & is 


ds*= (1+ gxc~*)*dt? —c~*(dx*+dy*+dz"), (14) 


7C. Mgller, Kgl. Danske Videnskab. Selskab, Mat.-fys. 
Medd., Bind XX, Nr. 19 (1943). 


6(t) = f ‘edt, 


0 


and g(t) is the acceleration of the origin of k as 
measured in the instantaneous proper system. 
From Eq. (13) it follows that the velocity in K 
of a point fixed in k is 


U=c tanhé. (15) 


According to Eq. (14), for a point fixed at the 
origin of k the coordinate time ¢ is equal to the 
proper time s. Since the rocket resides perma- 
nently at x=y=z=0, it follows that Eq. (15) 
is identical with Eq. (10), 0(¢) =6(s) is again de- 
termined by Eq. (11), and Eqs. (13) reduce to 
Eqs. (12). 

According to the viewpoint of the system of 
reference k, the rocket is situated in a gravita- 
tional field of intensity g(t) and keeps itself sta- 
tionary at the origin by use of its reaction engine. 


“since bodies, having but one common matter, can be differenced but by accidents, which seem all 
of them to be the effects and consequents of local motion, I see not why it should be absurd to think, 
that (at least among inanimate bodies) by the intervention of some very small addition or subtrac- 
tion of matter, (which yet in most cases will scarce be needed,) and of an orderly series of altera- 
tions, disposing by degrees the matter to be transmuted, almost of any thing, may at length be made 
any thing: as, though out of a wedge of gold one cannot immediately make a ring, yet by either wire- 
drawing that wedge by degrees, or by melting it, and casting a little of it into a mould, that thing 
may be easily effected. And so though water cannot immediately be transmuted into oil and much 
less into fire; yet if you nourish certain plants with water alone, (as I have done,) till they have 
assimilated a great quantity of water into their own nature, you may, by committing this trans- 
muted water . . . to distillation in convenient glasses, obtain, besides other things, a true oil, and 
a black combustible coal, (and consequently fire;) both of which may be so copious, as to leave no 
just cause to suspect, that they could be any thing near afforded by any little spirituous parts, which 
may be presumed to have been communicated by that part of the vegetable, that is first put into the 
water, to that far greater part of it, which was committed to distillation” .—Quoted by T. S. KuHn 
from The Works of the Honourable Robert Boyle, Volume 2, Page 375 (London, A. Millar, 1744). 
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Effect of Absorption by the Material of the Prism 
on Its Resolving Power 


MAHENDRA SINGH SODHA 
Allahabad University, Allahabad, India 


PARROW* has discussed the effect of absorption by the 
material of the prism on its resolving power on the 
basis of his criterion. As his criterion has been abandoned 
in favor of the Rayleigh criterion, expressed in terms of 
intensity as 


Tmin/ Imax - 8/r’, 


the author proposes to discuss the problem on the Rayleigh 
criterion. 

Considering the incidence of a plane wave front on the 
prism, Sparrow’s expression for the intensity of the emer- 
gent wave front, in a direction inclined at an angle 6 to 
that of the central maximum, can be put as 


I= (1—m)?+4m sin*’a 
e (logm)? +402 


where m=e-*’/2, ¢’ being the longest path traveled by any 
ray in the prism, & the intensity absorption coefficient and 
a=na sin6/\,a being the width of the emergent wave front. 

If the central maxima of two spectral lines (given by 
a=0) have an angular separation 5@ corresponding to Aa, 
the maximum and minimum of the resultant intensity 
pattern are given by 


(1) 


Tmax = Laut t-lamdta 
. (1—m)? . (1—m)?+4m sin?(Aq) 


. (logm)? (logm)?+4(Aa)? 


(1—m)?+4m sin®?(Ae/2)? 
(logm)?+ (Aa)? 

Rayleigh’s criterion for resolution requires 

8 (1—m)? (1—m)?+4m sin*Aa 

x (logm)?+4(Aa)? } 

(1—m)?+4m sin?(Aa/2) 
(logm)*+ (Aa)? 
This equation gives one negative value of (Aa). Taking the 

square root of the positive value we have, in general, 
Aa=2/f(m), (3) 
where f(m) is a function of m. We have 
a/f(m) = Aa=xa/d{sin(6+656) —sind}. 
= 7a/d(0+50— 0) = 1a50/X, 


since @ is small. This result may be expressed in terms of 
resolving power R by the equation 


R=)/dy\=f(m)ad6/dd = Rof (m). (4) 
From Eq. (2) we find that the maximum resolving power 
occurs when k=0, m=1, and f(m) =1, i.e., when there is 


no absorption. We also find that the resolving power 
decreases slowly with decreasing m. 


I min = 2Ta=Aa/2= 2Io 


n° (logm)? 


=21 >, (2) 
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An approximate method for the calculation of f(m) for 
values of m that are of practical importance (0.8 or more) 
has been developed by the author. 


When m2 0.8, (logm)? is very small compared to (Aa)? 
and may be neglected. Moreover, since Aa is only slightly 
different from + we can put sinda=0 and sinda/2=1. 
Introducing these approximations in Eq. (2) we get 

x (1-+m)? (1—m)?  (1—m)* 
4 (Aa)? (logm)? 4(Aa)? 


w(1+m)?—(1—m)? (1—m)? 
~— 4@a)t— logm)? 
Putting Aa=2/f(m) we get 
2(1—m) T 
= 2(1—m) 1 
2(1—m) (1—m)? 
. eee a nee). 6) 


where |logm| denotes the numerical value of logm. 
The author thanks Dr. K. Majumdar and Mr. Y. P. 
Varshni for their interest in the investigation. 


~ (logm)? 


f(m) = 


1C. M. Sparrow, Astrophys. J. 44, 76 (1916). 


Determination of the Coefficient of Surface 
Tension by the Bubble-Length Method* 


Joun C. SLonczEwskit AND RALPH HELLER 
Worcester Polytechnic Institute, Worcester 2, Massachusetts 


HEN a comparatively large air bubble is entrapped 

under a horizontal plate within the test liquid, an 

approximate value of the coefficient of surface tension, y, 

of the liquid may be obtained by measuring its density p 

and the distance h between the horizontal planes through 

the vertex of the bubble and its widest (equatorial) section, 

for then! 

2y/pg=h? (1) 

nearly. The left-hand member of Eq. (1) is termed the 
capillary constant of the liquid. 

In an interesting modification of this method of “‘sessile 
bubbles,’’ Euverard and Hurley? measured the over-all 
length / of an air bubble of given volume V in a horizontal 
tube of known internal diameter 27, filled with the test 
liquid. It is true that the mathematical relation which 
gives y in terms of p, 1, V, 2r (and possibly the angle of 
contact) is not known, but Euverard and Hurley were able 
to show that a smooth curve, approximated satisfactorily 
by one branch of an equilateral hyperbola, was obtained 
when the value of /p~! for each of 18 liquids was plotted 
against y (which, like p, was taken from tables). Stated 
mathematically, 

(ya) (lp-!—b) =c, (2) 
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Fic. 1. The reciprocal of the capillary constant, i.e., pg/2y, as a 
function of the over-all length, 1, of the bubble. (27 =1.283 cm; V =1.00 
cm!.) The experimental points represent the following liquids, from left 
to right: Water, glycerol, aniline, benzyl alcohol, nitrobenzene, benzene, 


methy! alcohol, -hexane, ethyl ether, D-C 200 silicon, carbon tetra- 
chloride, ethyl bromide. 


where a, b, and ¢ are constants determined experimentally 
for given V and 2r. With V and 2r fixed, it was, therefore, 
possible to make a fairly accurate estimate of the value 
of the coefficient y of a test substance by merely measuring 
the horizontal length, 1, of the air bubble with vernier 
calipers, which can be done far more easily and more 
quickly than the measurement of a small* vertical dimen- 
sion such as h. 

We have repeated the experiment of Euverard and 
Hurley, using 12 liquids (mostly organic and chemically 
pure) at room temperature, and our observations fit a 
hyperbolic curve of the form of Eq. (2) fairly well. It was 
possible to hold the volume, V, of the bubble (of the order 
of 1 cm?) to within 4 percent of a given value for all liquids 
by employing the following procedure. A graduated pipette 
was connected to one end of the tube, the tube was filled 
completely with the liquid, a short capillary tube was 
attached to the other end, and the desired volume of air 
was then admitted through the capillary tube by drawing 
an equal volume of liquid into the pipette. The length J of 
the bubble was usually measured by us with the aid of a 
small micrometer microscope. 

As to the mathematical relation among the various 
quantities involved, the variables in Eq. (2), with V and 
2r held constant, are y, p, and J. Dimensional reasoning 
suggests that the ‘‘capillary constant,” 2y/pg, should be 
a function of the measured length / (with coefficients 
depending on V and 2r), for 27/pg is of the dimension of 
the square of a length (as may be seen from Eq. (1)). 


Hence 2y/pg=f(l). When the reciprocal, pg/2y, of the 
capillary constant was plotted as a function of the length, 
1, of the bubble, a practically straight line of the form 


pg/2y=Al+B (3) 


was obtained. Our data, secured with a polarimeter tube of 
1.283 cm internal diameter and a bubble of 1.00 cm? 
volume, are shown in Fig. 1. The straight line found has 
the slope A =5.0cm™ and the intercept B= —6.1 cm=*. 
The data of Euverard and: Hurley likewise fall satis- 
factorily on a straight line of the form Eq. (3) when so 
plotted. Using their observations (with 2r=0.8cm and 
V=0.5 cm’) we find A=9.5cm-* and B=—10.4 cm. 
Once this straight line has been determined for a tube of 
given diameter and an air bubble of definite volume, the 
coefficient y of a test liquid (of known density) can, there- 
fore, be estimated with fair accuracy from the measured 
length of the bubble. 

Additional investigation of the bubble-length method, 
carried out both theoretically and experimentally, should 
prove worth while. 

* A portion of the thesis submitted to the Department of Physics, 
Worcester Polytechnic Institute, by John C. Slonczewski in partial 


fulfillment of the requirements for the degree of Bachelor of Science 
(June, 1950). 

+ Now National Science Foundation Fellow, Department of Physics, 
Rutgers University, New Brunswick, New Jersey. 

1 See, for instance, B. L. Worsnop and H. T. Flint, Advanced Practical 
Physics for Students (Methuen and Company, London, 1951), ninth 
edition, pp. 119-122. The experimental method is due to Georg Quincke, 
Pogg. Ann. 139, 7-8, 11-14 (1870). 

2M. R. Euverard and D. R. Hurley, Anal. Chem. 21, 1177-1180 
(1949). We wish to thank Professor K. L. Mayer and Professor E. D. 
Wilson for bringing this work to our attention. 

3 At room temperature, h as given by Eq. (1) ranges from 0.18 cm for 
ethyl bromide to 0.39 cm for water. 


Equivalent Networks of Electron Tubes 


SauL ARoNow 
Harvard University, Cambridge, Massachusetts 


HE four-terminal equivalent network of an electron 

tube is well known. However, the student is often 
confused by the variety and arbitrariness of the sign con- 
ventions used by different authors. Clarity in understand- 
ing equivalent networks was not very important when 
only the simplest methods of tube operation were in vogue. 
But now that the unconventional is ordinary, with new 
types of “‘tubes,’’ such as transistors, etc., and new tech- 
niques of circuitry, it becomes imperative that uniform, 
easily remembered conventions be adopted. 

The reason for the confusion is a conflict of intuition and 
mathematical necessity. Consider an active four-terminal 
network, arbitrarily assuming directions for the terminal 
currents and voltages as in Fig. 1. There is a general func- 
tional relation 


ip=f(Up, U9), (1) 


Fic. 1. Four-terminal “‘black box” for triode. 
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which leads to an identity in the partial derivatives, 


GE = (2) 


In setting up the small signal, linear, equivalent circuit for 
the system, we would like to consider the three derivatives 
as physically significant coefficients, and as such to have 
positive numerical values. But the mathematics requires 
that at least one of them be a negative number. Early 
circuit designers felt that it was “natural” for an amplifica- 
tion to be a positive quantity. Hence, they defined the 
amplification factor as the negative of the derivative, 


aie ~~ 
Hop= av 


since the static characteristic, in as negative grid region, 
has a negative slope. This convention leads to the common 


equivalent circuit having a distant-controlled generator 
with inverted polarity. 


If the tube draws grid current, we may set up an identity 
for the grid similar to Eq. (2): 


(=) (=) _ —(% 
Bigs vp\dvp/v, dvi, 


Now defining 


dup 
80,/ ip 


(3) 


= +{— d = +(=* 
= (a a =). 


Mpg appears, from the slope, as a positive and gy, as 
a negative number. Again, if the tube is operated in the 
dynatron region of its static characteristics,’ the slopes of 
the characteristics change and all three coefficients may 


be negative. To redefine the signs of the coefficients as the 
operating region changes would be hopelessly involved. 

Apparently in recognition of these difficulties, the Insti- 
tute of Radio Engineers Standards** have been revised to 
establish a logical, uniform procedure. The Standards 
define each coefficient as the positive partial derivative, 
i.e., with numerical value equal in magnitude and sign to 
the slope of the appropriate characteristic curve. The 
Standards specify the polarities of the equivalent network 
in terms of current generators as in Fig. 2(a), corresponding 
to the linear equations 


Th=yuVityi2Ve, 
I2=yoiVitye2 V2. 
This may be converted to voltage generators as in Fig. 2(b), 
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Fic. 2. Linear equivalent networks for electron tubes. 
(a) Current generators. (b) Voltage generators. 


corresponding to 
Vi=2iliturVe, 
Vo=uaVitzeele. , 


It should be noted that the coefficients are defined in a 
very flexible manner, with currents and voltages measured 
at whatever points it is convenient to designate as external 
terminals, so that with proper interpretation, the network 
may be used for many devices. The importance of a 
systematic procedure becomes more evident when it is 
noted that while the Standards, as always, have fallen 
behind recent technical developments, the basic principles 
still apply. The networks in Fig. 2 exhibit only 8 of the 24 
possible partial derivatives. In transistor circuits, for 
example, it is more convenient to define current amplifica- 
tion factors or to use systems where one current and one 
voltage are considered dependent variables. However, if 
the basic definitions and sign conventions of the Standards 
are followed, any arbitrary selection of variables may be 


made and the resulting equivalent network written down 
without hesitation. 


L. Chaffee, Theory of Thermionic Vacuum Tubes (McGraw-Hill 
Book Poon Inc., New York, 1933), p. 153. 
2 IRE Standards on Electron tubes; — of Terms, 50 IRE 
7.S1, Proc. Inst. Radio Engrs. 38, 426 (1950). 
3 IRE Standards on Electron Tubes, 50 IRE 7.S2, Proc. Inst. Radio 
Engrs. 38, 935 (1950). 
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An Extension of the Conical Pendulum Problem 
and Its Demonstration 


T is a simple exercise to show that the height & of the 
conical surface swept out by the string and bob in a 
conical pendulum motion is given by g/w*. That is, h 
depends only on the speed of rotation and not upon the 


weight of the bob or the length of the supporting string. If, 
then, a number of different masses, suspended from the same 
point by strings of different lengths, are given the same 
angular velocity about a vertical axis through the point of 
suspension, should not the masses all lie in the same hori- 
zontal plane? This is a necessary deduction, but not intui- 
tively clear. Indeed, this is exactly the kind of situation I 
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have referred to earlier where, as I said, I believe the result 
but would like to see it. 

The demonstration is easily accomplished by suspending 
bobs of different masses on strings of different lengths 
from a shaft held in the chuck of a rotator. A short piece of 
rod or metal tubing drilled along several diameters makes 
a good shaft. The method of fixing the string to the shaft 
is critical but not difficult to adjust. With different colored 
bobs the demonstration has a measure of beauty and 
lends itself nicely to shadow projection. 


Jutius SUMNER MILLER 
Ford Foundation Fellow 


10303 Mississippi Avenue 
West Los Angeles 25, California 


The Ultrasonic Interferometer as a Screw 
Calibrating Instrument 


RECISION screw calibration requires a small standard 

of length, precisely known, which can be easily 
compared with the linear advance of the screw being 
calibrated. The best standard of length is the wavelength 
of light, so that the optical interferometer may be con- 
sidered the ultimate screw calibrator despite the fact that 
the tediousness of the screw-to-standard comparison makes 
its use rather impracticable. But recent improvements in 
the acoustic interferometer, which uses the wavelength of 
sound as the standard of length, suggest that this instru- 
ment might be more successful in calibration applications 
where the ultimate in accuracy is not required. 

The principal advantage of the acoustic interferometer 
as compared with its optical counterpart is its tremendous 
superiority in speed and ease of the screw-to-standard 
comparison. The ultrasonic interferometer developed by 
Barthel and Nolle! for measurement of dispersion of 
ultrasonic sound in liquids is capable, with some modi- 
fication, of turning out a permanent record of a complete 
screw calibration within a minute. Preliminary experiments 
show that the instrument is capable of recording cumula- 
tive linear increments of 0.003+-0.00004 cm. For a screw 
with a one millimeter pitch, for example, the size of the 
increment thus corresponds to a record point about every 
11° of rotation of the screw. It seems likely that develop- 
mental research could still further reduce the size of the 
increment and increase the precision. 

The instrument referred to above is an interferometer 
of the two-crystal type in which the separation of the 
crystals, which are immersed in a liquid medium, is 
continuously varied. The recording process is equivalent 
to placing a mark on a rotating drum whenever the 
receiving crystal has moved a distance of one wavelength 
of sound in the liquid. This high-precision instrument 
utilizes low intensity traveling waves, rather than the high 
intensity standing waves used in most other interferom- 
eters, and thereby avoids the errors arising from local 
heating of the liquid. It surpasses other currently used 
interferometers in its rapidity of operation as an automati- 
cally recording instrument. For further details reference 
should be made to the original article describing the 
instrument. Much of the mechanical complexity of the 
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original instrument would be unnecessary in an instrument 
specifically designed for screw calibration. 


ROMARD BARTHEL 
St. Edward's University 
Austin, Texas 


1R. Barthel and A. W. Nolle, J. Acoust. Soc. Am. 24, 8 (1952). 


Concerning the Electric Charge on a 
Moving Vehicle 


T is common knowledge that moving vehicles acquire an 

electrostatic charge, and the mechanism appears to be 
frictional. Gasoline trucks are required by law to provide 
.a ground, which may be a dragging chain or strap, and 
some passenger cars are so equipped. At a ferry dock near 
here the approaching roadway carries an upright flexible 
steel strap which is intended to carry to ground any charge 
which an approaching vehicle may bear. The question as 
to how great a charge may accumulate is difficult to 
answer; it is not easily computed theoretically and practical 
measurements offer trouble. This note is intended to 
report an observation which point$ up further difficulties. 

On a dry summer day, not far from desert terrain, while 
driving at a distinctly slow rate, a friend overtook me in 
his car and came abreast of me on my left. His wife in the 
front seat reached out and handed me a package of gum. 
I reached for it with my left arm extended (a dangerous 
business, admittedly!), but came no closer than three 
inches, perhaps, if indeed as close as that. At this instant 
a terrific discharge took place which possessed the classical 
physiological effects. The shock was momentarily disabling, 
as a three-inch spark in air can well be. Both of us 
recovered promptly but I have since wondered what the 
mechanism of storage of charge might be. Simple frictional 
effects can produce a substantial charge but leakage from 
the many sharp edges of an automobile body must in- 
evitably arise, unless some capacitance effect exists. What 


can we say about this? 


JuLius SUMNER MILLER 
Ford Foundation Fellow 
University of California 
Los Angeles, California 


Demonstration of Subjective Harmonics 


ee harmonics produced by nonlinear re- 
sponse of the human ear have been known for a long 
time; yet they are not easily demonstrated under un- 
impeachable conditions. What is needed is a high intensity 
sinusoidal stimulus; and that is not likely to be furnished 
by any ordinary combination of oscillator, amplifier, 
and loud-speaker or headphones. 

Some years ago there was described in this Journal a 
thermionic tube circuit with which tuning forks can be 
driven to large and controllable amplitudes.! Tests with a 
high quality microphone and associated equipment do 
not show significant distortion in the air-borne sound. 
Yet, if an observer closely approaches one of our large 
128-cycle forks when it is thus driven, he hears what is 
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known in music as a common chord, with frequencies 
easily identified as 512, 640, and 768 cycles/sec. (With a 
little attention a few of the lower and higher harmonics 
may also be heard.) This effect promptly disappears when 
the observer moves away from the fork. 

The effect is so conspicuous that the observer notices it 
immediately; he does not need to be told in advance what 
to listen for. It cannot be due to overtones produced by 
the fork itself, for these are nonharmonic. The frequencies 
of the first two overtones of this fork have been measured? 
and were found to be 807 and 2154 cycles/sec. 

It seems reasonable to conclude that true subjective 
harmonics are being observed. 

CHARLES WILLIAMSON 


Carnegie Institute of Technology 
Pittsburgh, Pennsylvania 


1C. Williamson and W. J. Eisenbeis, Am. J. Phys. (Am. Physics 
Teacher) 4, 91 (1936). 

2 This measurement is conveniently and precisely made with our 
apparatus by exciting the driving amplifier with the output of an 
audiofrequency oscillator and checking the frequencies of the audible 
response peaks with a Stroboconn. 


Trigonometry before Physics? 


HYSICISTS today are offered many and varied 

opportunities in both the theoretical and applied 
fields, but strange as this may seem, there are insufficient 
qualified personnel to fill these positions. Therefore the 
question arises ‘‘why?” A possible answer may be found 
in the high mortality rate among physics majors in the 
colleges. Again the question, ‘“‘why?”’ A possible reason is 
that on application to college many find they have in- 
adequate background for physics. For example, today 
some are applying with little or no knowledge of trigonom- 
etry. The college offering the student a choice of taking 
general physics either in the freshman or sophomore year 
gives him an opportunity to remedy this deficiency. But 
to the college offering general physics in the freshman year 
only, it presents a serious problem. 

In an attempt to ascertain the attitude of colleges 
throughout the country and thus seek a solution to this 
problem, a questionnaire was circulated in May 1952 to 
145 institutions of higher learning that were representative 
of all important types found in the United States. Its 
primary purpose was to learn the attitude of other schools 
toward the question, ‘‘Do you regard trigonometry as a 
requisite or prerequisite for general college physics?” 
(Incidentally, a question involving trigonometry for the 
biology major was also inserted.) Most colleges solve that 
problem by giving general physics in the sophomore year, 
and those giving it in the freshman year either teach 
trigonometry in course or run a class concurrently. 

The concern, however, is with the physics major and 
his ability to survive his course. With inadequate back- 
ground he is always behind, for while new material is 
presented he does not have the necessary tools to handle 
it and, therefore, a high mortality rate exists. This is true 
not only in his physics courses; he also finds himself 
failing in his mathematics courses. Where, then, does the 
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solution to this problem lie? Judging by the requisites as 
indicated in many college catalogs it seems that colleges 
themselves are partly responsible for this condition, for 
it is difficult, if not impossible, to determine just what 
subjects are required of the high school graduate. This 
vagueness could very well lead a high school principal to 
say, as he did, “You won’t have to bother with trig, 
because you won't have much need of it.” 

The following are the results of the questionnaire. One 
hundred and forty-five were sent throughout the United 
States and one hundred and twenty-nine replies were 
received. 


i. Number requiring trigonometry as an entrance 
requirement for physics majors 
. Number requiring trigonometry as an entrance 
requirement for biology majors 
. If the entrant has not had trigonometry do you 
permit him to start the General Physics course? 
No 97 
Yes 5 
Pre-requisite or concurrently 17 
Impossible to determine 3 
Recommended $ 
Contradictory 4 
. Number not requiring trigonometry as a pre- 
requisite for physics majors 
. Number requiring trigonometry as a prerequi- 
site for physics majors 
Yes 
No 
Recommended 
Taught in Course 
Impossible to determine 


113 


Note that five colleges (two universities, one Eastern 
college and two Midwestern state colleges) do not consider 
trigonometry as a prerequisite, but their remarks indicate 
otherwise. One reported, ‘‘Some introduction to trigonome- 
try is given in the first semester of the first year in physics. 
In the second semester familiarity with elementary trigo- 
nometry is assumed, though no course is required as 
prerequisite.”” One wonders how the familiarity is obtained. 

Another reported, ‘‘Physics majors must take calculus 
before junior year. Doubtful if any student who has not 
had math through trig has ever tried to major in physics. 
It would be very difficult, but not impossible.” (!!) 

The strangest reason of all was given by a college from 
the Midwest, ‘‘We do not consider it the privilege of the 
university to dictate to the secondary school that they 
require or even offer trig.” 

Still another wrote, ‘‘This college has no high school 
math requirement. All students, however, are required to 
take a precollege math test.” 

The above are direct quotes from representative colleges 
in different sections of the country. 

On the positive side of the question, not all are in 
agreement. One wrote that he could “teach sufficient trig 
in two or two and one-half hours.” Another stated ‘We 
think we can teach sufficient trig in course for the student 
to get by.’’ Still another, “Trig is not an entrance require- 
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ment in any curriculum. All students in technical curricula 
are required to spend 75 hours (5 credits) in algebra and 
trigonometry.” 

But on the definitely positive side, a Western college 
wrote, ‘Mathematics instruction in the public schools of 
this state has reached the point where relatively few 
entering students have had trigonometry. We have, there- 
fore, been compelled to introduce this subject into a 
unified mathematics course for the freshman year which is 
required of all students taking physics.” 

Akin to this problem, Dr. Geer of the University of 
Southern California remarked! that “only 18 students out 
of 1800 complete trigonometry in a high I.Q. high school 
in Los Angeles.”” He points out the shortage of scientists 
and engineers and states that one of the contributing 
factors is ‘‘the omission of this course in many of our high 
schools, even those preparing a boy to pursue a science 
course in college.’’ The present survey entirely agrees with 
the above. 

The following points are offered as possible solutions to 
the problem: 
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Modern Physics. F. W. Van NAME, JR. Pp. 360+-vii. 
Prentice-Hall, Inc., New York, 1952. Price $5.50. 


The teaching of an undergraduate modern physics 
course usually presents something of a problem with regard 
to organization of material. Unlike mechanics or electricity 
and magnetism courses a logical, systematic approach is 
almost impossible, resulting both from the catchall nature 
of such courses in which everything except classical physics 
is included, and from the relatively sophisticated viewpoint 
necessary to unify the aspects of physical phenomena 
usually treated. Hence the available texts fall naturally 
into two classes, depending on whether the authors are 
primarily interested in presenting a great deal of miscel- 
laneous information as palatably as possible or in providing 
a sound and adequate prelude to further work; i.e., they are 
either topical discussions or introductions to quantum 
theory. There is much to be said for each of these ap- 
proaches, but in general the former is more useful for the 
less advanced student requiring a broad although not 
necessarily detailed knowledge of modern physics. 

Modern Physics attempts to fulfill this latter function 
by presenting a well-written survey of the field intended 
for a full year course. The book is divided into three 
sections, devoted to The Electron, Atoms and Molecules, 
and Nuclei and Nuclear Particles, with each of the chapters 
covering a specific topic. Chapters dealing with kinetic 
theory, low temperature phenomena such as superconduc- 
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. Greater unanimity among the professors in the 
field both on the college level and on the secondary 
level. 

. Colleges should state definitely, clearly and con- 
cisely their general entrance requirements and 
special requirements for scientific courses. Special 
provision should be made for the unusually capable 
applicant who does not have the necessary back- 
ground. 

3. A definite effort to integrate the high school curric- 
ula with the specific required college curricula. 


The deans of studies in the secondary schools would then 
know definitely the collegiate requirements and plan the 
secondary curriculum accordingly. The high school 
graduate would then be adequately prepared to embark 
upon his chosen major in college. Hence, he would not be 
continually ‘‘one year behind,’’ nor would he be forced to 
withdraw for what is weak background, not low mentality. 


W. A. MurtauGu, O. P. 
Providence College 


Providence 8, Rhode Island 
tAm. J. Phys. 20, 397-398 (1952). 
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tivity and the behavior of the He II, cosmic rays and 
mesons, and quantum mechanics supplement the orthodox 
ones on atomic spectra, x-rays, special relativity, etc. 
A discussion of nuclear fission and its applications is 
included. The author makes a distinct effort to encourage 
understanding by providing complete derivations of the 
more important equations and by the inclusion of worked- 
out numerical calculations, in many cases. However, there 
is a lack of uniformity about the levels of the various 
chapters, with the one on low temperature physics, for 
example, being almost entirely qualitative while the quan- 
tum mechanics part is rather a difficult morsel for the 
student to swallow at the stage in his education when he 
would be using a book of this type. It is likely that most 
instructors would prefer a more elementary exposition, 
although there is no reason why this material could not 
be covered adequately if a healthy amount of time were 
devoted to it. 

At first glance this book seems to consist of a larger 
ratio of words to formulas than one encounters elsewhere 
in literature of this kind, but it is necessary to consider 
that the problems appended to each chapter are an 
integral part of the text and their solution is necessary for 
full understanding of the material. Keeping this in mind, 
Modern Physics is welcome addition to the many works 
already in use, some of which it is certain to replace. 

ARTHUR BEISER 
New York University 
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Physics and Medicine of the Upper Atmosphere. Edited 
by Crayton S. WHITE AND Otis O. BENSON, JR. 
Pp. 611+xxiv, Figs. 150, plates 46. University of 
New Mexico Press, Alburquerque, New Mexico, 1952. 
Price $10.00. 


Rather than a treatise or a text in the usual sense, the 
book presents an edited group of reports giving the 
proceedings of a symposium on the field indicated by its 
title, held in San Antonio, Texas, on November 6, 7, 8, 
and 9, 1951, and sponsored by The Air University School 
of Aviation Medicine. The arrangements were made by 
The Lovelace Foundation for Medical Education and 
Research of Albuquerque, New Mexico. 

Many of the thirty-eight contributors are well known 
and outstanding physicists or medical men, and the list 
includes high ranking officers of the United States Air 
Force, which officially supported the project. The subject 
matter deals primarily with the physics of the extreme 
upper atmosphere with emphasis on those aspects which 
present hazards to aviators who now, or eventyally, may 
operate in this region. The medical contributors to this 
text have attempted to appraise the seriousness of these 
hazards from the standpoint of the aviator. The writers do 
not question the ability of the aeronautical engineer to 
produce craft capable of reaching and traveling in the 
aeropause, as the border region between outer space and 
the atmosphere of the earth is called. They take it for 
granted that the limitations to such travel are human, 
not mechanical. Certainly the current developments of 
rockets and of rocket-type craft of various designs support 
this point of view. Actually much of the material in the 
book is based on scientific data secured recently by 
instruments carried aloft by rocket propulsion. Many 
important but only lately discovered facts are quite at 
variance with the notions long held and taught in physics 
classes, and physicists generally would do well to read 
carefully this new book before using again any old “lecture 
notes”’ dealing with the upper atmosphere. This is partic- 
ularly true with respect to air temperatures, air composi- 
tion, and to the nature and the intensity of the various 
radiations of importance in the aeropause. The present 
tidal wave of interest in nuclear physics finds added 
strength in the new observations made in high flying 
planes and in results obtained through the use of rockets. 
But the medical contributors are genuinely concerned as to 
the possible effect of particle radiations, as well of meteor- 
oids and the like, on any human venturing into the space- 
atmosphere border regions, or even into the highest layers 
in the atmosphere now reached by manned, rocket-driven 
planes. The enthusiasm of youth as well as of many 
uninformed adults for “space travel” would be cooled by 
reading in this book scientifically reliable reports on the 
many difficult problems which must be solved before a 
human may try to reach the altitudes now attained by 
unmanned rockets. 

The subjects stressed in the more than forty papers 
included deal largely with the problems of temperature, 
pressure, and composition of the air, and electromagnetic 
and particle radiations. Only a few chapters are heavily 


319 


laden with theory. Both the physicists and the medical 
scientists may find indicated new types of experimental 
techniques with which they have had no experience in the 
laboratory. The question as to whether human beings can 
tolerate the cosmic radiation found to exist at extreme 
altitudes is discussed by several of the contributors, each 
expressing his own opinions on the subject. Whether 
fortunate or otherwise, the contributions on any subject 
may not fully agree, but this might well be expected on a 
subject of which the element of speculation is so largely 
involved. This does suggest what some may criticize—that 
there is a great deal of overlapping, and disproportionate 
space may be given to subjects of minor importance. But 
the reviewer considers this inevitable in a book of collected 
papers unless the editors should choose to completely 
revamp much of the material. To have done this would 
have not only dimmed the individualities of the contribu- 
tors but greatly delayed the publication of this up-to-the- 
minute material. The authors, editors, and the publishers 
are all to be congratulated on having published such a 
volume in less than a year from the time the material 
became available. The book presents not only a great deal 
of new material directly but includes at the ends of the 
chapters extended bibliographies. These should be ex- 
tremely valuable to research workers and to all others 
concerned with the problems discussed. A most interesting 
feature of the book is the large collection of plates included 
at the end of the text proper. Without doubt this book 
will find a place in the libraries of physics departments and 
of medical schools generally. It is well printed. 
E. L. HARRINGTON 
Riverside, California 


Textbook on Sound. J. W. WINSTANLEY. Pp. 239, Figs. 154, 
14X22 cm. Longmans, Green and Company, London, 
1952. Price $2.60. 


This book is designed, in the English market, for what 
amounts to high school senior or college freshman level. 
At this level the author treats the usual gamut of acoustical 
topics: wave motion, audible acoustic phenomena, ultra- 
sonics, musical and architectural acoustics, instrumenta- 
tion. The subject treatment is rather uniform in difficulty, 
and the text would probably be sufficiently challenging 
for students at these levels. However, since acoustics 
courses in this country are generally giverf to third- or 
fourth-year college students, the book is probably at too 
low a level to be practicable as a text. 

Among the better features of this book are the rather 
large number of worked problems, as well as a good 
collection of meaningful problems for homework assign- 
ments. 

One may also note what appears to be a standard 
British textbook practice, namely, the description of a 
large number of interesting experiments which could be 
performed in the laboratory or at the lecture bench. It is 
this reviewer’s opinion that if some of these experiments 
are actually performed, the recital can be of value, but 
that otherwise they tend to appear as so much superfluous 
material. 
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In several places, the author has omitted important 
modern developments. Thus the chapters on ultrasonics do 
not contain any account of modern pulse techniques, nor 
of the use of barium titanate elements as transducers. 
Also, the discussion of sound refraction could well use some 
mention of the use of plastic lenses. In the discussion on 
Kundt’s tube (p. 151) it would be appropriate to discuss 
measurements of velocity of compressional waves in rods 
by applying dust powder directly to the rod, as was done 
by Réhrich and many subsequent observers. 

It is the reviewer's task to point out that there exist a 
number of incorrect or misleading statements on ultrasonics 
and wave motion. The statement that CO, becomes 
‘“‘opaque’’ to sound waves at higher frequencies (p. 45) is 
certainly misleading, while the present upper limit of high 
frequency sounds is 500 times higher than the value of 
two megacycles given on p. 35. Finally, few people in 
underwater sound work will believe that “‘sound ranging 
in the sea also has the advantages that the sea is homogene- 
ous and there is no need for a meteorological correction” 
(p. 50). 

Because of statements such as these, it is difficult for 
the reviewer to recommend this book with much enthusi- 
asm. It was not, we believe, designed for the American 
market. 

RoBERT T. BEYER 
Brown University 


Comets and Meteor Streams. J. G. Porter. Pp. 123+-xi. 
Figs. 31, 16X24 cm. John Wiley and Sons, Inc., 
New York, 1952. Price $5.25. 


This small book is No. 2 in the International Astro- 
physics Series, edited by M. A. Ellison, The Royal Ob- 
servatory, Edinburgh, and A. C. B. Lovell, Professor of 
Radio Astronomy, University of Manchester. The editor’s 
note states that the aim of this series ‘‘is to provide a 
collection of authoritative volumes dealing with the main 
branches of Astrophysics and Radio Astronomy.” The 
reviewer's duty is to appraise the contents of the book on 
the basis of this clearly expressed intention. 

The book from this point of view is a dismal failure. 
It is nothing but a very short, even if excellently written, 
monograph on the orbital theory of comets and meteors. 
Of astrophysics, as this term is generally understood, the 
author is quite innocent. Spectra of comets are treated in 
less than a page, spectra of meteors are not even mentioned. 
The brightness of comets deserves somewhat more atten- 
tion—two whole pages are devoted to it; the brightness of 
meteors is again completely overlooked. The internal 
structure of the comet is simply mentioned in a half-page 
discussion of the problem of light pressure. One wonders 
what the author and editors mean by astrophysics. 

To any person familiar with the subject it is clear that 
comets are intimately connected with meteor streams, 
meteorites, and diffused matter in space. Precise observa- 
tion of all these phenomena is difficult, owing mostly to the 
transitory nature of most of these objects. We have, 
therefore, a large collection of observational data, often 
mutually exclusive and contradictory. Nevertheless, con- 
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siderable progress has been made in recent years in the 
collection of better data and theoretical interpretation of 
the spectra and luminosity of comets. Practically nothing 
of this is mentioned in the reviewed book. So far as meteors 
are concerned, the new and very promising method of 
radar observation is discussed in considerable detail (pp. 
72-74 and 99-103), but the contribution here is again 
toward the dynamical rather than physical properties of 
meteors. 

The quarrel is with the title and purpose of the book. 
If the author had called it ‘Orbits of Comets and Meteor 
Streams,” he might have been congratulated in carrying 
out his task in so short a space. There are many useful 
tables in this book such as a list of short period comets 
observed at more than apparition (p. 37), comet groups 
(pp. 48-49), cometary meteor radiants (p. 92), etc. The 
subject is extremely condensed, and it is doubtful that a 
person without a previous study of the theory of perturba- 
tions will derive much of an idea of this involved subject 
from the few pages devoted to it (pp. 54-65). The last 
14 pages of the book are occupied with instructions for the 
calculation of a comet’s ephemeris and of meteor orbits. 
If we subtract about eight pages of bibliography appended 
to various chapters, we find that the whole story of comets 
and meteors is told in about 96 pages. This much space 
would be necessary to treat fully and coherently the subject 
of the spectra of comets alone. 

There are six very good plates showing comets and 
meteors entirely unconnected with the text. 

The price $5.25 for such a small book seems to be out 
of all reason. How can a student or even a research worker 
be expected to invest so much in a book dealing with the 
restricted aspect of a very restricted problem? On the 
basis of pages (text, illustrations, formulas, etc., being 
comparable) the recent book Astrophysics, 703 pages 
(editor J. A. Hynek) published by the McGraw-Hill Co. 
should cost nearly $30.00. In reality its price is only $12.00. 

N. T. BoBROVNIKOFF 
Perkins Observatory 


A University Text-Book of Physics, Volume III, Heat. 
11th Edition. J. H. Awpery. Pp. 439+-x, Figs. 197, 
16X23 cm. Charles Griffin & Company, Ltd., London, 
1952. Price $5.50. 


Many physicists undoubtedly enjoyed reading at least 
one of the volumes of Poynting and Thomson’s well-known 
Text-Book of Physics at some time or other. This was an 
elementary but exact presentation of physics by distin- 
guished physicists. But textbooks, unlike creative ideas, 
rarely outlive their authors. This, I am afraid is true of 
Poynting and Thomson’s ‘‘Heat,’’ now partly rewritten 
by Mr. Awbery, and available in a new edition, printed 
and bound in the best English tradition. 

While modernized in many respects, the book is certainly 
not up-to-date. With the exception of the references in 
some sections added by Mr. Awbery, most references are 
to old sources and texts no longer easily available. This 
is in good order when important original work is referred to, 
but little purpose is served when several methods of meas- 
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uring a single quantity are described, some only differing 
in detail, and every one of them antiquated. Thus, seven 
methods of measuring the solar constant are outlined, yet 
none later than 1910. The chapter on water in the atmos- 
phere (a field where progress has been particularly rapid 
since the advent of aviation, sondes and radio) similarly 
summarizes the knowledge and experimental equipment of 
the 1900’s and includes the canard (dated 1905), now 
famous among geophysicists, that the masses of individual 
raindrops in the same rainfall are in the ratio 1:2:4:8. 

Considerable space is devoted to thermodynamics. The 
first law is presented as the result of a number of experi- 
ments on the ‘‘mechanical equivalent of heat’’; altogether 
too much attention is devoted to a discussion of J, which 
is logically nothing but a measure of the heat capacity of 
water. The discussion of the second law is good and includes 
the usual topics—Clausius’ inequality, entropy, and defini- 
tion of temperature. There is also a chapter on the thermo- 
dynamics of radiation. 

An excellent, but very elementary, account is given of 
the kinetic theory of gases, including transport phenomena; 
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this is supplemented at the end of the book by a chapter 
on statistical mechanics and quantum theory. Even though 
the more difficult mathematical details are avoided, it is 
doubtful whether a student who turns to this book for 
instruction in the elementary topics can benefit much 
from it. More advanced readers, on the other hand, will 
more likely than not be put off by the elementary character 
of the remainder of the book. 

This, as some other features, suggest that this book will 
probably not be found suitable as a text for any one course 
on heat, though it might serve as a companion for the 
student from his first contact with the subject of heat 
through to his graduate work. The student would miss, 
however, a more adequate coverage of recent develop- 
ments and a treatment of the more philosophical parts of 
the subject (such as the relation of experimental heat to 
thermodynamics to kinetic theory). He might also deplore 
the absence of problems. 

WALTER HITSCHFELD 
McGill University 
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Illinois Section 


The fall meeting of the Section was held at Rockford 
College, Rockford, Illinois, on Friday and Saturday, 
November 7 and 8, 1952, with an attendance of approxi- 
mately 40. A varied program, organized by the local 
chairman, PROFESSOR A. FRANCES JOHNSON, was presented. 


1. Preparation of specimens for use in an electron 
microscope. MILDRED BULLIET, Rockford College—A dem- 
monstration was made by allowing a drop of a solution of 
parlodion in amyl acetate to spread into a thin film over 
the surface of water in a wide dish. Small round screens of 
fine mesh which fit into the specimen holder of the micro- 
scope were placed on the film. The film and screens were 
then picked up on a cover glass and allowed to dry. A 
drop of a suspension in water of particulate matter to be 
viewed in the microscope was placed on each screen. 
When dry the specimen was ready to be used. A variation 
of this method was made in showing how the small particles 
might be dispersed in a parlodion solution before it was 
dropped on the water. 

The preparation of Formvar and silica replicas of sur- 
faces was described. Pictures of particles and replicas made 
at Rockford College were shown. 


2. Two new physics films, circular motion, Carnot cycle 
and Kelvin scale. 


3. Modern applications of an old theory. LEONARD R. 
INGERSOLL, Emeritus Professor of Physics, University of 
Wisconsin.—While heat-conduction theory is, by and large, 


a century old, it has many very modern and up-to-date 
applications. These include thermal insulation in all its 
various aspects, with special graphical solutions for cases 
like corners and edges. Subterranean temperatures can be 
treated by Fourier analysis, and this has been applied to 
a study of the “‘thermal history”’ (postglacial time calcula- 
tions) of northern Michigan on the basis of special tempera- 
ture measurements in the deep Calumet and Hecla copper 
mines. 

One of the most interesting and recent applications is 
in connection with the heat pump. When this rapidly 
growing method of year-round air conditioning uses a 
ground coil as its source (or sink) of heat, line-source 
conduction theory can be used to good advantage, since 
it is the only practical theoretical way of avoiding the 
various pitfalls which have so far retarded progress with 
this phase of heat-pump design. P 


4. Adding two distributions. Epwarp C. VARNUM, 
Barber-Colman Company.—By adding random pairs se- 
lected from two independent distributions, a third distribu- 
tion is formed having a mean and variance equal, respec- 
tively, to the sum of the means and variances of the 
original distributions. Assuming that the distribution of 
measuring-instrument errors is independent of the distri- 
bution of actual data in a physical test, the above property 
may be applied to the study of an observed data distribu- 
tion. When the differences of the random pairs are taken 
in a specified order, the variances are also additive. This 
result has been applied to the study of overlapping toler- 
ances of mating parts. A chart for determining the proba- 
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bility of interference of two overlapping normal distribu- 
tions was given to each member of the audience, and the 
use of the chart was illustrated by means of a numerical 
example expressed in terms of the distributions of shaft 
and hole diameters in an assembly process. 


5. Impressions of great German scientists as recalled 
by a student of physics in Munich in the early 1930's. 
MARGARETE BrucH, Winnebago.—The speaker, daughter 
of the late professor of chemistry, Dr. Richard Willstatter, 
had the opportunity both at home and as a student at the 
University in Munich to know Professor Wilhelm Conrad 
Réntgen, tall and impressive, with a soft voice and tiny 
handwriting; his successor Wilhelm Wien, known for his 
work on blackbody radiation; and Wien’s successor, 
Walter Gerlach. She also had the good fortune to study 
under Arnold Sommerfeld, well known for his explanation 
of the fine structure of spectra and one of the first to 
accept Bohr’s theory. An inspiring teacher, he developed, 
in his lectures before his students, his theory of the electron 
gas obeying Fermi statistics. She heard beautifully clear 
lectures by Max Planck and challenging lectures by Albert 
Einstein, whose interest in the most varied problems was 
shown in his talking with the same verve on unified field 
theory and on the forces that make moist sand cling 
together. (This talk was illustrated by portraits and 
autographed letters.) 


6. An operational approach for the classroom develop- 
ment of basic concepts and Newton’s laws. Cari J. 
RIGNEY, Northern Illinois State Teachers’ College-—In this 
instructional method, care is taken to employ an explicitly 
outlined scientific procedure which involves operational 
definitions only. The procedure of physics is formulated as 
“the development of our concepts of what is basic in 
nature by a systematic use of what is basic in our experience 
with the physical world.’’ Length, time, and force are 
agreed on as mechanical concepts that are basic in our 
experience. Weights are used to calibrate a spring scale in 
arbitrary force units. The findings of an imaginary scien- 
tific expedition show that the weight of an object at any 
location is directly proportional to the free-fall acceleration 


RECENT MEETINGS 


at that location. The constant of proportionality for a 
chosen object is selected arbitrarily as the unit mass, and 
the force unit—no longer arbitrary—follows. 


7. What can we do for the less-than-mediocre students 
now crowding our colleges? CHALMER N. PATTERSON, 
Bradley University—The seriousness of this problem is 
widely acknowledged. Quick elimination of such students 
is rejected because (a) many weak students have, by 
properly guided hard work, achieved considerable success; 
(b) even failing students cam learn much of value to them- 
selves when theories are properly presented with careful 
illustration from daily life. 

Three possible solutions are suggested: (1) Easy in- 
formational courses, credited toward graduation at only 
half value, with D grade, and only after the student has 
earned the required grade point average in other classes. 
Objection: Such students are already labeled ‘‘Failure.” 
(2) Survey courses with one extra period, open to all but 
required of weak students, scheduled each week. Only D 
grade permitted for student required to attend extra 
period unless he shows notable improvement as semester 
progresses. (3) Similar extra periods for students entering 
courses for Arts students with minimum prerequisites. 

The writer recognizes that there may be many other 
ways of trying to improve a bad situation. 

Group activities included a demonstration of an electron 
microscope in operation, a visit to the Air Distribution 
Development Laboratory of the Barber-Colman Company, 
a discussion and demonstration of a Philbrick analog 
computer, and a tour of the laboratories of the Woodward 
Governor Company. At the business meeting of the 
Section, Proressor O. L. RatLtsBack, University of 
Illinois, Navy Pier, Chicago, reported upon recent national 
meetings, and PROFESSOR GLEN Q. LEFLER, Eastern 
Illinois State College, presented a report of the Committee 
on Constitution Revision. A dinner, a luncheon, and a 
special program for the ladies offered pleasant occasions 
for social contacts. 

O. L. RAILSBACK 
A. FRANCES JOHNSON 


Symposium on Molecular Structure 


The Symposium on Molecular Structure and Spectroscopy will be held at the Department 
of Physics and Astronomy, The Ohio State University, as usual in June, from June 15 to June 
19, 1953. There will be discussions of the interpretation of molecular spectroscopic data as well 
as methods for obtaining such data. In addition, there will be sessions devoted to those phases 
of spectroscopy of current interest. A dormitory will be available for those who wish to reside 
on the campus during the meeting. For further information, or for a copy of the program 
when it becomes available, write to Professor Robert A. Oetjen, Department of Physics and 
Astronomy, The Ohio State University, Columbus 10, Ohio. 
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Important New Books for Classroom Use 


MICROWAVE THEORY 
AND TECHNIQUES 


By H. J. REICH, P. F. ORDUNG, 
H. L. KRAUSS, and J. G. SKALNIK 


The comprehensive treatment of the theory and 
applications in this new text is carried from field 
theory through transmission lines, waveguides, 
radiation, and microwave amplifiers and oscilla- 
tors. Thus, the textual material plus the gener- 
ous references to papers and books, make this a 
comprehensive text and a valuable reference book. 
Ready in June. 


PRINCIPLES 
OF MODERN ACOUSTICS 


By GEORGE W. SWENSON, Jr. 


This new text develops a number of important 
analytical techniques. It will appeal to the scien- 
tist because of its thorough explanation and use 
of equivalent circuits for mechanical and acousti- 
cal systems. It concentrates almost entirely on 
methods and theory, rather than details of the 
myriad of modern applications and devices. It 
is spiced with a number of example applications. 
Ready in April. 


ELECTRONICS 


PHYSICAL PRINCIPLES & APPLICATIONS 
By A. O. WILLIAMS, Jr. 


This new text gives a rounded, physical introduc- 
tion to the subject, without exhaustive enumera- 
tion of special circuits, but with constant attention 
to the associated phenomena of electricity, mag- 
netism, radiation (including optics), and atomic 
physics that make electronics an interesting and 
useful tool to the physicist. Especially suited for 
physics, chemistry and other majors. 306 pages. 
$4.50 


GENERAL 
COLLEGE PHYSICS 
By ERTLE LESLIE HARRINGTON} 


The outstanding feature of this text is the selec- 
tion of illustrations and applications. Here the 
student will find, repeatedly, applications of 
physics with which he has become acquainted but 
without recognizing them as such nor even under- 
standing the physical principles involved. Author 
took great care to see that his text would support 
the standard laboratory course. 787 pages. $6.75 


REVISED 2nd EDITION 


MODERN COLLEGE PHYSICS 


By HARVEY E. WHITE 


THE OUTSTANDING CHANGES IN THIS NEW TEXT are the considerable expansion of alternating 
current theory and the use of the mks system of units in electricity and magnetism. Chapters on éosmic 
rays, nuclear physics, and atomic energy have been largely rewritten to bring them up to date. All new 
problems are given at the end of each chapter. 823 pages. $6.75 


EXAMINATION COPIES SENT ON APPROVAL 


D. VAN NOSTRAND COMPANY, INC. 


250 FOURTH AVENUE 


PUBLISHERS SINCE 1848 





NEW YORK 3, NEW YORK 


CONTENTS 


AMERICAN JOURNAL OF PHysICcs VoLuME 21, NuMBER 4 APRIL, 


Construction of a Diffraction Grating Spectrograph 
On a Routine Analytic Method for the Solution of Problems in Statics A. W. Sumon 
Constant Battery Wheatstone Bridge T. TowNsENpD SMITH 


Spin Waves in Ferromagnetic and Antiferromagnetic Materials 
F. Kerrer, H. KAPLAN, AND: Y. YAFET 


H. Morz anp L. I. Scuirr 
Behavior of Magnetic Materials R. M. Bozort# 


Rocking Experiment with Two Degrees of Freedom Joun SATTERLY 
Geometrical Solution for a Nuclear Moment in a Magnetic Field . . Roatp K. WANGSNEsSS 
Physics and Physiology in Diving Decompression HILBERT SCHENCK, JR. 


Calculus of Measures and Its Application to the Relation Between Rationalized and Unration- 
alized Electromagnetic Systems J. R. BARKER AND R. O. Davies 


Longitudinal Vibrations of a Vertical Column by the Method of Laplace Transform . L. H. HALL 
Proposed Procedure for Selecting and Using Symbols for Physical Units . . DUANE ROLLER 


Teaching a Philosophy of Experimentation in a Course in Electrical Measurements ~ 
Francis T. WORRELL 

Experimental Reactor Physics Course at the Oak Ridge School of Reactor Technology 
E. D. Kiema, R. J. STEPHENSON, AND S. TAYLOR 
Surface Energy, A Mode for Energy Absorption During Impact Joun S. RINEHART 


Undergraduate Experiment for the Determination of Alpha-Particle Range 
WALTER C. MICHELS AND Eva WIENER 
Relativistic Rocket Theory 


Notes and Discussion: 
Effect of Absorption by the Material of the Prism on Its Resolving Power . MAHENDRA SINGH SODHA 
Determination of the Coefficient of Surface Tension by the Bubble-Length Method 
Joun C. SLONCZEWSKI AND RALPH HELLER 
Saut ARonow 


Letters to the Editor: 


Extension of the Conical Pendulum Problem and its Demonstration . . . Juttus SUMNER MILLER 
Ultrasonic Interferometer as a Screw Calibrating Instrument RomMARD BARTHEL 
Concerning the Electric Charge of a Moving Vehicle Juttus SUMNER MILLER 
Demonstration of Subjective Harmonics CHARLES WILLIAMSON 

W. A. MurtTAuGH 


Announcements and News: 


Practical Aids for Physics Teachers 
Tape Recordings of Important Speeches 
To Secretaries of Regional Sections of AAPT 


Symposium on Molecular Structure 
Recent Meetings: 





